快速幂+逆元求组合数

原创 于 2024-12-08 22:24:46 发布 · 498 阅读 · 11 ·
CC 4.0 BY-SA版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
文章标签：

#c++ #算法 #开发语言

在计算组合数 C(n, k) = \frac{n!}{k!(n-k)!} 时，直接暴力计算既慢又容易溢出。今天我们来揭开 快速幂 和 模逆元 的神秘面纱，带你一边学习理论，一边轻松解决实际问题！

什么是快速幂？

快速幂是一种高效计算 a^b \mod p 的方法。它利用指数的二进制表示，巧妙地减少了乘法次数，把原本复杂的幂运算大幅加速。

快速幂的原理

假设你要计算 a^b，比如 2^{13}，可以写成：

2^{13} = 2^{8} \cdot 2^{4} \cdot 2^{1}

这里的 8,4,1 是 13 的二进制表示 1101。快速幂的核心思想就是通过分解指数，把幂运算拆分成多个平方和相乘的过程。

具体分解步骤如下：

  1. 如果指数 b 是奇数，取出当前底数 a 的值乘到结果中。
  2. 把底数 a 平方（相当于处理下一位二进制），指数 b 除以 2。</
最低0.47元/天 解锁文章
ZCMU1934: ly的二叉树（卡特兰数算法+快速幂+逆元的求法）
凡心所向，素履所往，生如逆旅，一苇以航
07-23 907
ly的二叉树： Description 某一天，ly正在上数据结构课。老师在讲台上面讲着二叉树，ly在下面发着呆。 突然ly想到一个问题：对于一棵n个无编号节点的有根二叉树，有多少种形态呐？你能告诉她吗？ Input 多组输入，处理到文件结束 每一组输入一行，一个正整数n(1≤n≤1000000),意义如题目所述。 Output 每组数据输出一行，包含一个正整数表示答案，由于数字可能非...
[AcWing]876. 快速幂求逆元（C++实现）快速幂求逆元模板题
cloud的博客
12-06 1930
[AcWing]875. 快速幂（C++实现）快速幂模板题1. 题目2. 读题（需要重点注意的东西）3. 解法4. 可能有帮助的前置习题5. 所用到的数据结构与算法思想6. 总结 1. 题目 2. 读题（需要重点注意的东西） 思路： 快速幂的作用 快速地求出a的k次方模上p的结果 快速幂的主要思想 预处理出logk个数，然后用这logk个数中的若干个数在O(1)的时间组合出a的k次方模上p的结果。 3. 解法 -------------------------------------------
快速幂 逆元
benTuTuT的博客
03-11 693
快速幂，逆元模板
快速幂求逆元
城南的花
04-18 1936
费马小定理 在模为**素数**p的情况下，有费马小定理 a^(p-1)=1（mod p） 那么a^(p-2)=a^-1(mod p) 也就是说a的逆元为a^(p-2) 而在模**不为素数**p的情况下，有欧拉定理 a^phi(m)=1（mod m） (a⊥m) 同理a^-1=a^(phi(m)-1) 因此逆元x便可以套用快速幂求得了x=a^(phi(m)-1) 但是似乎还有个问题...
快速幂_逆元_求组合数
qq_36706625的博客
04-22 487
快速幂：递归形式：static long pow_mod(long a, long n) { if (n == 0) { return 1; } long x = pow_mod(a, n / 2); long ans = x * x % MOD; if ((n &amp; 1) == 1) { ans = ans * a % MOD; } return ...
快速幂求组合数
zhe_zhenbuhui_a的博客
09-26 260
ll quick(ll a,ll b){ //计算a^b并且对mod取余 ll ans = 1; while(b){ if(b & 1) ans = ans * a % mod; a = a * a % mod; b >>= 1; } return ans % mod; } ll calc(ll n,ll m) {//求组合数 if(m > n) return 0; return (fac[n] .
用快速幂的方法求组合数
名字不能相同
11-14 706
题目链接???? 需要用到的知识 组合数 其中：C(r,n) = r! / (n!*(m-n)!) //默认n>r 判断一个数是否是奇数：if( n&1 ) 其中：n为奇数 // 奇数2进制末位为1，1与1进行与运算结果为1； 拓展： n>>=1 右移1位 除以2； n<<=1 左移1位 乘以2； 费马小定理 费马小定理：假如p是质数，且gcd(a,p)=1，那么 a^(p-1)≡1（mod p） 所以a*a^(p-2)≡1（mod p）
求组合数 ：DP、逆元+快速幂
qq_40586164的博客
08-30 452
要求：求 Cmn%pC_{m}^n \% pCmn​%p ，p为素数（经典p=1e9+7） 方法一：DP C[n][m]=CnmC_{n}^mCnm​ 公式Cnm=Cn−1m−1+Cn−1mC_{n}^m=C_{n-1}^ {m-1}+C_{n-1}^{m}Cnm​=Cn−1m−1​+Cn−1m​ 证明：班级中有n个人，选出m个人开除，有两种选法： 班长不开除，从剩下的人中开除m个人 开除班长，从剩下的人中开除m-1个人 const int MOD = 1e9 + 7; int C[1005][1005
快速入门之——乘法逆元求组合数
serendipity的博客
04-25 3632
逆元+快速幂+阶乘求组合数（快速入门） 文章目录逆元+快速幂+阶乘求组合数（快速入门）前言：何为逆元？怎么求逆元？怎么求阶乘逆元？相关代码实现： 前言： 大家基本上应该都知道用杨辉三角法求组合数C(n,m)（n为下标） 也就是直接暴力打表求法O(n*m)： for(int i=0;i<=n;i++){ c[i][0]=c[i][i]=1; for(int j=1;j<i;j++) ...
组合数计算（快速幂）
喵不理
11-06 1800
用快速幂更高效。。。 #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; const ll mod=1000000009; const ll maxn=100; ll c[12345
快速幂与组合数
m0_53726614的博客
08-10 375
#include<iostream> using namespace std; const int N=2010,mod=1e9+7; int c[N][N]; void init() { for(int i=0; i<N; i++) { for(int j=0; j<=i; j++) { if(!j) c[i][j]=1; else
快速幂 + 逆元
zhuibushixi的博客
08-10 358
逆元 乘法逆元： 有 a * x = 1 (mod c) 称 x 为 a 与 c 的乘法逆元。 理论依据：（mod 对于加减乘除都可进行分配） f / a mod c= ? (f mod c / a mod c = ?) 如果存在 a*x = 1(mod c) 那么两边同时乘起来，得到 f*x = ?(mod c) 成立条件： 1 模方程 a*x= 1(mod c)存在解 ...
求组合数（快速幂+逆元）
Hello Will
04-18 354
/** * 快速幂取模算法 * 分析：a^b mod c = (a mod c)^b mod c * a^b mod c = (a^2)^b/2 mod c , b为偶数 * a^b mod c = ((a^2)^b/2·a) mod c , b为奇数 * 每次让指数折半 */ private static final l...
快速幂 逆元 求组合数
Minoz的博客
08-12 627
我们知道 如果 a*x = 1 那么x是a的倒数，x = 1/a 但是a如果不是1，那么x就是小数 那数论中，大部分情况都有求余，所以现在问题变了 a*x  = 1 (mod p) 那么x一定等于1/a吗 不一定 所以这时候，我们就把x看成a的倒数，只不过加了一个求余条件，所以x叫做    a关于p的逆元，a和p互质，a才有关于p的逆元。 a的逆元，我们用inv(a)来表示 ...
快速求组合数
weixin_40136732的博客
05-13 5934
摘自https://www.jianshu.com/p/718a5ac26238 逆元+快速幂解法 （一）基本概念 上面两种方法都使用了递归方法，递归方法有个缺陷，就是在数据较大时效率较低。所以这里要介绍一个种新的求组合算法。在了解此算法之前，要先了解一些概念。 1 同余 同余是数论中的重要概念。 给定一个正整数m，如果两个整数a和b满足a-b能够被m整除，即(a-b)/m得到一个整数，那么就称整...
快速幂求逆元c++模版
最新发布
08-24
### 快速幂算法在C++中求解逆元的方法 在模运算中，当模数 $ p $ 为质数时，可以使用**费马小定理**来求解乘法逆元。根据费马小定理，若 $ a $ 与 $ p $ 互质，则 $ a^{p-2} \mod p $ 即为 $ a $ 关于模 $ p $ 的逆元[^2]。为了高效计算 $ a^{p-2} \mod p $，通常采用**快速幂算法**。 快速幂算法的核心思想是利用**分治法**和**二进制分解**来减少幂运算的计算次数。其基本原理是将指数 $ b $ 分解为二进制形式，并在每一步中将底数平方，从而在 $ O(\log b) $ 的时间复杂度内完成计算。 以下是一个通用的 C++ 模板代码，用于通过快速幂算法求解逆元： ```cpp #include <iostream> using namespace std; typedef long long LL; LL qmi(LL a, LL b, LL p) { LL res = 1; while (b) { if (b & 1) res = res * a % p; a = a * a % p; b >>= 1; } return res; } int main() { int n; cin >> n; while (n--) { int a, p; cin >> a >> p; if (a % p == 0) { cout << "impossible" << endl; } else { cout << qmi(a, p - 2, p) << endl; } } return 0; } ``` 该代码中，`qmi` 函数实现了快速幂算法，用于计算 $ a^b \mod p $。在求逆元时，调用 `qmi(a, p - 2, p)` 来计算 $ a^{p-2} \mod p $。若 $ a $ 与 $ p $ 不互质（即 $ a \mod p = 0 $），则说明逆元不存在，输出 `impossible` [^4]。 ### 应用场景 - **密码学**：在 RSA 加密算法中，常需要计算模逆元来生成密钥。 - **组合数学**：在计算组合数的模运算时，常使用快速幂求逆元以避免除法。 - **竞赛编程**：快速幂是处理大指数幂模运算的标准工具，广泛应用于各类算法竞赛中。 ### 注意事项 - 快速幂算法适用于模数 $ p $ 为质数的情况，若 $ p $ 不是质数，则需使用扩展欧几里得算法或线性筛法求解逆元。 - 在使用快速幂算法时，必须确保 $ a $ 与 $ p $ 互质，否则逆元不存在。 - 若 $ a $ 与 $ p $ 不互质，则输出 `impossible`，表示不存在逆元。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值