nyoj 998(欧拉定理的运用)
Sum
- 描述
-
给你一个数N,使得在1~N之间能够找到x使得x满足gcd( x , N ) >= M,
求解gcd(x,N)的和
- 输入
-
多组测试数据
每行输出两个数N,M(N,M不超int)
- 输出
- 输出sum
- 样例输入
-
5 3
- 样例输出
-
5
-
-
解题思路:假设gcd(x,n) =k >= m,那么k*gcd(x/d,n/d) = k。也就是说,x/d与n/d是互质的,它们的gcd是1,再乘以k那当然就是gcd(x,n)啦。。那么首先就是枚举n的因子,在利用欧拉定理求出小于等于n/d且与之互质的数的个数有多少。有点点绕,但是想明白还是比较容易的
-
AC:
-
#include<iostream> #include<cstdio> using namespace std;typedef long long LL; LL Euler(LL n) {LL ans = n;for(int i = 2; i * i <= n; i++){if(n % i == 0){ans = ans / i * (i-1);while(n % i == 0)n /= i;}}if(n > 1) ans = ans / n * (n-1);return ans; }int main() { LL n,m;while(cin>>n>>m){LL ans = 0;for(int i = 1; i * i <= n; i++){if(n % i == 0){if(i >= m){int d = i;ans += d*Euler(n/d);}if(i * i != n && n / i >= m){int d = n / i;ans += d*Euler(n/d);}}}cout<<ans<<endl;}return 0; }
nyoj 998(欧拉定理的运用)相关推荐
- NYOJ 998 解题报告
Sum 时间限制:1000 ms | 内存限制:65535 KB 难度:3 描述 给你一个数N,使得在1~N之间能够找到x使得x满足gcd( x , N ) >= M, 求解gcd(x, ...
- 欧拉定理 费马小定理
前言 学基础数论的时候看过证明,然而很快就忘了,最近在学习高深一点的数论,于是再复习一下欧拉定理和费马小定理. 欧拉定理 内容 若正整数 \(a,n\) 互质,则 \(a^{\varphi(n)}\e ...
- 「欧拉定理」学习笔记(费马小定理)
欧拉定理:对于互质的两个正整数$a, n$,满足$a^{φ(n)} ≡ 1\ (mod\ n)$ 证明: 设集合$S$包含所有$n$以内与$n$互质的数,共有$φ(n)$个:$$S = \{ x_1 ...
- 费马定理中值定理_数论-欧拉函数、欧拉定理
欧拉函数 积性函数 满足 ( 互质) 定义 对于正整数 ,欧拉函数是小于等于 的所有数中与 互质的数的 个数. 欧拉函数是积性函数(这个证明不是很显然,这个链接里面有很多种证明方法) 记作: 公式 , ...
- UVA1342 That Nice Euler Circuit(ACM - ICPC 2004 Asia - Shanghai)(计算几何、欧拉定理)
整理的算法模板合集: ACM模板 欧拉定理:设平面图的顶点数.边数和面数分别为V,E,F,则V+F-E=2. #include<bits/stdc++.h> using namespace ...
- 关于素数常用结论--威尔逊定理、欧拉定理、费马小定理、米勒罗宾算法
再需要判定的数比较大时,用枚举法肯定不行的,但目前数学界也没有任何一种又快又准确的判定素数的方法,并且也证明了素数不存在任何一种通项表达式.但作为初等数论中最大的一部分内容,数学家们对素数性质进行了大 ...
- NYOJ 30 Gone Fishing JAVA+解析
Gone Fishing这道题目运用的多次折合成一次这种思想我首次见,我想的一个思路是,每次算一下鱼量和时间代价比,这个代码我没有敲,下面的代码是一位仁兄敲得,我研读了一下,做了一个注释,应该有利于后 ...
- 中gcd函数_欧拉函数φ(n)的计算及欧拉定理
1 欧拉函数定义 在数论中,对正整数n,欧拉函数φ(n)是小于或等于n的正整数中与n互质的数的数目.此函数以其首名研究者欧拉命名,它又称为φ函数(由高斯所命名)或是欧拉总计函数(totient fun ...
- Manacher算法 , 实例 详解 . NYOJ 最长回文
51 Nod http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1089 Manacher 算法 定义数组 p[i]表示以i为 ...
最新文章
- 为什么说C语言和linux是分不开的?
- 如何使用React.js和Heroku快速实现从想法到URL的转变
- 按摩师-总预约时间最长
- Scrapy框架中的crawlspider爬虫
- Safari 不能播放Video ,Chrome等可以 问题解决。
- intellij idea 1314 插件推荐及快速上手建议
- (解题思路)Entity Framework 如动态创建表或者列
- 微信重大更新!这特么是为上班摸鱼开发的吧.....(附内测地址)
- Hi,我们是-MobileNet-家族
- 关于 Win 8 RT 你应该知道的 15 件事
- HDOJ 1698 Just a Hook(线段树成段更新)
- 老白聊数据-关于销售预测的那些事
- eplan2022启动报错存储空间不够或者提示没有settings权限
- 疯狂语音红包 5.2.0 微信红包引流吸粉
- 宿命java_昆仑奇侠传外传之宿命篇
- 推荐凸优化经典书籍,来自boyd
- 7月书讯 | 12本新书如期而至!
- 创建vue3项目(vite+js)
- C/C++找工作这点事
- 百度地图离线开发2.0