HYSBZ 2705: [SDOI2012]Longge的问题【欧拉函数】

2705: [SDOI2012]Longge的问题

Time Limit: 3 Sec   Memory Limit: 128 MB
Submit: 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;
}