bzoj 2705 [SDOI2012]Longge的问题

分析来自：

http://blog.csdn.net/cjk_cjk/article/details/46324349

http://blog.csdn.net/xaphoenix/article/details/50677870

题目大意：给定一个整数N，你需要求出∑gcd(i, N)(1<=i <=N)。0<N<=2^32

按gcd的不同，题目可以转化为：sigma(g*与n的最大公约数为g的数的个数)
先考虑与n的gcd=1的i，它的数量为 φ(n)
假设 gcd(i, N)=g,i=g*a,n=g*b，那么i与n同除以g后互质，也就是说：与n的最大公约数为g的数的个数为 φ(n/g)

可以在sqrt(n)范围内枚举n的因子，并求出其欧拉函数值

举个例子 比如100 那么对于10来说从1到10只要与10互质的数字都是可以乘以10与100成为一对，那么就好办了 累加所有的情况（分N/i和N/(N/i)）

AC代码：

#include<cmath>
#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<cstdlib>
#include<set>
#define LL long long
using namespace std;
LL N,M,T;

LL getphi(LL x){
	LL phi=1;
	LL i;
	for(i=2;i*i<=x;i++){
		if(x%i==0){
			x/=i;phi*=i-1;
			while(x%i==0){
				x/=i;phi*=i;
			}
		}
	}
	if(x>1)	phi*=x-1;
	return phi;
}	

int main()
{
	scanf("%lld",&N);
	LL s=sqrt(N);
	LL ans=0;
	for(LL i=1;i<=s;i++)
	{
		if(N%i==0)
		{
			if(i*i<N)ans+=i*getphi(N/i);
			ans+=getphi(i)*N/i;
		}
	} 
	
	printf("%lld\n",ans);
	return 0;
}