P2568 线性筛 + 欧拉函数

题意

传送门 P2568 GCD

题解

$p$ 代表 $[1,n]$ 范围内的质数，求解 $$\sum _{1\le p\le n}\sum _{1\le x\le n}\sum _{1\le y\le n}[gcd(x,y)==p]$$ 即 $$\sum _{1\le p\le n}\sum _{1\le x\le n/p}\sum _{1\le y\le n/p}[gcd(x,y)==1]$$ 考虑到 $x,y$ 值域的对称性，$x\le y$ 限制下互质的数对个数与 $y\le x$ 限制下的一致，重复计数的情况只有 $(1,1)$。对于 $x\le y$ 的情况，可以枚举 $y$，那么对于固定的 $y$，互质的数对个数为 $\varphi (y)$，则在固定质数的情况下，问题转化为求解欧拉函数的前缀和。设欧拉函数 $[1,i]$ 的前缀和为 $su{m}_i$，则答案为 $$\sum _{1\le p\le n}(2\times sum[n/p]-1)$$

#include <bits/stdc++.h>
using namespace std;
#define maxn 10000005
typedef long long ll;
int n, np, prime[maxn], minp[maxn], phi[maxn];
ll sum[maxn];

void sieve(int n)
{
    phi[1] = 1;
    for (int i = 2; i <= n; ++i)
    {
        if (minp[i] == 0)
            prime[np++] = i, minp[i] = i, phi[i] = i - 1;
        for (int j = 0; j < np && i * prime[j] <= n; ++j)
        {
            int k = i * prime[j];
            minp[k] = prime[j];
            if (minp[i] == prime[j])
            {
                phi[k] = prime[j] * phi[i];
                break;
            }
            else
                phi[k] = (prime[j] - 1) * phi[i];
        }
    }
    for (int i = 1; i <= n; ++i)
        sum[i] = sum[i - 1] + phi[i];
}

int main()
{
    cin >> n;
    sieve(n);
    ll res = 0;
    for (int i = 0; i < np; ++i)
        res += (sum[n / prime[i]] << 1) - 1;
    cout << res << endl;
    return 0;
}