AcWing 220. 最大公约数 （欧拉函数）

AcWing 220. 最大公约数

题意

给一个整数 $n(1\le n\le 10^7)$ ，问 $1\le x,y\le n$ 且 $gcd(x,y)$ 是素数的数有多少对。

解法

$$\begin{array}{rl}ans& =\sum _p\sum _{x=1}^n\sum _{y=1}^n[gcd(x,y)==p]\\ & =\sum _p\sum _{x=1}^n\sum _{y=1}^n[gcd(x/p,y/p)==1]\\ & =\sum _p\sum _{x=1}^{n/p}\sum _{y=1}^{n/p}[gcd(x,y)==1]\\ & =\sum _p(\sum _{x=1}^{n/p}2(\sum _{y=1}^x[gcd(x,y)==1])-1)\\ & =\sum _p(\sum _{x=1}^{n/p}2\phi (x)-1)\end{array}$$

筛出欧拉函数和 $n$ 以内的质数，处理出前缀和，然后暴力加即可。

代码

#pragma region
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <map>
#include <queue>
#include <set>
#include <vector>
using namespace std;
typedef long long ll;
#define tr t[root]
#define lson t[root << 1]
#define rson t[root << 1 | 1]
#define rep(i, a, n) for (int i = a; i <= n; ++i)
#define per(i, a, n) for (int i = n; i >= a; --i)
#pragma endregion
const int maxn = 1e7 + 5;

int cnt = 0, prime[maxn];
bool mp[maxn];
void phi_table(int n, int* phi) {
    for (int i = 2; i <= n; i++) phi[i] = 0;
    phi[1] = 1;
    for (int i = 2; i <= n; i++) {
        if (!phi[i]) {
            prime[++cnt] = i;
            for (int j = i; j <= n; j += i) {
                if (!phi[j]) phi[j] = j;
                phi[j] = phi[j] / i * (i - 1);
            }
        }
    }
}
int phi[maxn];
ll pre[maxn];
int main() {
    int n;
    scanf("%d", &n);
    phi_table(n, phi);
    rep(i, 1, n / 2) pre[i] = pre[i - 1] + 2 * phi[i];
    ll ans = 0;
    rep(i, 1, cnt) {
        int p = prime[i];
        if (p > n) break;
        ans += pre[n / p] - 1;
    }
    printf("%lld\n", ans);
}