乘法逆元
当求解答案过程中对p取余同时又涉及除法运算时,需要把除法转成乘法才能保证经过模法运算后答案是正确的。
一、线性递推求逆元
输入
n
n
n ,
p
p
p 。
求
1
,
2
,
3
,
.
.
.
,
n
1,2,3,...,n
1,2,3,...,n 在模
p
p
p 意义下的逆元。
inv[1] = 1;
for(int i = 2; i <= n; i++) inv[i] = (p - p/i) * inv[p%i]%p;
二、费马小定理求逆元
输入
a
a
a ,
p
p
p
求
a
a
a 在模
p
p
p 意义下的逆元
逆元为 : a p − 2 a^{p-2} ap−2 / / 通 过 快 速 幂 实 现 //通过快速幂实现 //通过快速幂实现
三、拓展欧几里德求逆元
输入
a
a
a ,
p
p
p
求
a
a
a 在模
p
p
p 意义下的逆元
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
#define ll long long
void extend_gcd(ll a, ll b, ll &x, ll &y) {
if (!b)
{
x = 1, y = 0;
return;
}
else
{
extend_gcd(b, a % b, y, x);
y -= x * (a / b);
return;
}
}
ll inv(ll a, ll n)
{
ll x, y;
extend_gcd(a,n,x,y);
x = (x % n + n) % n;
return x;
}
int main()
{
ll a,p;
scanf("%lld%lld",&a,&p);
printf("%lld\n",inv(a,p));
return 0;
}