P3811 【模板】乘法逆元

线性求逆元

逆元定义：若$a*x\equiv1 (\bmod {b})$，且$a$与$b$互质，那么我们就能定义: $x$为$a$的逆元，记为$a^{-1}$，所以我们也可以称$x$为$a$的倒数，

所以对于$\frac{a}{b} (\bmod {p})$ ，我们就可以求出$b$在$\bmod {p}$下的逆元，然后乘上$a$，再$\bmod {p}$，就是这个乘法逆元的值了。

一、exgcd求逆元(O(l$og_n$))

这个就是利用拓欧求解线性同余方程$a*x \equiv c (\bmod {b})$

的c=1的情况。我们就可以转化为$a*x + b*y = 1$，求解这个方程的解。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>#define ll long long
using namespace std;void exgcd(ll a,ll b,ll &x,ll &y){if(!b){x=1,y=0;return;}exgcd(b,a%b,x,y);ll tmp=x;x=y;y=tmp-a/b*y;
}int n,p;int main()
{scanf("%d%d",&n,&p);for(int i=1;i<=n;i++){ll x,y;exgcd(i,p,x,y);x=(x%p+p)%p;printf("%lld\n",x);}return 0;
}

Exgcd

二、快速幂+费马小定理

费马小定理：若$p$为素数，$a$为正整数，且$a,p$互质。则有$a^{p-1} \equiv 1 (\bmod {p})$。

$a*x\equiv 1(\bmod {b})$

$a*x\equiv a^{p-1} (\bmod {b})$
$x \equiv a^{p-2} (\bmod {b})$

有点儿慢。。。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>#define ll long long
using namespace std;int n,p;ll pow(int a,int b,int mod){ll s=1,t=a%mod;for(;b;b>>=1,t=1ll*t*t%mod)if(b&1) s=1ll*s*t%mod;return s;
}int main()
{scanf("%d%d",&n,&p);for(int i=1;i<=n;i++){printf("%lld\n",pow(i,p-2,p));}return 0;
}

快速幂

三、线性递推求逆元

如何递推呢？

大佬说：推一下就好了嘛

蒟蒻（我）：。。。=_=

首先$1^{-1}\equiv1(modp)$

设$p=k\times i+r$，$r < i$，$1< i < p$

那么$k\times i+r \equiv 0 (mod p)$

方程同乘$i^{-1},r^{-1}$得

$k\times r^{-1}+i^{-1} \equiv 0 (mod p)$

移项得：$i^{-1}\equiv -k\times r^{-1} (mod p)$

即$i^{-1}\equiv\lfloor{\frac{p}{i}}\rfloor \times (p mod i)^{-1} (mod p)$

也就是下面这一行代码：

inv[i]=(p-p/i)*inv[p%i]%p;

奉上代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>#define ll long long
#define N 10101010
using namespace std;int n,p;
ll inv[N];
int main()
{scanf("%d%d",&n,&p);inv[1]=1;printf("1\n");for(int i=2;i<=n;i++){inv[i]=(p-p/i)*inv[p%i]%p;printf("%lld\n",inv[i]);}return 0;
}

转载于:https://www.cnblogs.com/song-/p/9724555.html

洛谷——P3811 【模板】乘法逆元相关推荐

【训练题22：线性求逆元】【模板】乘法逆元 | 洛谷 P3811
[模板]乘法逆元难度普及/提高\color{yellow}普及/提高普及/提高板子题,给出两种方法题意给 n,pn,pn,p 求出所有 i∈[1,n]i\in[1,n]i∈[1,n]在模 p ...
乘法逆元（洛谷-P3811）
题目描述给定n,p求1~n中所有整数在模p意义下的乘法逆元. 输入输出格式输入格式: 一行n,p 1≤n≤3×10^6 ,n<p<20000528,输入保证 p 为质数. 输出格式: ...
专题·树链剖分【including 洛谷·【模板】树链剖分
初见安~~~终于学会了树剖~~~ [兴奋]当初机房的大佬在学树剖的时候我反复强调过:"学树剖没有前途的!!!" 恩.真香. 一.重链与重儿子所谓树剖--树链剖分,就是赋予一个链的 ...
洛谷·【模板】点分树 | 震波【including 点分树
初见安-这里是传送门:洛谷P6329 [模板]点分树 | 震波一.点分树其实你会点分治的话,点分树就是把点分治时的重心提出来重新连城一棵树. 比如当前点是u,求出子树v的重心root后将root与 ...
洛谷 P1919 模板】A*B Problem升级版（FFT快速傅里叶）
https://www.luogu.com.cn/problem/P1919 题目背景本题数据已加强,请使用 FFT/NTT,不要再交 Python 代码浪费评测资源. 题目描述给你两个正整数 a ...
洛谷 P3803 多项式乘法
题目背景这是一道FFT模板题题目描述给定一个n次多项式F(x),和一个m次多项式G(x). 请求出F(x)和G(x)的卷积. 输入输出格式输入格式: 第一行2个正整数n,m. 接下来一行n+1 ...
洛谷.4897.[模板]最小割树(Dinic)
题目链接最小割树模板.具体见:https://www.cnblogs.com/SovietPower/p/9734013.html. ISAP不知为啥T成0分了.. Dinic: //1566ms ...
强连通分量：洛谷P3387 模板：缩点
传送门顾名思义,模板awa #include <cstdio> #include <cstring> #include <cmath> #include < ...
【后缀数组】洛谷P3809模板题
题目背景这是一道模板题. 题目描述读入一个长度为 n n n 的由大小写英文字母或数字组成的字符串,请把这个字符串的所有非空后缀按字典序从小到大排序,然后按顺序输出后缀的第一个字符在原串中的位置. ...

洛谷——P3811 【模板】乘法逆元

P3811 【模板】乘法逆元

一、exgcd求逆元(O(l$og_n$))

二、快速幂+费马小定理

三、线性递推求逆元

洛谷——P3811 【模板】乘法逆元相关推荐

最新文章

热门文章