- 综合排序
- 最热优先
- 最新优先
- 来自专栏花猪的学习记录
数论基础（部分计算方法）
正文 模重复平方计算 例1：计算 68879 mod 3337 例2：计算 97263533 mod 11413 所以 97263533 mod 11413 = 5761 扩展欧几里得计算 例：计算
1.3K10编辑于 2022-02-22 - 来自专栏Coggle数据科学
数论及数论四大定理
请勿转载 数论及数论四大定理hankin2015.github.io ? 数论 （数学分支） 数论是纯粹数学的分支之一，主要研究整数的性质。整数可以是方程式的解（丢番图方程）。 按研究方法来看，数论大致可分为初等数论和高等数论。初等数论是用初等方法研究的数论，它的研究方法本质上说，就是利用整数环的整除性质，主要包括整除理论、同余理论、连分数理论。 高等数论则包括了更为深刻的数学研究工具。它大致包括代数数论、解析数论、计算数论等等。 2、欧拉定理 在数论中，欧拉定理（Euler Theorem，也称费马-欧拉定理或欧拉函数定理）是一个关于同余的性质。欧拉定理得名于瑞士数学家莱昂哈德·欧拉，该定理被认为是数学世界中最美妙的定理之一。 i++) { if(gcd(i, n) == 1) cnt++; } cout << cnt << endl; //根据通式计算
3.9K10发布于 2019-09-12 - 来自专栏blog(为什么会重名，真的醉了)
数论-素数
文章目录 判断素数 筛法求素数 例题 HDU-1262 HDU-3792 判断素数 ---- 枚举 [2 ， x ] bool prime(int x) { if (x <= 1)return false int sieve(int x) { for (int i = 0; i <= x; i++)vis[i] = false;//初始化 for (int i = 2; i * i <= x; i+ if (a % 2 == 0)a--; int b = m - a; while (prime(a)==false||prime(b)==false) { a -= 2; b Twin Prime Conjecture states that “There are infinite consecutive primes differing by 2”. Sample Input 1 5 20 -2 Sample Output 0 1 4 若两个素数相差2则称为一对孪生素数，求区间[1,n]内的孪生素数个数。
99930发布于 2020-09-15 - 来自专栏陈黎栋的专栏啦
数论——快速幂算法 快速计算a^b mod c的值
// 快速计算 (a ^ p) % m 的值 __int64 FastM(__int64 a, __int64 p, __int64 m){ if (p == 0) return 1;
1K40发布于 2020-02-18 - 来自专栏全栈程序员必看
数论 欧拉函数_数论欧拉函数
欧拉函数的通式：φ(n)=n*(1-1/p1)(1-1/p2)(1-1/p3)*(1-1/p4)……(1-1/pn) 其中p1, p2……pn为n的所有质因数，n是不为0的整数。 所以，根据通式我们可以打出以下代码： ll eular(ll n) { ll ans = n; for(int i=2; i*i <= n; ++i) { void euler() { for(int i=2;i<maxn;i++){ if(! （特别地，当p为质数时，phi（p）=p-1,此时逆元为x^(p-2)，即费马小定理） ④ 当n为奇数时，phi(2n)=phi(n) ⑤ 若x与p互质，则p-x也与p互质，因此小于p且与p互质的数之和为 phi(x)*x/2; ⑥N>1，不大于N且和N互素的所有正整数的和是 1/2 *N *eular(N)。
68520编辑于 2022-09-23 - 来自专栏奇妙的算法世界
HDOJ 1018（数论）
Sample Input 2 10 20 Sample Output 7 19 思路 刚看到题的一瞬间，果断写了个高精度交了上去，TLE了，再看数据范围，嚯！ 1e7，然后就去翻了题解，发现是数论问题，求阶乘位数有两种方法： 1.10m<n!<10(m+1) 若求得M，则M+1为答案。对方程两边以10为底求对数，得M<log10(n!) <M+1，通过循环求值即可 2.斯特林公式：n! ≈ sqrt(2* n * pi)* (n/e)^n，则 M+1=(int)(0.5 * log(2.0 * n *PI)+n * log(n)-n)/(log(10.0)) )+1; AC代码 方法 for(int i=1;i<=n;i++){ ans+=log(i)/log(10); } cout<<(int)ans+1<<endl; } return 0; } 方法2
52820发布于 2020-10-23 - 来自专栏CSDN旧文
数论--模板整理
数论–康托展开与逆康托展开模板 数论–组合数（卢卡斯+扩展卢卡斯）模板 数论–Miller_Rabin判断素数 数论–中国剩余定理模板 数论–逆元（拓展欧几里得）模板 数论–逆元（费马小定理）模板 数学–数论–因子和线性筛 （模板) 数学–数论–随机算法–Pollard Rho 大数分解算法（纯模板带输出） 数学–数论–快速幂–最大公约数–位运算模板 线性筛求积性函数的模板 数学–图论–莫比乌斯线性筛模板 数学–数论—欧拉筛 模板 数学–数论–素数
41210发布于 2020-10-28 - 来自专栏码神随笔
素数判断——数论
GetUglyNumbur(1500); return 0; } 但是很遗憾没有拿满分，翻开我那几乎积灰的剑指offer，看到这个题是放到了用空间换时间的算法中，又想了想，之所以会超时，是因为上面的题解中计算了许多不是丑数的数据 ，再看丑数的定义是，应该是另一个丑数乘以2，3，5的结果，进行优化 int GetUglyNumber2(int n) { if (n <= 0) return 0; int *pUglyNumber * 2, *pMulitiply3 * 3, *pMulitiply5 * 5); pUglyNumber[nextUglyNumber] = mmin; while (*pMulitiply2 * 2 <= pUglyNumber[nextUglyNumber]) ++pMulitiply2; while (*pMulitiply3 * 3 <= pUglyNumber[nextUglyNumber int number3) { int min = (number1 < number2) ?
47720编辑于 2022-12-13 - 来自专栏CSDN旧文
数学--数论--素数
定义判断： bool isPrime (int n) { for(int i=2;i*i<=n;i++) { if(n%i==0) return false; } else return false memset(bprime,false,sizeof(bprime)); bprime[0]=true; bprime[1]=true; for(int i=2; bprime[i]){ prime[cnt++]=i; for(LL j=i*2;j<=n;j+=i) bprime[j] cnt; bool bPrime[N]; void getPrimes(int n){ memset(bPrime,false,sizeof(bPrime)); for(int i=2; namespace std; bitset<100000010>v; int prime[6000001]; int m=0; void primes(int n) { for(int i=2;
60810发布于 2020-11-05 - 来自专栏bigsai
基础数论总结
这个到除了前面几次的计算量多一点，到后面随着数据增大对整个复杂度相加是比较小的，算法复杂度为O(nloglogn);别小瞧多的这个logn，数据量大一个log可能少不少个0，那时间也是十倍百倍甚至更多的差距 计算方法： 计算n的分解方式。主要是通过数的自身对从最小的质数开始整除除一直到不能整除，直到跳出限制条件。 你可以从2到n;逐个遍历判断，满足条件的话就在数组中添加对应的count。 当然，每被计算一次的时候，这个数就要被除一次。 上面方法对于大的数据显然复杂度太高。 这里不是按照次幂计算的，而是按照实打实的一个一个数判断的。 根据Xzhila的传统， 竹子的分数=Φ（竹子的长度） （Xzhilans非常喜欢数论）。对于您的信息，Φ（n）=小于n的数字，它们相对于素数（除了1之外没有公约数）到n。
1.1K30发布于 2019-09-24 【HPUoj】Bet（数论）
比如样例 N = 2 的时候，赔率比分别是1, 2，你有1000块钱，那么买 第一个 2000/3，后一个 1000/3，这样最坏情况下你的回报是 666.67 。 输入 多组数据。 样例输入 2 1 2 1000 样例输出 666.67 正好是2016年ICPC-ECfinal做过的一道题，比赛的时候卡了一个精度，用数组模拟就行了。 如果我们想让最低的收益最大，那么就让a1*m1 = a2*m2 = …… = ai * mi。 这样每一项的收益相同，最低的就最大，满足条件。 然后我们可以得到： a1/a2 = m2 / m1 a1/a3 = m3 / m1 a1/a4 = m4 / m1 …… a1/an = mn / m1 然后我们把左右相加： a1/a2 + a1/a3 + a1/a4 + …… + a1/an = (m2+m3+m4+…+mn)/ m1 继续变形 a1/a2 + a1/a3 + a1/a4 + …… + a1/an = (M-m1)/ m1 我们发现未知数只有一个
26510编辑于 2025-08-27- 来自专栏月亮与二进制
C++函数论
函数之前声明一个函数原型，表明函数的返回值类型、函数名、参数类型、参数名，如下： void swap (int a, int b); int main (void) { int a = 1; int b = 2; 因为在调用时你要写参数，不可能参数空一个不写，而写完了你要设置的，剩下的就都是右边默认的了： int func1 (int n, int m = 4, int j = 5); // 有效 int func2 (int n, int m = 6, int j); // 无效 // 不允许这样调用： func(1, , 2);// 无效 函数重载 c++允许有多个同名函数，只要其参数的类型或者数量不一样就可以了
67810编辑于 2022-01-07 - 来自专栏奇妙的算法世界
codeforces 573A （数论）
给你几个数字，你可以对这些数字进行无限次的乘2和乘3操作，询问你最后是否能将它们变为相同 思路 假设最后相同的数字为k，则得到式子：2x3xz=k，我们只需要将每个数所有的2和3去除，最后判断数字是否相同即可 int> PII; typedef pair<long,long> PLL; typedef pair<char,char> PCC; typedef long long LL; const int N=2* 10; const int M=150; const int INF=0x3f3f3f3f; const int MOD=998244353; int a[N]; int ff[4],ft[4]={1,2,3,6 int n;cin>>n; for(int i=1;i<=n;i++) cin>>a[i]; for(int i=1;i<=n;i++){ while(a[i]%2= =0) a[i]/=2; while(a[i]%3==0) a[i]/=3; } string ans="Yes"; for(int i=1;i<n;i++) if
42810发布于 2020-10-23 - 来自专栏C++
【算法】数论与数学
1、最大公约数（GCD）、最小公倍数（LCM） 2、素数判定与筛法 3、快速幂算法 4、数位操作 数字统计 BC153 [NOIP2010]数字统计 很惭愧，第一次做这个题没做出来。 ; i <= r; i++) { int tmp = i; while (tmp) { if (tmp % 10 == 2)
16900编辑于 2025-03-15 - 来自专栏叶子的开发者社区
三数论大小（引用）
输入 第一行输入t表示有t个测试实例 第二行起，每行输入三个整数 输入t行 输出 每行按照从大到小的顺序输出每个实例，三个整数之间用单个空格隔开 输入样例1 3 2 4 6 88 99 77 111 333 222 输出样例1 6 4 2 99 88 77 333 222 111 思路分析 对于这么一道简单的题目
22610编辑于 2023-07-28 - 来自专栏编程驿站
C++初等数论
除了理解数论概念，更重要能融会贯通。把对数论相关知识的认知运用到编程领域。 2. 同余式 概念 如果两个整数a,b 的差值除另一个整数(m)的值为一个整数，同称a,b对模m同余数。 余数判别法 基本思想：求N被m除的余数，先找到一个较简单的数R，使得N与R对于除数m同余．由于R是一个较简单的数，所以可以通过计算R被m除的余数来求得N被m除的余数。 同余关系是数论中的一种等价关系。数学上使用符号≡表示。同余类指模 m同余的所有整数的集合称为同余类。如所有偶数模2的余数都为0，可称所有偶数为同余类。剩余类是同余类的另一种叫法。 7.模运算意义下的逆元 在信息学竞赛中，当答案过于庞大的时候，我们经常会使用到模运算(Modulo Operation)来缩小答案的范围，以便输出计算得出的答案。 是数论中一个重要定理。又称中国余数定理。
94300编辑于 2024-03-11 - 来自专栏mythsman的个人博客
数论基础专题小结
int>prime; bool p[10000004]; int main(){ int tt; prime.clear(); memset(p,0,sizeof p); for(int i=2; i<=10000000;i++){ if(p[i]==0){ prime.push_back(i); for(int j=2;i*j<=10000000;j++){ p[j*i] ("%d",&tt); for(int t=1;t<=tt;t++){ int n; scanf("%d",&n); int ans=0; for(int i=0;prime[i]*2<
30710编辑于 2022-11-14 - 来自专栏owent
数论模板（个人模板）
} return tp; } 1.欧拉函数 Ψ（n） = 小于n且与n互质的数的个数 int eular(int n) { int res = 1, i; for(i = 2; [i]; } return i; } // 排列 // 参数: A(m,n),m >= n , p[] 传出数的数组表示指针 // 返回值：结果包含的素数个数 int A_Cache[2] return r; } // 取模计算:参数: mod为取模的值，其他参数同上 // 全排列取模 int Arrangement_Mod(int n, int p[], int mod) { [i]; mat[r2][i] = mat[r1][i] - mat[r2][i]; mat[r1][i] = mat[r1][i] - mat[r2][ i]; } } //高斯消元（整数） //返回0为有无穷解或无解，返回1有唯一解并计算答案，返回-1无解 bool gauss(int n) {
3.3K40发布于 2018-08-01 - 来自专栏全栈程序员必看
数论——欧拉函数
定义 小于n的正整数中与n互质的数的数目（φ(1)=1） 通式 证明： 设p是N的质因子，1~N中p的倍数有p,2p,3p,…,(N/p)*p，共N/p个。 2.根据算数基本定理，N=p1C1*p2C2*…*pkCk 。设N的最小质因子为p，当p的指数为1时，phi(N)=(p-1)*phi(N/p)。 3. 当p的指数不为1时，同2可证得phi(N)=p*phi(N/p)。 2的证明： 根据欧拉函数通式， phi(N)=N*(p1-1)/p1*(p2-1)/p2*…*(pk-1)/pk， phi(N/p1)=N/p1*(p2-1)/p2*…*(pk-1)/ 0)x/=i; } } if(x>1)res=res/x*(x-1); return res; } 线性筛法 根据前面的欧拉线性筛质数的算法（可参考本人博客：数论
63510编辑于 2022-09-06 - 来自专栏CSDN旧文
数学--数论--鸽巢原理
输出的第一行是选择元素的个数M，接着M行分别是选择的元素的值 由鸽笼原理可知此题一定有解，不存在输出0的结果 分析： 我们可以依次求出a[0],a[0]+a[1],a[0]+a[1]+a[2],…,a[ 0]+a[1]+a[2]…+a[n]； 假设分别是sum[0],sum[1],sum[2],…,sum[n] 如果在某一项存在是N的倍数，则很好解，即可直接从第一项开始直接输出答案 但如果不存在，则sum [i]%N的值必定在[1,N-1]之间，又由于有n项sum，有抽屉原理： 把多于n个的物体放到n个抽屉里，则至少有一个抽屉里有2个或2个以上的物体。
1.1K10发布于 2020-11-06
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号