2017.9.19 Gcd 思考记录

一开始想错了    

    有时候还是要用欧拉函数的。。

    用反演的话需要提取质数，比较麻烦，

先利用初中奥数分成gcd和两个质数，然后剩下的就是互质了

设i=gcd*a ,  j=gcd*b      , b>a所以

ans=   phi（2/2）+phi（4/2）+phi（6/2）+...phi（n/2）

     +     phi（3/3）+phi（6/3）+phi（9/3）+...phi（n/3）

...

前缀和即可

码：

#include<iostream>
#include<cstdio>
using namespace std;
#define N 10000007
long long su[N],phi[N],n,i,ans,tot,lin,j;
bool he[N];
void eular(int n)
{
	phi[1]=1;
	for(i=2;i<=n;i++)
	{
	if(!he[i])	
		{
			su[++tot]=i;
			phi[i]=i-1;			
		}
		for(j=1;j<=tot&&su[j]*i<=n;j++)
		{
			he[su[j]*i]=1;
			if(i%su[j]==0)
			{
				phi[i*su[j]]=phi[i]*su[j];
				break;
			}else phi[i*su[j]]=phi[su[j]]*phi[i];	
		}	
	}	
	for(i=2;i<=n;i++)phi[i]+=phi[i-1];
}
int main()
{
	scanf("%lld",&n);
	eular(n);
for(i=1;i<=tot;i++)
{//cout<<su[i];
	ans+=(phi[n/su[i]]*2-1);
}
	printf("%lld",ans);
}