51nod 1040 最大公约数之和

#include<bits/stdc++.h>
using namespace std;

long long oula(long long n)
{
	long long rea=n;
	for(long long i=2; i*i<=n; i++)
	{
		if(n%i==0)
		{
			rea=rea-rea/i;
			while(n%i==0)
				n/=i;
		}
	}
	if(n>1)
		rea=rea-rea/n;
	return rea;
}

int main()
{
	long long n,ans,i,limit;
	while(~scanf("%lld",&n))
	{
		limit=sqrt(n);
		ans=0;
		for(i=1;i<=limit;i++)
		{
			if(n%i==0)
				ans+=i*oula(n/i)+(n/i)*oula(i);
		} 
		if(limit*limit==n)
			ans-=limit*oula(limit);
		printf("%lld\n",ans);
	}
}