[2021CCPC威海热身赛] Number Theory (打表 推公式)
题意
求
∑ k = 1 n ∑ i ∣ k ∑ j ∣ i λ ( i ) λ ( j ) \sum_{k = 1}^{n}\sum_{i \mid k} \sum_{j \mid i} \lambda(i) \lambda(j) k=1∑ni∣k∑j∣i∑λ(i)λ(j)
对 998244353 998244353 998244353 取模
其中 λ ( x ) = ( − 1 ) ∑ i e i , x = ∏ i p i e i \lambda(x) = (-1)^{\sum\limits_{i}e_i},x=\prod\limits_{i}p_i^{e_i} λ(x)=(−1)i∑ei,x=i∏piei
分析:
λ ( x ) \lambda(x) λ(x) 为刘维尔函数,可以打表发现 ∑ d ∣ n λ ( d ) = [ n = a 2 , a ∈ N + ] \sum_{d \mid n}\lambda(d) =[n = a^2,a \in N^+] d∣n∑λ(d)=[n=a2,a∈N+]
也就是 n n n 是否为完全平方数
把式子中的 λ ( i ) \lambda(i) λ(i) 提到前面
∑ k = 1 n ∑ i ∣ k λ ( i ) ∑ j ∣ i λ ( j ) \sum_{k = 1}^{n}\sum_{i \mid k} \lambda(i)\sum_{j \mid i} \lambda(j) k=1∑ni∣k∑λ(i)j∣i∑λ(j)
那么就变为
∑ k = 1 n ∑ i ∣ k λ ( i ) [ i = a 2 , a ∈ N + ] \sum_{k = 1}^{n}\sum_{i \mid k} \lambda(i)[i= a^2,a \in N^+] k=1∑ni∣k∑λ(i)[i=a2,a∈N+]
那么完全平方数的刘维尔函数为 1 1 1,再设 f ( x ) = [ i = a 2 , a ∈ N + ] f(x)=[i= a^2,a \in N^+] f(x)=[i=a2,a∈N+] 得
∑ i = 1 n ∑ d ∣ i f ( d ) \sum_{i = 1}^{n}\sum_{d \mid i}f(d) i=1∑nd∣i∑f(d)
交换求和次序
∑ d = 1 n f ( d ) ⌊ n d ⌋ \sum_{d = 1}^{n}f(d) \lfloor\frac{n}{d}\rfloor d=1∑nf(d)⌊dn⌋
这样直接枚举平方数即可,时间复杂度 O ( n ) O(\sqrt{n}) O(n )
代码:
#include <bits/stdc++.h>
#define int long long
using namespace std;
const int mod = 998244353;
int n, res;
signed main() {cin >> n;for (int i = 1; i * i <= n; i ++) {res = (res + n / (i * i)) % mod;}cout << res << endl;
}
[2021CCPC威海热身赛] Number Theory (打表 推公式)相关推荐
- hdu 4215 Number Theory?(打表)
题目: E(N) = |{i | gcd(N, i) = 1, 1 <= i <= N}| F(N) = |{i | N % i = 0, 1 <= i <= N}| 求有多少 ...
- 2016级算法第二次上机-F.ModricWang's Number Theory II
891 ModricWang's Number Theory II 思路 使得序列的最大公约数不为1,就是大于等于2,就是找到一个大于等于2的数,它能够整除序列中的所有数. 考虑使得一个数d整除数组中 ...
- concrete maths ch4 number theory
ch4 number theory 数论研究正数的性质 1.整除 gcd lcm 扩展欧几里得. 整除求和\(\sum_{n|m}\)的几个公式.ch2的知识会很有用. 2.质数 Fundamenta ...
- Number Theory Problem(The 2016 ACM-ICPC Asia China-Final Contest 找规律)
题目: Mr. Panda is one of the top specialists on number theory all over the world. Now Mr. Panda is in ...
- NUMTRYE - Number Theory (Easy)
NUMTRYE - Number Theory (Easy) Hard 版本就是用 pollard_rho 分解质因子. f(n)=∏(pi2ei+1+1)f(n) = \prod(p_i ^{2e_ ...
- Codeforces Round #628 (Div. 2) E. Ehab‘s REAL Number Theory Problem 巧妙的质因子建图
传送门 文章目录 题意: 思路: 题意: 给你nnn个数,每个数的因子个数不超过777个,选出最少的数使其乘积为平方数. n≤1e5n\le 1e5n≤1e5 思路: 由于因子不超过777个,所以由约 ...
- Segment 2:Introduction Number Theory——Fermat and Euler【费马定理和欧拉定理】
Segment 2:Introduction Number Theory--Fermat and Euler[费马定理和欧拉定理] 这是整个数论简介的内容,下面是这个的主目录,其链接为:https:/ ...
- Elliptic Curves Number Theory And Cryptography——pairing learning
<Elliptic Curves Number Theory And Cryptography 2n>中Example 11.5, magma脚本: clear; q:=7; Fq:=GF ...
- CodeForces - 466C Number of Ways(推公式/dp)
题目链接:点击查看 题目大意:给出一个长度为 n 的数列,现在要求出满足条件的 ( i , j ) 的匹配数量,满足: 题目分析:训练时推的公式,简单说一下吧,维护前缀和 sum,则确定两个断点 ( ...
最新文章
- gis 数据框裁剪_BIM+GIS的八大挑战!大挑战,见未来
- 火币矿池参选EOS超级节点, 不只为分一杯羹
- 多分辨率适配常用目录
- TiDB 源码阅读系列文章(十五)Sort Merge Join
- 《卓有成效的程序员》----读书笔记二
- 什么是透传通道?(透明传输通道,就是当中继使)
- python僵尸进程和孤儿进程_python中多进程应用及僵尸进程、孤儿进程
- logo自动旋转 html,HTML/CSS3 Logo选择器
- python mysqldb 查询不到最新记录_python – MySQLdb是否缓存SELECT结果?
- 《linux c编程指南》学习手记2
- GoLang笔记—容器篇
- Unity资源加载闪退问题深度分析
- 各大云服务厂商 轻量应用服务器 性能评测对比,阿里云、腾讯云、华为云、Ucloud
- 刚上大一,C 语言压根听不懂怎么办?
- 数据类型,栈内存、堆内存
- 2023五一旅游必备物品清单!快记到手机待办APP里
- 2022计算机软件毕业设计大全
- html 设置请求头,http请求头和响应头设置
- 限时免费 | 学习成本超低的Python数据分析手册
- 英文Essay写作诀窍讲解