比赛 |
20120718 |
评测结果 |
AAAAAAAAAA |
题目名称 |
最大公约数和 |
最终得分 |
100 |
用户昵称 |
Czb。 |
运行时间 |
0.007 s |
代码语言 |
C++ |
内存使用 |
1.40 MiB |
提交时间 |
2012-07-18 10:48:44 |
显示代码纯文本
#include<stdio.h>
typedef long long Int;
Int n,t,ans,prime[20000];
bool flag[1000000];
Int ou(Int k)
{
Int i,s=k,tmp=k;
for(i=1;i<=t&&prime[i]*prime[i]<=k;i++)
{
if(k%prime[i]==0)s=s/prime[i]*(prime[i]-1);
while(k%prime[i]==0)k/=prime[i];
}
if(k>1)s=s/k*(k-1);
return s;
}
int main()
{
freopen("gcdsum.in","r",stdin);
freopen("gcdsum.out","w",stdout);
Int i,j;
scanf("%lld",&n);
for(i=2;i*i<=n;i++)
{
if(!flag[i])
{
prime[++t]=i;
for(j=i+i;j*j<=n;j+=i)
{
flag[j]=true;
}
}
}
for(i=1;i*i<=n;i++)
{
if(n%i==0)
{
if(i*i==n)
{
ans+=ou(i)*i;
break;
}
ans+=ou(i)*n/i;
ans+=ou(n/i)*i;
}
}
printf("%lld\n",ans);
return 0;
}