质数-luogu2568/bzoj2818GCD

给定整数N，求1<=x,y<=N且Gcd(x,y)为素数的
数对(x,y)有多少对.

因为gcd x,y为素数
所以(x/g,y/g)=1
假设x比y小
只需枚举小于n的素数g
用欧拉函数求x的个数
最后*2再加上x=y的情况

可以用前缀和优化

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define maxn 10000005
#define LL long long
using namespace std;
int n,phi[maxn],prime[maxn],tot;
LL ans,sum[maxn];
bool vis[maxn];

int main(){
    scanf("%d",&n);
    for(int i=2;i<=n;i++){
        if(!vis[i]) prime[++tot]=i,phi[i]=i-1;
        for(int j=1;j<=tot && i*prime[j]<=n;j++){
            vis[prime[j]*i]=1;
            if(i%prime[j]==0) {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }
    for(int i=1;i<=n;i++) sum[i]=sum[i-1]+phi[i];
    for(int i=1;i<=tot;i++){
        ans+=sum[n/prime[i]];
    }
    printf("%lld\n",ans*2+tot);
    return 0;
}