还是逆元之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 版权协议，转载请附上原文出处链接和本声明。

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

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

除草

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

然后发现了可以这样

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 1778
求对阶乘 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
阶乘逆元+巧妙解法
qq_36921652的博客
06-05 1159
B - RGB ColoringTime limit時間制限 : 2sec / Memory limitメモリ制限 : 1024MB配点 : 700 点問題文高橋君はタワーを 1 つ持っており、それは N 個のブロックが縦一列に重なって構成されています。 はじめすべてのブロックは無色ですが、高橋君はいくつかのブロックを赤色、緑色、青色のいずれかの色で塗ることで、 タワーを美しくしようとしています。...
阶乘和阶乘逆元
weixin_30483013的博客
08-08 147
扫盲。今天做题才知道这玩意。。（那你之前是怎么算阶乘的哇。。只会暴力暴力暴力嘛。。。。） 阶乘逆元，原理是费马小定理，有些意思。 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 1959
线性求逆元， 阶乘逆元
阶乘逆元 三种 快速 方法 费马小定理 递推 线性
mrcrack的博客
03-31 2745
费马小定理 #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) { ...
逆元的线性递推求解方法及阶乘逆元
ComingLiu的博客
03-11 4656
我们已经学过了利用扩展欧几里得和费马小定理求解乘法逆元的方法，但是对于这道题，这两种方法都失效，那么我们需要找到一种线性的算法来求解，时间复杂度才合格，如果对于逆元概念存在问题可以看一下第一个链接 可以这样考虑， ...
记一下阶乘的逆元
qq_41081096的博客
07-24 396
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 + ...
hdu 5651 (组合数学 + 阶乘求逆元)
__小明003的博客
03-28 1397
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 1124
有时候做题会遇到，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)
[xdoj1158]阶乘求逆元(常用于求组合数)
weixin_30294021的博客
05-04 250
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^...
求逆元、阶乘逆元、线性求逆元
TIANFENG_123456的博客
09-16 341
目录 什么是逆元 如何求逆元 阶乘逆元 线性求逆元 本文章内，若无特殊说明，数字指的是整数，除法指的是整除。 什么是逆元 我们称\(a\)是\(b\)在模\(p\)情况下的逆元，则有\(a \times b \equiv 1 ( mod\,\,p)\)。 ...
O(n)递推求乘法逆元和阶乘逆元
yhhy666的博客
07-21 326
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​。
逆元在阶乘上的运用
QAQ
03-12 885
what’s 逆元 原题：xiaoxin juju needs help 题意： 给一字符串，问用这个字符串里面的字符可以拼成多少个不同的回文串 解析： 把回文串分成左右两边，即字母x出现次数除2，求排列C 假设左半边有x个字母a，y个字母b，z个字母c，C==(x+y+z)! / x! / y! / z! 普通的做法有个bug，(x+y+z)!当然可以做，但是%mod后再除x!呢？...
求组合数、求逆元、求阶乘 O(n)
凝望的博客
08-15 927
在O(n)的时间内求组合数、求逆元、求阶乘。·。· #include &lt;bits/stdc++.h&gt; #define pi acos(-1) #define se second #define fi first #define fastio ios_base::sync_with_stdio(false);cin.tie(NULL);cout.tie(NULL); using nam...
阶乘逆元 记一下
Tc的专栏
04-20 1351
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;
乘法逆元（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 )。
线性阶乘逆元
northsky0307的博客
03-11 246
ll fac[maxn]; ll inv[maxn]; ll qpow(ll a, ll b){ ll res = 1, base = a; while(b) { if(b & 1) res *= base; b >>= 1; base *= base; if(res >= mod) res %= mod; if(base >
乘法逆元
Daniel
03-02 1208
满足x与y的乘积对z取模等于1时，x和y对于模数z来说互为逆元 即 x*y=1(mod z) 为什么要求逆元呢？ 比如求 (b/a)mod p 时，如果b大到不可直接求出，我们可以先求出在模数p下与a互为逆元的数x，那么 (b/a)mod p 转化成 (b*x)mod p证明略。。。求逆元： 方法一（a和p互质，扩展欧几里算法）： ：a * x + p * y = gcd( a , p )
关于数论乘法逆元及相关知识点
yo_bc的博客
05-10 9901
在求解a/b%m时，可以转化为(a%(b*m))/b，转化过程如下 令k = (a/b)/m(向下取整), x = (a/b)%m; a/b = k*m + x (x a = k*b*m + b*x; a%(b*m) = b*x; a%(b*m)/b = x; 得证: a/b%m = a%(b*m)/b; 所有的除法取模问题都可以用这种方法，但是当b很大的时候，则会出现爆精度问题，所以引出乘法逆元，将除法取模转换为乘法取模。 b存在乘法逆元的充要条件是b与模数m互质。设令c为b的逆元，即b*c ≡ 1(
如何在C++中用O(n)的时间复杂度预处理排列组合（不要用乘法逆元）
最新发布
07-19
因此，在模运算下，不使用逆元，且要求O(n)预处理（这里n是指组合数的参数n，即我们要预处理C(n,0)到C(n,n)）的方法，实际上没有标准的O(n)方法（因为除法需要逆元，而不用逆元则无法在模意义下做除法）。...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值