乘法逆元

最新推荐文章于 2024-11-12 12:29:51 发布

weifeng2356

最新推荐文章于 2024-11-12 12:29:51 发布

阅读量136

点赞数

分类专栏：算法数论文章标签：逆元算法数论

本文链接：https://blog.csdn.net/qq_45799024/article/details/115267639

版权

算法同时被 2 个专栏收录

15 篇文章 0 订阅

订阅专栏

数论

6 篇文章 0 订阅

订阅专栏

本文介绍了三种在模意义下计算乘法逆元的方法：拓展欧几里得算法、费马小定理结合快速幂以及线性递推。这些方法都在O(logp)的时间复杂度内完成，适用于大规模质数下的计算。同时，提供了具体的C++代码实现，便于理解和应用。

摘要由CSDN通过智能技术生成

定义

如果 $\equiv 1 \pmod p$ ，则 $a$ 和 $b$ 在模 $p$ 意义下互为乘法逆元，记作 $\mathrm{inv} (b)$ .

求法

拓展欧几里得（exgcd）

时间复杂度为 $O(\log p)$

#define ll long long
ll exgcd(ll a, ll b, ll &x, ll &y) {
	if(b == 0) {
		x = 1;
		y = 0;
		return a;
	}
	ll d = exgcd(b, a % b, x, y);
	ll t = x;
	x = y;
	y = t - a / b * x;
	return d;
}
ll inv(ll a, ll p) {
	ll x, y;
	if(exgcd(a, p, x, y) != 1) return -1;
	return ((x % p) + p) % p;
}

费马小定理 + 快速幂

若 $p$ 是质数，且 $\gcd (a, p) = 1$ ，则有 $\equiv 1 \pmod p$ .
即 $\cdot \mathrm{inv}(a) \equiv 1 \equiv a ^ {p-1} \pmod p$ ，于是有 $\mathrm{inv} (a) \equiv a ^ {p-2} \pmod p$ .

时间复杂度为 $O(\log p)$

#define ll long long
ll ksm(ll a,ll b,ll mod){
	ll ans=1;
	for(;b;b>>=1,a=a*a%mod)
		if(b&1) ans=ans*a%mod;
	return ans%mod;
}
ll inv(ll a, ll p) {
	return ksm(a, p - 2, p);
}

线性递推

首先有 $\mathrm{inv} (1) = 1$ .
设 $\lfloor \frac {p} {i} \rfloor = k$ ， $\bmod i = r$ ，即 $\cdot i + r$ .
在模 $p$ 意义下，有 $\cdot i + r \equiv 0 \pmod p$ .
等式两边同乘 $\mathrm{inv} (i) \cdot \mathrm{inv} (r)$ ，得 $\cdot \mathrm{inv} (r) + \mathrm{inv} (i) \equiv 0 \pmod p$ .
移项并化简得 $\mathrm{inv} (i) \equiv - \lfloor \frac {p} {i} \rfloor \cdot \mathrm{inv} (p \bmod i) \pmod p$ .

时间复杂度为 $O (n)$
Luogu - P3811 【模板】乘法逆元

//Luogu P3811 乘法逆元
#include<bits/stdc++.h>
#define N 3000006
#define ll long long
using namespace std;
ll n, p, inv[N];
int main() {
	scanf("%lld%lld", &n, &p);
    inv[1] = 1;
	printf("%lld\n", inv[1]);
	for(ll i = 2; i <= n; ++i) {
		inv[i] = inv[p % i] * (p - p / i) % p;
		printf("%lld\n", inv[i]);
	}
	return 0;
}