转换成 即可。
先算出欧拉函数,再算出互质对的前缀和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");
}