POJ2480

POJ2480

#include<cstdio>
#define ll long long
ll phi(ll x)
{
	ll ans=x;
	for(ll i=2;i*i<=x;i++)
	{
		if(x%i==0)
		{

			ans=ans/i*(i-1);
			while(x%i==0)
				x/=i;
		}
	}
	if(x>1)
		ans=ans/x*(x-1);
	return ans;
}
int main()
{
	ll n;
	while(scanf("%lld",&n)!=EOF)
	{
		ll ans=0;
		for(ll i=1;i*i<=n;i++)
		{
			if(n%i==0)
			{
				ans+=(i*phi(n/i));
				if(i*i!=n)
				{
					ans+=(n/i*phi(i));
				}
			}
		}
		printf("%lld\n",ans);
	}
	return 0;
}