AcWing220-最大公约数

给定整数 $N$ ，求 $1\le x,y\le N$ 且 $GCD(x,y)$ 为素数的数对 $(x,y)$ 有多少个

考虑 $1$~ $N$ 内的每个质数 $p_i$, 都有 $$\gcd (x,y)=p_i$$
对该式化简有 $$\gcd (\frac{x}{p_i},\frac{y}{p_i})=1$$
不妨设 $1\le x<y\le N$，而 $y<x$ 的情况仅需要乘 $2$ 即可
那么对每个 $p_i$ 有 $$1\le \frac{x}{p_i}<\frac{y}{p_i}\le \frac{N}{p_i}$$
也就是说，对每个 $\frac{y}{p_i}$ 我们要求与他互质的个数，即对每个 $p_i$ 求 $$\sum _{j=1}^{\lfloor \frac{N}{p_i}\rfloor }\phi (j)$$
再去想当 $j=1$ 的时候,也就是 $y=p_i$，此时也必然有 $x=p_i$，而两个相等的数只能算一种，因此最后答案$$\sum _{i=1}^{cntp}(2\times \sum _{j=1}^{\lfloor \frac{N}{p_i}\rfloor }\phi (j)-1)$$
对于后一个和式可用前缀和在 $O(n)$ 内算出

CODE

#include <iostream>
#define ll long long

using namespace std;

const int maxn = 1e7 + 10;

int npri[maxn];
int pri[maxn], eurl[maxn];
ll sum[maxn];
int cntp = 0;


void geteurl(int n){
    eurl[1] = 1;
    for (int i = 2; i <= n; i++){
        if (!npri[i]){
            pri[++cntp] = i;
            eurl[i] = i - 1;
        }
        for (int j = 1; i * pri[j] <= n; j++){
            npri[pri[j] * i] = 1;
            if (i % pri[j] == 0){
                eurl[i * pri[j]] = eurl[i] * pri[j];
                break;
            }
            eurl[i * pri[j]] = eurl[i] * (pri[j] - 1);
        }
    }
}


int main(){
    geteurl(maxn);
    for (int i = 1; i <= maxn; i++) sum[i] = sum[i - 1] + eurl[i];
    int n;
    cin >> n;
    ll ssum = 0;
    for (int i = 1; i <= cntp; i++){
        if (pri[i] > n) break; 
        ssum += 2 * sum[n / pri[i]] - 1;
    }
    cout << ssum;
    return 0;
}