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;