51nod1040
2015-09-22 20:39
127 查看
链接:点击打开链接
题意:给出一个n,求1-n这n个数,同n的最大公约数的和。比如:n = 6。1,2,3,4,5,6 同6的最大公约数分别为1,2,3,2,1,6,加在一起等于15
代码:
题意:给出一个n,求1-n这n个数,同n的最大公约数的和。比如:n = 6。1,2,3,4,5,6 同6的最大公约数分别为1,2,3,2,1,6,加在一起等于15
代码:
#include <iostream> #include <cstdio> #include <stdlib.h> #include <cstring> #include <algorithm> using namespace std; long long phi(long long x) { long long i,ans=x; for(i=2;i*i<=x;i++) { if(x%i==0) ans=ans-ans/i; while(x%i==0) x/=i; } if(x>1) ans=ans-ans/x; return ans; } int main(){ //因为1~N中的数与N的最大公因数,一定都是N的因子 long long i,n,sum; //所以问题变为求∑xy(x代表N的因子的大小,y代表x的个数) while(scanf("%lld",&n)!=EOF){ sum=0; for(i=1;i<=n;i++){ if(i*i>n) break; if(n%i==0){ if(i*i==n) sum+=(phi(i)*i); //因为1~N中与N最大公因子为f的数的个数,与1~N/f中与N/f互质 else{ //的个数相等,所以f的个数为N/f的欧拉函数值 sum+=(phi(n/i)*i); sum+=(phi(i)*n/i); } } } printf("%lld\n",sum); } return 0; }
相关文章推荐
- 第4周—项目3.(1)结点逆置
- mycd命令
- 替代sql in 性能优化
- 关于FFTW3和NE10在树莓派2上执行FFT和IFFT的耗时比较
- C++函数的默认值问题
- c#特性加反射实现解耦
- 结对项目之需求分析与原型模型设计
- systemUI之statusbar
- Android 位置服务——用户定位User Location
- More Effective C++ 条款16 谨记80-20法则
- 关于网络字节序和主机字节序
- c/c++ 获取当前程序(EXE)所在的路径
- opencv在C++窗体程序中 导包出错
- Web.xml配置详解
- bzoj-3307 雨天的尾巴
- jqPlot,一个 jQuery这个 JavaScript 框架的绘图插件
- 如何在代码中动态为TextView/Button设置drawableRight
- MFC获取指针.doc
- hihoCoder 1231 Border Length(ACM-ICPC国际大学生程序设计竞赛北京赛区(2015)网络赛 E)
- 关于FFTW3和NE10源代码库的编译