码迷,mamicode.com
首页 > 其他好文 > 详细

BZOJ 2705 [SDOI2012]Longge的问题(欧拉函数)

时间:2017-05-18 01:13:15      阅读:182      评论:0      收藏:0      [点我收藏+]

标签:problem   bsp   代码   ble   for   names   函数   target   php   

 

【题目链接】 http://www.lydsy.com/JudgeOnline/problem.php?id=2705

 

【题目大意】

  求出∑gcd(i,N)(1<=i<=N)

 

【题解】

  $∑_{i=1}^{N}gcd(i,N)$
  $=∑_{i=1}^{N}∑_{d|gcd(i,N)}\phi(d)$
  $=∑ \phi(d)∑ _{1=<i<=N \land d|i \land d|N}1$
  $=∑_{d|N}\phi(d)\lfloor\frac{i}{d}\rfloor$

 

【代码】

#include <cstdio>
#include <algorithm>
using namespace std;
int Euler(int n){
    int t=1,i;
    if(!(n&1))for(n>>=1;!(n&1);n>>=1,t<<=1);
    for(i=3;i*i<=n;i+=2)if(n%i==0)for(n/=i,t*=i-1;n%i==0;n/=i,t*=i);
    if(n>1)t*=n-1;
    return t;
}
int main(){
    int n;
    long long ans=0;
    while(~scanf("%d",&n)){
        for(int i=1;i*i<=n;i++){
            if(n%i==0){
                ans+=1LL*i*Euler(n/i);
                if(i*i<n)ans+=1LL*(n/i)*Euler(i);
            }
        }printf("%lld\n",ans);
    }return 0;
}

BZOJ 2705 [SDOI2012]Longge的问题(欧拉函数)

标签:problem   bsp   代码   ble   for   names   函数   target   php   

原文地址:http://www.cnblogs.com/forever97/p/bzoj2705.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!