题目大意:给定一个整数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的因子,并求出其欧拉函数值
【代码】
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
typedef long long LL;
LL phi(LL m)
{
LL ans=m,sq=(LL)sqrt(m),i;
for(i=2;i<=sq;i++)
{
if(m%i==0) ans=ans/i*(i-1);
while(m%i==0) m/=i;
}
if(m>1) ans=ans/m*(m-1);
return ans;
}
int main()
{
LL n,i,j,m,ans=0;
scanf("%lld",&n);
for(i=(LL)sqrt(n);i>=1;i--)
{
if(n%i==0)
{
ans+=i*phi(n/i);
if(i*i<n) ans+=(n/i)*phi(i);
}
}
printf("%lld",ans);
return 0;
}
这里有与程序竞赛有关的数学知识点:https://blog.csdn.net/ctsas/article/details/55000318