Description
Longge的数学成绩非常好,并且他非常乐于挑战高难度的数学问题。现在问题来了:给定一个整数N,你需要求出∑gcd(i, N)(1<=i <=N)。
Input
一个整数,为N。
Output
一个整数,为所求的答案。
Sample Input
6
Sample Output
15
HINT
【数据范围】对于60%的数据,0<N<=2^16。
对于100%的数据,0<N<=2^32。
枚举n的每一个约数,我们假设a(i)为1~n中与n的最大公约数为i的个数,则每找到一个约数,我们就把ans加上(i*a(i))(到这里还能自己搞出来但是到下面就不会了QAQ还是感谢HZWER神犇的帮助QUQ)
因为gcd(m,n)=k,所以gcd(m/k,n/k)=1,于是s(k)=phi(n/k)
最重要的一条啊。。。(欧拉函数什么的,就是不会用╭(╯^╰)╮)
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
int n;
long long m;
long long ans;
long long phi(long long x)
{
long long t=x;
for (long long i=2;i<=m;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);
m=sqrt(n);
for (int i=1;i<=m;i++)
{
if (n%i==0)
{
ans+=i*phi(n/i);
if (i*i<n) ans+=n/i*phi(i);
}
}
printf("%lld",ans);
}