还是逆元之O(n)阶乘逆元。。。

最新推荐文章于 2022-07-30 21:06:45 发布
原创 最新推荐文章于 2022-07-30 21:06:45 发布 · 5.4k 阅读 · 3 ·
CC 4.0 BY-SA版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文分享了一种高效计算阶乘逆元的方法，利用快速幂求最末项逆元，再递归计算整个序列，提高了组合数计算的速度。

除草

做一个题发现了一个逆元的知识盲点，就是阶乘的逆元

然后发现了可以这样

fac[0]=fac[1]=1;

for(int i=2;i<=MAXN;i++)fac[i]=fac[i-1]*i%mod;

inv[MAXN]=quipow(fac[MAXN],mod-2);

for(int i=MAXN-1;i>=0;i--)inv[i]=inv[i+1]*(i+1)%mod;

发现这个以后惊呆了

就是那种这明明没有什么特别难理解的道理可是就是忘了那些最基本的道理的感觉

继续努力

窝嘞割草
关注
阶乘的逆元
qq_52142857的博客
10-29 1880
求对阶乘 1! 到 n! 的逆元 一般写法 typedef long long ll; ll quick_mi(ll m,ll k) { ll res=1; while(k) { if(k&1)res=res*m%mod; m=m*m%mod; k>>=1; } return res; } infact[0]=1; for(int i=1;i<=n;i++) { infact[i]=infac
阶乘逆元 三种 快速 方法 费马小定理 递推 线性
mrcrack的博客
03-31 2868
费马小定理 #include<cstdio> typedef long long LL; const LL MOD = 1e9 + 7; LL fac[1000000+5]; //阶乘 LL inv[1000000+5]; //逆元 LL quickMod(LL a,LL b) { LL ans = 1; while (b) { ...
O(n)递推求乘法逆元和阶乘逆元
yhhy666的博客
07-21 388
i−1∗i≡1(modp)p=k∗i+r(k=⌊ip​⌋)(r=pmodi)k∗i+r≡0(modp)k∗r−1+i−1≡0(modp)i−1≡−⌊ip​⌋∗pmodi(modp)i−1≡(p−p/i)∗inv[pmodi](modp)阶乘逆元inv[i+1]=(i+1)!1​。
阶乘和阶乘逆元
weixin_30483013的博客
08-08 213
扫盲。今天做题才知道这玩意。。（那你之前是怎么算阶乘的哇。。只会暴力暴力暴力嘛。。。。） 阶乘逆元，原理是费马小定理，有些意思。 1 #include<cstdio> 2 typedef long long LL; 3 const LL MOD = 1e9 + 7; 4 LL fac[1000000+5]; //阶乘 5 LL inv[10000...
【阶乘逆元】【线性求逆元】【组合计数】牛妹的数学难题
行码棋
02-13 2047
线性求逆元， 阶乘逆元
乘法逆元(阶乘逆元+排列组合).cpp
08-14
c/c++相关代码（cpp），乘法逆元有关模板整合（附有注释），顺便整合了了阶乘逆元，以及排列组合计算和Lucas的模板代码。
c++ 乘法逆元
KonjakuLAF的博客
01-26 1415
主要参考：OI-WIKI 为什么要逆元 当一个题目让你求方案数时常要取余，虽然\((a+b)\% p=(a\% p+b\% p)\%p\)\((a-b)\% p=(a\% p-b\% p)\%p\)\((a\times b)\% p=(a\%p\times b\%p)\%p\) 但是\((\dfrac{a}{b})\%p\ne(\dfrac{a\%p}{b\%p})\%p\) 由于会出现这种情况...
逆元的线性递推求解方法及阶乘逆元
ComingLiu的博客
03-11 5059
我们已经学过了利用扩展欧几里得和费马小定理求解乘法逆元的方法，但是对于这道题，这两种方法都失效，那么我们需要找到一种线性的算法来求解，时间复杂度才合格，如果对于逆元概念存在问题可以看一下第一个链接 可以这样考虑， ...
记一下阶乘的逆元
qq_41081096的博客
07-24 414
fac[0] = 1; for(int i = 1; i &lt;= MAX; i++) fac[i] = (fac[i - 1] * i) % MOD; inv_fac[MAX] = qpow(fac[MAX], MOD - 2); for(int i = MAX - 1; i &gt;= 0; i--) inv_fac[i] = (inv_fac[i + 1] * (i + ...
阶乘逆元+巧妙解法
qq_36921652的博客
06-05 1175
B - RGB ColoringTime limit時間制限 : 2sec / Memory limitメモリ制限 : 1024MB配点 : 700 点問題文高橋君はタワーを 1 つ持っており、それは N 個のブロックが縦一列に重なって構成されています。 はじめすべてのブロックは無色ですが、高橋君はいくつかのブロックを赤色、緑色、青色のいずれかの色で塗ることで、 タワーを美しくしようとしています。...
阶乘逆元 记一下
Tc的专栏
04-20 1380
fac[0] = 1; for(int i = 1; i     fac[i] = (fac[i - 1] * i) % MOD; inv_fac[MAX] = qpow(fac[MAX], MOD - 2); for(int i = MAX - 1; i >= 0; i--)     inv_fac[i] = (inv_fac[i + 1] * (i + 1)) % MOD;
hdu 5651 (组合数学 + 阶乘求逆元)
__小明003的博客
03-28 1416
xiaoxin juju needs help Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 1149    Accepted Submission(s): 329 Problem Description
阶乘逆元
hpuspring的博客
08-20 1143
有时候做题会遇到，MOD一般为1e9+7;   代码： #include typedef long long LL; const LL MOD = 1e9 + 7; LL fac[1000000+5]; //阶乘 LL inv[1000000+5]; //逆元 LL quickMod(LL a,LL b) { LL ans = 1; while (b) { if (b&1)
阶乘逆元线推
在AC与WA间徘徊
11-09 450
今天因为昨天晚上左神出了道数论题，虽然是大神界的水题，但让学生们个个爆了零，从中吸取了一个教训，数据大时相乘容易爆int涨了一个姿势，一种求阶乘逆元的打表新方法。
乘法逆元（inverse element）及四大相关求法详解（含证明）
热门推荐
qq_63786973的博客
07-30 1万+
若整数 a，m 互质，并且对于任意的整数 b，如果满足 a | b ( a 能整除 b )，则存在一个整数 x，使得 b/a ≡ bx ( mod m ) ，则称 x 为 a 模 m 的乘法逆元，记为 a^{−1}( mod m ) 或者 inv（ a )。
算法笔记——递推求阶乘逆元（倒推）&线性求逆元
kai_wei_的博客
03-21 1380
求阶乘的逆元，需要求每个阶乘的逆元，一般求最大的逆元，然后倒推剩下的全部。 void init() { fact[0] = 1;//0的阶乘 for (int i = 1; i < maxn; i++) { fact[i] = fact[i - 1] * i %mod; } inv[maxn - 1] = power(fact[maxn - 1], mod - 2);//快速幂 for (int i = maxn - 2; i >= 0; i--) { inv[i] = inv
线性递推阶乘逆元
陈颜的博客
12-31 4378
在求组合数时，我们需要用到公式(n)!(n−m)!⋅(m)!\frac{(n)!}{(n-m)!·(m)!}(n−m)!⋅(m)!(n)!​，其中阶乘可以OnOnOn的递推预处理，除法取模需要用到逆元，逆元又可以用费马小定理和exgcdexgcdexgcd求，不过两者的时间复杂度都是OlognOlognOlogn，事实上，我们可以先求出最后一个阶乘的逆元，再根据公式inv(n!)=inv((n+1)!×(n+1))inv(n!)=inv((n+1)!\times (n+1))inv(n!)=inv((n+1
[xdoj1158]阶乘求逆元(常用于求组合数)
weixin_30294021的博客
05-04 271
http://acm.xidian.edu.cn/problem.php?id=1158 解题关键：此题注意将$\sum\limits_{i = 0}^x {C_x^iC_y^{i + k}}$转化为$C_{x + y}^{x + k}$ 利用二项式定理， 一方面， ${(1 + a)^y}{(1 + \frac{1}{a})^x}$的${a^k}$项的系数，第一个二项式的${a^...
如何O(n)预处理逆元
最新发布
08-16
### 逆元预处理算法的实现方法 在编程中，若需频繁计算 `1, 2, ..., n` 模 `p` 的逆元，可以使用 **O(n)** 时间复杂度的逆元预处理算法。该方法基于以下递推公式： $$ \text{inv}[i] = (p - p / i) \times \text{inv}[p \% i] \% p $$ 其中 `p` 是一个质数，且 `inv[1] = 1`。通过这个公式，可以递推计算出 `1` 到 `n` 的所有逆元。 #### 算法实现步骤 1. 初始化 `inv[0] = 1`（虽然 `inv[0]` 无实际意义，但为了方便计算）。 2. 设置 `inv[1] = 1`，因为 `1` 的逆元是它本身。 3. 使用上述递推公式，从 `2` 到 `n` 依次计算每个数的逆元。 #### 示例代码 ```cpp #include <bits/stdc++.h> using namespace std; const int MOD = 1e9 + 7; // 假设模数为一个大质数 const int MAXN = 1e5 + 10; // 最大需要计算的逆元数量 int inv[MAXN]; // 存储逆元 void precompute_inv(int n) { inv[1] = 1; // 1的逆元是1 for (int i = 2; i <= n; ++i) { inv[i] = (long long)(MOD - MOD / i) * inv[MOD % i] % MOD; } } int main() { int n = 10; // 假设需要计算1~10的逆元 precompute_inv(n); for (int i = 1; i <= n; ++i) { cout << "inv[" << i << "] = " << inv[i] << endl; } return 0; } ``` #### 代码说明 - `MOD` 是模数，通常是一个大质数。 - `precompute_inv(n)` 函数用于预处理 `1` 到 `n` 的逆元。 - `inv[i]` 表示 `i` 的逆元。 #### 算法复杂度分析 - **时间复杂度**：`O(n)`，因为每个逆元的计算只需要一次递推。 - **空间复杂度**：`O(n)`，用于存储 `1` 到 `n` 的逆元。 #### 应用场景 - **组合数计算**：在需要频繁计算组合数 `C(n, k)` 时，可以通过预处理阶乘和阶乘逆元，快速计算出结果[^4]。 - **模运算优化**：在需要多次进行模运算的场景中，预先计算好逆元可以显著提升性能。 #### 注意事项 - 该算法要求模数 `p` 必须是质数。 - 确保 `i` 与 `p` 互质，否则逆元不存在。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值