[学习笔记]乘法逆元

逆元是什么以及为什么要使用逆元

逆元的作用：在模意义下做乘法的逆运算。

在实数运算中，除以一个数，等于乘上这个数的倒数。而现在在模意义下，逆元就充当了“倒数”的角色，模意义下乘上逆元就相当于除以了一个数。

逆元定义：如果一个线性同余方程 $ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$，则 $x$ 称为 $a\mathrm{m}\mathrm{o}\mathrm{d}b$ 的逆元，记作 $a^{-1}$。所以求逆元实际上就是求如下方程的解：
$$ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

费马小定理求逆元

费马小定理

若 $p$ 为素数，且 $\gcd (a,p)=1$，即 $a,p$ 互质，那么 $a^{p-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$。

---

如何用于求逆元

当 $p$ 是质数时，

$\because a^{p-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

$\therefore a\cdot a^{p-2}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

$\therefore x=a^{p-2}\mathrm{m}\mathrm{o}\mathrm{d}p$

用个快速幂就结束了。

注意限制条件：$p$ 必须是质数。

但是大多数情况下OI题目中的模数都是质数，所以这一方法很常用。

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

欧拉定理求逆元

例题：NOIP提高组2012 同余方程

欧拉函数

欧拉函数 $\phi (x)$ 的定义：小于等于 $x$ 的数中与 $x$ 互质的数的个数。

如何求欧拉函数：
$$\phi (x)=(\prod _{i=1}^k\frac{p_i-1}{p_i})x$$

其中 $k$ 是 $x$ 中的质因子个数，$p_i$ 是 $x$ 的一个质因子。注意，这里忽略质因数分解中的次数，如
$$20=2^2\times 5$$

但是我们计算时只考虑 $2$ 和 $5$ 这两个质因子。

---

欧拉定理

若 $\gcd (a,p)=1$，那么 $a^{\phi (p)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$。

---

如何用于求逆元

$\because a^{\phi (p)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

$\therefore a\cdot a^{\phi (p)-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

$\therefore x=a^{\phi (p)-1}\mathrm{m}\mathrm{o}\mathrm{d}p$

质因数分解求 $\phi (p)$ 加上一个快速幂即可。

最坏时间复杂度 $O(\sqrt{p}+\log p)$。

---

和费马小定理的关系

有没有发现二者长得很像？
其实费马小定理是欧拉定理的一种特殊情况：
当 $p$ 为质数时，易得 $\phi (p)=p-1$，此时欧拉定理就成了费马小定理。

---

例题代码

#include<iostream>
using namespace std;
long long a, b;

long long phi(long long x)
{
	long long res = x;
	for(long long i = 2; i * i <= x; i++)
		if(x % i == 0)
		{
			res = res / i * (i - 1);
			while(x % i == 0)
				x /= i;
		}
	if(x > 1)
		res = res / x * (x - 1);
	return res;
}

long long qpow(long long a, long long b, long long c)
{
	long long res = 1;
	while(b)
	{
		if(b & 1)
			res = res * a % c;
		b >>= 1;
		a = a * a % c;
	}
	return res;
}

int main()
{
	cin >> a >> b;
	cout << qpow(a, phi(b) - 1, b) << endl;
	return 0;
}

扩展欧几里得算法求逆元

例题：NOIP提高组2012 同余方程

裴蜀定理

对于任意整数 $a,b$，若 $a,b$ 不全为 $0$，则存在整数 $x,y$，满足 $ax+by=\gcd (a,b)$。

扩展欧几里得算法（简称扩欧，英文名称 exgcd）就是用来求一组满足条件的 $x,y$ 的。

---

扩展欧几里得算法

当 $b=0$ 时，有 $ax=\gcd (a,0)=a$，易得 $x=1,y=0$。

我们知道 $\gcd (a,b)=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$，那么我们在知道

$$b{x}^{\prime }+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y^{\prime }=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$$

的一组特解 $x^{\prime },y^{\prime }$ 的情况下如何推出 $x,y$ 呢？

$$a=\lfloor \frac{a}{b}\rfloor \cdot b+a\mathrm{m}\mathrm{o}\mathrm{d}b$$

将其代入 $ax+by=\gcd (a,b)$ 中整理可得

$$b(\lfloor \frac{a}{b}\rfloor \cdot x+y)+(a\mathrm{m}\mathrm{o}\mathrm{d}b)x=\gcd (a,b)$$

所以

$$\{\begin{array}{l}\lfloor \frac{a}{b}\rfloor \cdot x+y=x^{\prime }\\ x=y^{\prime }\end{array}$$

所以

$$\{\begin{array}{l}x=y^{\prime }\\ y=x^{\prime }-\lfloor \frac{a}{b}\rfloor \cdot x\end{array}$$

---

如何用于求逆元

求方程 $ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 的解本质上就是求二元一次不定方程

$$ax+py=1$$

的解中的整数解 $x$。使用扩展欧几里得算法即可求出 $x$ 的一个特解。

而我们发现，

$$a(x+kp)+p(y-ka)=1$$

所以我们可以不断调整 $x$ 得到通解。

也可看出，逆元存在的充分必要条件是 $a$ 与 $p$ 互质。

时间复杂度 $O(\log (a+p))$。

---

例题代码

注意要求 $x$ 的最小正整数解。

#include<iostream>
using namespace std;

int exgcd(int a, int b, int &x, int &y)
{
	if(b == 0)
	{
		x = 1;
		y = 0;
		return a;
	}
	int d = exgcd(b, a % b, y, x);
	y -= a / b * x;
	return d;
}

int main()
{
	int a, b, x, y;
	cin >> a >> b;
	exgcd(a, b, x, y);
	x %= b;
	if(x < 0)
		x += b;  //注意x可能为负数
	cout << x << endl;
	return 0;
}

线性求连续正整数的逆元

例题：Luogu P3811 【模板】乘法逆元

首先易得 $1^{-1}=1$。

然后因为 $p=\lfloor \frac{p}{i}\rfloor \cdot i+p\mathrm{m}\mathrm{o}\mathrm{d}i$，令 $k=\lfloor \frac{p}{i}\rfloor$，$j=p\mathrm{m}\mathrm{o}\mathrm{d}i$，放到模意义下就会得到
$$ki+j\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

两边同时乘 $i^{-1}{j}^{-1}$ 得到
$$k{j}^{-1}+i^{-1}\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

所以有
$$i^{-1}\equiv -k{j}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

即
$$i^{-1}\equiv -\lfloor \frac{p}{i}\rfloor \cdot (p\mathrm{m}\mathrm{o}\mathrm{d}i{)}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

代码实现（来自OI-wiki）：

inv[1] = 1;
for(int i = 2; i <= n; i++)
    inv[i] = (p - p / i) * inv[p % i] % p;

使用 p - p / i 来防止出现负数，且对结果不会产生影响。

线性求阶乘逆元

$$\because n!=(n-1)!\times n$$

$$\therefore (n!{)}^{-1}\equiv [(n-1)!{]}^{-1}\times n^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

$$\therefore (n!{)}^{-1}\times n\equiv [(n-1)!{]}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

最终我们得到
$$[(n-1)!{]}^{-1}\equiv (n!{)}^{-1}\times n(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

其实不看推导过程这个也很好理解。而 $(n!{)}^{-1}$ 我们可以用扩欧、欧拉定理或者费马小定理算一下就好了。

例题：

1. SDOI2016 排列计数（这题需要其它的知识，求解的过程中要用到此方法）
2. NOI Online #2 入门组 建设城市（这题也需要其它的知识）