BZOJ 2818 GCD

思路：因为要求GCD（x,y)=m，m为素数，有多少有序对。则可以转化为求GCD（a,b)=1，x=a*m，y=b*m，有多少有序对。假设x<=y,则a<=b，因为1<=x<=y<=n，则b<=(n/m).通过枚举素数m，求小于b，与b互斥的有多少个。通过处理欧拉函数的前缀和，每次会多求出一个a=b=1，GCD（a,b)=1的解，又因为是有序对，所以结果为sum[b]*2-1.

#include<bits/stdc++.h>
using namespace std;
typedef long long int ll;
const int MAX=1e7;
int tot,n,pri[MAX+10],phi[MAX+10];
ll sum[MAX+10];
bool mark[MAX];
void getphi()
{
    phi[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!mark[i])
        {
            phi[i]=i-1;
            pri[++tot]=i;
        }
        for(int j=1;j<=tot;j++)
        {
            if(i*pri[j]>n)break;
            mark[i*pri[j]]=1;
            if(i%pri[j]==0){phi[i*pri[j]]=phi[i]*pri[j];break;}
            else phi[i*x]=phi[i]*phi[pri[j]];
        }
    }
}
int main()
{
    tot=0;
    ll ans=0;
    while(~scanf("%d",&n))
    {
        ans=0;
        getphi();
        for(int i=1;i<=n;i++)sum[i]=sum[i-1]+phi[i];
        for(int i=1;i<=tot;i++)ans+=sum[n/pri[i]]*2-1;
        printf("%lld\n",ans);
    }
    return 0;
}