F(n)=∑ni=1∑nj=1ϕ(ϕi,ϕj)
其中
ϕ
表示欧拉函数。欧拉函数
ϕn
是不超过n的数中与n互质的数的数目。
ϕ(ϕi,ϕj)
表示i,j欧拉函数值的最大公约数的欧拉函数值.
给出n,求F(n)的值。
Input
一行,包含正整数T,表示数据组数。T<=5 接下来T行每行一个正整数n。n<=2*10^6
Output
共T行,每行包含一个答案。
Input示例
1 1
Output示例
1
莫比乌斯反演~
推荐讲解:http://blog.csdn.net/worldwide_d/article/details/53427483~
不能先预处理再更新答案;记得清空数组!
#include<cstdio>
#include<cstring>
#define N 2000001
#define ll long long
int t,n,m,a[N>>3],phi[N],mu[N];
ll s[N],f[N],ans;
bool b[N];
void sol()
{
ans=0;scanf("%d",&n);
mu[1]=phi[1]=1;a[0]=0;
memset(f,0,sizeof(f));
memset(s,0,sizeof(s));
for(int i=2;i<=n;i++)
{
if(!b[i]) a[++a[0]]=i,mu[i]=-1,phi[i]=i-1;
for(int j=1;j<=a[0] && i*a[j]<=n;j++)
{
b[i*a[j]]=1;
if(!(i%a[j]))
{
phi[i*a[j]]=phi[i]*a[j];break;
}
phi[i*a[j]]=phi[i]*(a[j]-1);mu[i*a[j]]=-mu[i];
}
}
for(int i=1;i<=n;i++) s[phi[i]]++;
for(int i=1;i<=n;i++)
for(int j=i;j<=n;j+=i) f[i]+=s[j];
for(int i=1;i<=n;i++) f[i]=f[i]*f[i];
for(int i=1;i<=n;i++) if(mu[i]) for(int j=1;i*j<=n;j++) ans+=mu[i]*phi[j]*f[i*j];
printf("%lld\n",ans);
}
int main()
{
scanf("%d",&t);
while(t--) sol();
return 0;
}