P2398 GCD SUM 欧拉函数

qq_38232157

已于 2022-10-13 14:54:16 修改

阅读量165

点赞数

分类专栏：洛谷数论文章标签：数论

于 2022-10-13 14:50:33 首次发布

本文链接：https://blog.csdn.net/qq_38232157/article/details/127301403

版权

1、先打印素数表 prime[i] ，在计算欧拉函数 phi[i]
2、gcd(x, y), 1 <= x <= n, 1 <= y <= n, 题目的意思相当于有 n * n 个 gcd(x, y)
求和。假如 x < y,且 gcd(x, y) = 1, 则小于 y 并且与 y 互素的数有 phi[y]
个。这 phi[y] 对数 (x, y) 满足： gcd(k * x, k * y) = k， k <= n / phi[y]即可。
其相关的 gcd(x, y)总和 S
= phi[y] * (n / phi[y] + n / phi[y] - 1 + n / phi[y] - 2 + … + 2 + 1)
3、1 ~ n 这n个数，按照上面式子算出 S，在累加得到 sum(S)。因为 gcd(x, y) = gcd(y, x)
, 结果需要乘以2， sum(S) * 2
4、两数相等时，所有的 gcd(x, x)总和为 (1 + n) * n / 2.
最后的结果： sum(S) * 2 + (1 + n) * n / 2. 题目需要注意 int 溢出的情况

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
using namespace std;
const int MaxN = 100010;
int phi[MaxN];	//欧拉函数
int prime[MaxN];
int m;