题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2705
题目大意:给定一个N,要求sigma(gcd(i,N))的值,1<=i<=N。
题目分析:我们枚举gcd(i,N)的值,即,枚举N的约数,时间复杂度为sqrt(N),然后对于每一个枚举值p,我们计算phi(N/p),由于N的值比较大,我们不能预处理,所以我们就暴力在sqrt(N)的时间内求出phi(N/p),那么gcd(i,N)为p的那一部分的和就等于p*phi(N/p)。我们把这些结果累加起来即可。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#define maxn 11000000
using namespace std;
typedef long long LL;
LL d;
int n;
int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
LL phi(LL x)
{
LL t=x;
for(LL i=2;i<=d;i++)
if(x%i==0)
{
t=t/i*(i-1);
while(x%i==0)x/=i;
}
if(x>1)t=t/x*(x-1);
return t;
}
int main()
{
scanf("%d",&n);
LL ans=0;
d=sqrt(n);
for (int i=1;i<=d;i++)
{
if (n%i==0)
{
ans+=phi(i)*n/i;
if (i!=n/i) ans+=phi(n/i)*i;
}
}
cout<<ans<<endl;
return 0;
}