2705: [SDOI2012]Longge的问题
Time Limit: 3 Sec Memory Limit: 128 MBSubmit: 2045 Solved: 1260
[ Submit][ Status][ Discuss]
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。
思路:
欧拉函数个一个重要性质,gcd(m,n)=k,所以gcd(m/k,n/k)=1,于是s(k)=euler(n/k),如果m,n的最大公约数是k,那么在n中,与n公约数为k的数字个数是euler(n/i),
因为如果一个数x跟n/i互质的话,说明x再乘以其他k时,跟n的公约数最大也就还是k,只要我们控制范围,就可以确定这些x的范围,在计算的时候我们应该考虑一个双向
性,因为如果 i 可以整除n,那么n/i也是可以整除n的,所以计算euler的时候顺带计算下民euler(i)*n/i,
AC-code:
#include<iostream>
#include<cstdio>
#include<math.h>
long long euler(long long x)
{
long long t=x;
for(int i=2;i<=sqrt(x);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()
{
long long n,m,ans;
int i;
scanf("%lld",&n);
m=sqrt(n);
ans=0;
for(i=1;i<=m;i++)
{
if(n%i==0)
{
ans+=i*euler(n/i);
if(i*i<n)
ans+=n/i*euler(i);
}
}
printf("%lld\n",ans);
return 0;
}