bzoj2818 欧拉函数

题目大意:求x,y<=n的gcd(x,y)为素数的共有多少对
思路:不妨简化这个问题,gcd(x,y)=1的共有多少对.
即求对于y的欧拉函数值.设(x<=y)
解决了这个问题之后得到:gcd(px,py)=p(p为素数)
那么y<=n/p;
所以我们枚举1~n的每一个素数,计算它对于当前答案的贡献值.(与简化问题相同)
注意:gcd(x,y)=gcd(y,x) 且当x=y=1时只算一遍

#include<iostream>
#include<cstdio>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=1e7+10;
typedef long long ll;
ll prime[N],fai[N],tot,n,sum[N],ans;
bool check[N];
void Euler_prime()
{
    memset(check,0,sizeof(check));
    fai[1]=1; tot=ans=0;
    for(int i=2;i<=n;i++) {
        if(!check[i]) {
            prime[++tot]=i;
            fai[i]=i-1;
        }
        for(int j=1;j<=tot;j++) {
            if(i*prime[j]>n) break;
            check[i*prime[j]]=1;
            if(i%prime[j]==0) {
                fai[i*prime[j]]=fai[i]*prime[j];
                break;
            }else {
                fai[i*prime[j]]=fai[i]*(prime[j]-1);
            }
        }
    }
}
int main()
{
//  freopen("prime.in","r",stdin);
//  freopen("prime.out","w",stdout);
    scanf("%lld",&n); Euler_prime();
    for(int i=1;i<=n;i++)
        sum[i]=sum[i-1]+fai[i];
    for(int i=1;i<=tot;i++)
        ans+=sum[n/prime[i]]*2-1;
    printf("%lld",ans);
}