51nod 1594 Gcd and Phi

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;
}