博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
51nod 1040最大公约数和(欧拉函数)
阅读量:4307 次
发布时间:2019-06-06

本文共 1323 字,大约阅读时间需要 4 分钟。

题目来源: 
基准时间限制:1 秒 空间限制:131072 KB 分值: 80 
 收藏
 关注
给出一个n,求1-n这n个数,同n的最大公约数的和。比如:n = 6
1,2,3,4,5,6 同6的最大公约数分别为1,2,3,2,1,6,加在一起 = 15
 
Input
1个数N(N <= 10^9)
Output
公约数之和
Input示例
6
Output示例
15
思路:
目的是求∑(i= 1,n) gcd( i , n );
gcd( i , n ) = x,表示x是n的因子。稍作变形gcd( i / x , n / x) = 1,
看到这个式子可以想到欧拉函数,也就是求比n/x小的与其互质的个数。
因为这些书和n/x互质,乘上x后与n的最大公约数只有x。
也就是说我们先求出每个因子,然后计算每个因子有多少贡献即可。
 
 
/* * Author:  sweat123 * Created Time:  2016/6/27 14:01:46 * File Name: main.cpp */#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define INF 1<<30#define MOD 1000000007#define ll long long#define lson l,m,rt<<1#define rson m+1,r,rt<<1|1#define pi acos(-1.0)using namespace std;const int MAXN = 1000000;int n;ll ef(int n) { ll cnt = n; int i; for(i = 2; i * i <= n; i++){ if(n % i == 0) { cnt -= cnt / i; // (x-x/p1) *(1-1/p2)*(1-1/p3)*(1-1/p4)..... while(n % i == 0) n /= i; } } if(n > 1) cnt -= cnt / n; return cnt; } int main(){ while(~scanf("%d",&n)){ ll ans = 0; for(int i = 1; i <= (int)sqrt(n); i++){ if(n % i == 0){ ans += ef(n / i) * i; if(n / i != i){ ans += ef(i) * (n / i); } } } printf("%lld\n",ans); } return 0;}

 

转载于:https://www.cnblogs.com/sweat123/p/5620086.html

你可能感兴趣的文章
Android中热修复框架AndFix原理解析及案例使用
查看>>
python3安装scrapy
查看>>
python正则表达式入门一
查看>>
python正则表达式入门二
查看>>
scrapy运行
查看>>
XPATH入门
查看>>
python爬虫 CSS选择器
查看>>
正常关闭java程序
查看>>
查看linux核心数
查看>>
数据结构与算法三: 数组
查看>>
Activiti工作流会签二 启动流程
查看>>
Activiti工作流会签三 撤销,审批,驳回
查看>>
Oauth2方式实现单点登录
查看>>
CountDownLatch源码解析加流程图详解--AQS类注释翻译
查看>>
ES相关度评分
查看>>
我们一起做一个可以商用的springboot脚手架
查看>>
idea在搭建ssm框架时mybatis整合问题 无法找到mapper
查看>>
java设计基本原则----单一职责原则
查看>>
HashMap的实现
查看>>
互斥锁 synchronized分析
查看>>