专题·线性同余方程相关各项定理及操作【including 中国剩余定理，扩展中国剩余定理，费马小定理，欧拉定理，乘法逆元

最新推荐文章于 2023-06-04 16:44:42 发布

樱狸❀

最新推荐文章于 2023-06-04 16:44:42 发布

阅读量511

点赞数

分类专栏：数论文章标签：基础数论中国剩余定理 &amp;amp;amp; 扩展CRT 费马小定理欧拉定理线性求解逆元

本文链接：https://blog.csdn.net/qq_43326267/article/details/89425106

版权

数论专栏收录该内容

47 篇文章 2 订阅

订阅专栏

初见安~本专题是基于读者对于线性同余方程的理解~建议先走走：传送门

这里是关于数论的专题（4）【详见数论专栏：）

1.中国剩余定理（CRT）

在n1， n2……nk两两互质的情况下，求解如下方程：

$\large \left\{\begin{matrix} x \equiv a1(mod n_{1})& & & & & & & & & \\ x \equiv a2(mod n_{2}) & & & & & & & & & \\ x \equiv a3(mod n_{3})& & & & & & & & & \\ ...... & & & & & & & & & \\ x \equiv ak(mod n_{k})& & & & & & & & & \end{matrix}\right.$

我们就需要用到中国剩余定理。

定理内容及求解过程如下：

设 $\large N = n_{1} * n_{2} * n_{3} * ...*n_{k}$ 。并且有 $\large m_{1} = N / n _{1}, m_{2} = N / n){2}, ...,m_{k} = N / n_{k};$

这时我们就可以得到一个新的同余方程：

$\large \left\{\begin{matrix} m_{1}x_{1} \equiv 1(mod n_{1})& & & & & & & & & \\ m_{2}x_{2} \equiv 1(mod n_{2})& & & & & & & & & \\ m_{3}x_{3} \equiv 1(mod n_{3})& & & & & & & & & \\ ...... & & & & & & & & & \\ m_{k}x_{k} \equiv 1(mod n_{k})& & & & & & & & & \end{matrix}\right.$

由此我们也可以解出x1，x2……xk。

所以总的方程解为：

$\large x = \sum_{i = 1}^{n} a_{i}m_{i}x_{i} + cN(c \in Z)$

这个定理是很好用的，至于解最后为什么是这么个形式，也是在反复利用gcd的特性，可以手动算一下试试：）（~~说白了就是自己证不来~~

*扩展中国剩余定理（exCRT）

扩展crt就是处理当n不一定为质数的情况，两两合并成x equiv r(mod lcm(n1, n2))，继续联立即可。【lcm：最小公倍数

$\large x = k_1 * n_1 + a_1 = k_2 * n_2 + a_2$ ，用扩欧定理可以解出k1和k2，带入即可。

举个例子：

$\large \left\{\begin{matrix} x \equiv 2 mod 6 & & & & & & & & & \\ x \equiv 4 mod 10& & & & & & & & & \\ x \equiv 8 mod 12& & & & & & & & & \end{matrix}\right.$

对于前两个方程，有2 + 6 * x1 = 4 + 10 * x2，解得x1 = 2, x2 = 1，所以方程可以合并为：

$\large x \equiv 2 + 6 * 2(mod lcm(6,10))$ ，即 $\large x \equiv 14 (mod 30)$

再与后面的方程合并求解即可。最后扩欧出来的就是答案。

在CRT上加几步操作即是了。

2.费马小定理

若p为一素数，且a、p互质（a 不等于p），则有： $\large a^{p - 1} \equiv 1(mod p)$ .。这个定理在求原根（传送门建设中）的时候会用到。

我们可以举个例子看看：求 $\large 3^{2046} \equiv ?(mod 5)$ ，我们可以通过定理 $\large 3^{4} \equiv 1(mod 5)$ 转化为：

$\large 3 ^{2046} = 3^{4 * 511 + 2} \equiv 3^{2} \equiv 4(mod 5)$

所以在应用时我们还可以用到一个结论： $\large a^{b} mod p \equiv a^{b mod(p - 1)} mod p$

3.欧拉定理

费马小定理是欧拉定理的一个分支，也就是p为质数的情况。欧拉定理中：

如果a、n为正整数且互质，则有 $a^{\phi(n)} \equiv 1(mod n)$ 。在费马小定理中，φ（n） = n - 1。

所以欧拉定理总结起来有以下三种情况：

a、n互质， $a^b \equiv a^{b mod\phi(n)}(modn)$
a、n不互质，并且 $b>=\phi(n)$ ， $a^b \equiv a^{b mod \phi(n)+\phi(n)}(modn)$
a、n不互质，并且 $b<\phi(n)$ ，直接暴力求即可。

对于第二种情况的证明，建议多多参考别的博文……因为事实证明第二种和第三种是不能混用一个公式的，但是证着证着却又是可以共用的了。这里就不多说了。不过如果只存在第一种和第二种情况的话，是可以直接套用第二种的公式的。：）古代猪文就是一个很好的例子。

4.乘法逆元

当有a、n互质时，若有 $\large ax \equiv 1(mod n)$ ，则x是a关于n的逆元，记作 $\large a^{-1}$ 。注意，这里的-1并不是-1次方，而只是逆元的符号。作为-1次方时， $\large a * a^{-1} = 1$ ，但是作为逆元，是 $\large a * a^{-1} \equiv 1(mod n)$ ，同余和等式的区别还是很大的：）

有一个结论：在（0， n ] 的范围内，逆元是唯一确定的。这个定理手算一下大概也都能明白。而设a在（0， n ] 内的逆元为x，a的逆元因为有乘上a有同余的性质，所以有 $\large a * x \equiv a * (x + n) \equiv 1(mod n)$ 。这样就可以表示出所有的逆元了。

那么如何求解逆元呢？有以下几种方法：

1）和求解线性同余方程一样，转化为求解一元二次方程ax + ny = 1，用扩展欧几里得定理求解。

2）用到欧拉定理可以得出： $\large \left\{\begin{matrix} ax \equiv 1(mod n) & & & & & & & & & \\ a^{\phi(n)} \equiv 1(mod n) & & & & & & & & & \end{matrix}\right.$ 两个基本式子，得到： $\large x = a^{\phi(n) - 1}$ 。

3）这样求逆元的话如果遇上数据量大了就很容易TLE，因为计算欧拉函数需要耗费大量时间。那么如何在线性时间里求出逆元呢？

以这么一个例题为背景——求解1~N中每个数的逆元。（mod p）

很明显N大一点儿时间复杂度就爆了。下面开始推导：

由取模这个操作的定义取余我们可以尝试这么来表示一下：

设 $\large k = \left \lfloor \frac{p}{a} \right \rfloor$ （向下取整）， $\large r = p mod a$ ，则可得 $\large p = k * a + r, r \in (0, a)$ 。

转化为同余方程可得： $\large k * a + r \equiv 0(mod p)$ 。

所以在此基础上再乘上任何数都可以得到0，我们就可以构造： $\large (k * a + r) * a^{-1} * r^{-1} \equiv 0(mod p)$

由乘法分配律展开： $\large k * r^{-1} + a^{-1} \equiv 0(mod p)$

所以再替换回去可得： $\large a ^{-1} \equiv -k *r^{-1}(mod p)= (p - \left \lfloor \frac{p}{a} \right \rfloor) * (p mod a) ^{-1} mod p$

这就是线性求解的递推转移方程啦！代码实现如下：

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

4）其实更多的时候，我们用到逆元都是用来求组合数，也就是某个数的阶乘的逆元。对于这个，我们有个更直观的求解方法：

ll fac[maxn], inv[maxn];
fac[0] = inv[0] = 1;
for(int i = 1; i <=maxn; i++) fac[i] = fac[i - 1] * i % mod;
inv[maxn] = pow(fac[maxn], mod - 2);
for(int i = maxn - 1; i > 0; i--) inv[i] = inv[i + 1] * (i + 1) % mod;

很明显，因为inv[i]相当于是除以i的阶乘，所以就是 $\frac{i+1}{(i+1)!}=\frac{1}{i!}$ 。

第四行对于inv[maxn]的求法，其实是根据2）得到的。

以上就是这次讲解的全部内容啦！！！

迎评：）
——End——

樱狸❀

关注

0
点赞
踩
2

收藏

觉得还不错? 一键收藏
0
评论
专题·线性同余方程相关各项定理及操作【including 中国剩余定理，扩展中国剩余定理，费马小定理，欧拉定理，乘法逆元

初见安~本专题是基于读者对于线性同余方程的理解~建议先走走：传送门这里是关于数论的专题（4）【详见数论专栏：）1.中国剩余定理（CRT）在n1， n2……nk两两互质的情况下，求解如下方程：我们就需要用到中国剩余定理。定理内容及求解过程如下：设。并且有这时我们就可以得到一个新的同余方程：由此我们也可以解出x1，x2……xk。所以总的方程解为：...
复制链接

扫一扫