分析来自:
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;
}