Description

Longge的數(shù)學成績非常好，并且他非常樂于挑戰(zhàn)高難度的數(shù)學問題。現(xiàn)在問題來了：給定一個整數(shù)N，你需要求出∑gcd(i, N)(1<=i <=N)。

Input

一個整數(shù)，為N。

Output

一個整數(shù)，為所求的答案。

Sample Input

6

Sample Output

15

HINT

?

【數(shù)據(jù)范圍】

對于60%的數(shù)據(jù)，0<N<=2^16。

對于100%的數(shù)據(jù)，0<N<=2^32。

———————————————————————

這道題如果一個數(shù)x gcd(n,x)==y 那么gcd(b/y,x/y)==1

所以我們枚舉因數(shù)d 求一下1-n/d有多少個數(shù)和n/d的gcd為1 這個可以用歐拉函數(shù)

#include<cstdio> #include<cstring> #include<algorithm> #define LL long long const int M=1e3+7; LL read(){LL ans=0,f=1,c=getchar();while(c<'0'||c>'9'){if(c=='-') f=-1; c=getchar();}while(c>='0'&&c<='9'){ans=ans*10+(c-'0'); c=getchar();}return ans*f; } LL n,v,ans; int p[M],cnt; LL f(LL x){for(int i=1;i<=cnt;i++)if(x%p[i]==0) x=x/p[i]*(p[i]-1); return x;} int main(){n=read(); v=n;for(LL x=2;x*x<=v;x++)if(v%x==0){p[++cnt]=x;while(v%x==0) v/=x;}if(v!=1) p[++cnt]=v;for(LL x=1;x*x<=n;x++)if(n%x==0){LL y=n/x;ans=ans+y*f(x);if(x!=y) ans=ans+x*f(y);}printf("%lld\n",ans);return 0; } View Code

?

轉載于:https://www.cnblogs.com/lyzuikeai/p/7664991.html

總結

以上是生活随笔為你收集整理的bzoj 2705: [SDOI2012]Longge的问题——欧拉定理的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。