BZOJ2818 Gcd

轉換成 即可。

先算出歐拉函數，再算出互質對的前綴和pre[n]。

注意除了a=1&&b=1的情況，其它的a！=b算是兩對。（gcd(a,b)與gcd(b,a)）

gcd(a,b)=d的個數與gcd(n/d,n/d)=1個數是一樣的。

注意爆int，要用longlong。

#include <bits/stdc++.h>
#define ms(a,b) memset(a,b,sizeof(a))
#define mc(a,b) memcpy(a,b,sizeof(b)) 
using namespace std;
typedef long long ll;
const int N=1e7+5;
const ll mod=1e7+7;
int pri[N],vis[N],phi[N];
ll pre[N];
void prime(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(!vis[i]) 
        {
            pri[++pri[0]]=i;
            phi[i]=i-1;
        }
        for(int j=1;j<=pri[0]&&pri[j]*i<=n;j++)
        {
            vis[pri[j]*i]=1;
            if(i%pri[j]==0)
            {
                phi[pri[j]*i]=pri[j]*phi[i];
                break;
            }
            else phi[pri[j]*i]=(pri[j]-1)*phi[i];
        }
    }
}
int main()
{
    int n;
    scanf("%d",&n);
    prime(n);
    phi[1]=1;
    ll ans=0;
    pre[1]=1;
    for(int i=2;i<=n;i++) pre[i]=2*phi[i]+pre[i-1];
    for(int i=1;i<=pri[0];i++) ans+=pre[n/pri[i]];
    printf("%lld",ans);
    //system("pause");
}