数论1.4 同余问题

一、同余概述

1.基本定义和定理

定义1：给定$m\in Z$，若$a\mathrm{m}\mathrm{o}\mathrm{d}m=b\mathrm{m}\mathrm{o}\mathrm{d}m(a,b\in Z)$，
则称 $a$ 和 $b$ 对模 $m$ 同余，记作 $a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$ ，并称该式子为同余式；否则称 $a$ 和 $b$ 对模 $m$ 不同余。

定义2：$Z$ 被分为 $m$ 个不同的集合 $A_{1...m}$，这些集合被称为模 $m$ 剩余类（同余类），如果 $\forall a,b\in A_i,i\in [1,m]$, 则 $a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$。

费马小定理(费小)：若 $p$ 是质数，则对于 $\forall a\in Z$，有 $a^p\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}p)$

证明：对 $a$ 使用数学归纳法。

Base step:当 $a=0$ 时，$0\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$。

Inductive step:假设$a^p\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}p)$，考虑 $(a+1{)}^p$，根据二项 式定理展开，$(a+1{)}^p={\sum }_{k=0}^p{C}_p^k{a}^k$。

由于 $p$ 是质数，$C_p^k=\frac{p!}{k!(p-k)!}$，则$p|C_p^k,k\in [1,p-1]$,则$(a+1{)}^p\equiv a^p+1(\mathrm{m}\mathrm{o}\mathrm{d}p)$。

又因为 $a^p\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}p)$，则 $(a+1{)}^p\equiv a^p+1\equiv a+1(\mathrm{m}\mathrm{o}\mathrm{d}p)$。

根据数学归纳法，对于$\forall a\in Z$，都成立 $a^p\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}p)$，证毕。

欧拉定理：若 $a,n$ 互质，$a,n\in N^{\ast }$，则 $a^{\phi (n)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)$，
其中 $\phi (n)$ 为 欧拉函数。

欧拉定理的推论：若 $a,n$互质，$a,n\in N^{\ast }$，则对于 $\forall b\in N^{\ast }$，有
$a^b\equiv a^{b\mathrm{m}\mathrm{o}\mathrm{d}\phi (n)}(\mathrm{m}\mathrm{o}\mathrm{d}n)$。

证明：

设 $b=q\phi (n)+r$， 其中 $r\in [0,\phi (n))$，$q,r\in N$，即 $r=b\mathrm{m}\mathrm{o}\mathrm{d}\phi (n)$。

于是：$a^b\equiv a^{q\phi (n)+r}\equiv (a^{\phi (n)}{)}^q{a}^r\equiv 1^q{a}^r\equiv a^r\equiv a^{b\mathrm{m}\mathrm{o}\mathrm{d}\phi (n)}(\mathrm{m}\mathrm{o}\mathrm{d}n)$。

证毕。

2.同余的性质

对于 $a,b,c,d,m,n,x,p,q(a,b,c,d,x\in Z,m,n,p,q\in N)$，

对模 $m$ 同余满足：

1. 自反性：$a\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}m)$；
2. 对称性：若 $a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$，则 $b\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}m)$；
3. 传递性：若 $a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$，$b\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}m)$，则 $b\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}m)$；
4. 同加性：若 $a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$，则 $a+c\equiv b+c(\mathrm{m}\mathrm{o}\mathrm{d}m)$；
5. 同乘性：若 $a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$，则 $ac\equiv bc(\mathrm{m}\mathrm{o}\mathrm{d}m)$;
   一般情况下，$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$，$c\equiv d(\mathrm{m}\mathrm{o}\mathrm{d}m)$，则 $ac\equiv bd(\mathrm{m}\mathrm{o}\mathrm{d}m)$；
6. 同幂性*：若 $a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$，则 $a^n\equiv b^n(\mathrm{m}\mathrm{o}\mathrm{d}m)$；
7. 若 $a\equiv x(\mathrm{m}\mathrm{o}\mathrm{d}p)$，$a\equiv x(\mathrm{m}\mathrm{o}\mathrm{d}q)$，其中 $p,q$ 互质，则 $a\equiv x(\mathrm{m}\mathrm{o}\mathrm{d}pq)$；

证明：因为$a\equiv x(\mathrm{m}\mathrm{o}\mathrm{d}p)$，$a\equiv x(\mathrm{m}\mathrm{o}\mathrm{d}q)$，其中 $p,q$ 互质， 则一定 $\exists s,t\in Z$，使 $a=sp+x,a=tq+x$，所以 $sp=tq,q|s,p|t$， 则一定 $\exists r\in Z$，使 $s=rq$， 所以：$a=rpq+x$，得出 $a\equiv x(\mathrm{m}\mathrm{o}\mathrm{d}pq)$。

8. 但是，同余不满足同除性，及不满足 $a/n\equiv b/n(\mathrm{m}\mathrm{o}\mathrm{d}m)$。

二、扩展欧几里得算法（扩欧）

1.欧几里得算法（辗转相除法）

$\forall a,b\in N,b\ne 0$，则$$gcd(a,b)=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$$

2.推导与证明

定理一

设 $a,b\in Z$，则 $\exists x,y\in Z$，使得 $ax+by=gcd(a,b)$。

证明：因为 $gcd(a,b)=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$;

所以 $ax+by=b{x}^{\prime }+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y^{\prime }$；
$\Rightarrow ax+by=b{x}^{\prime }+(a-\lfloor \frac{a}{b}\rfloor b)y^{\prime }$；

$\Rightarrow ax+by=b{x}^{\prime }+a{y}^{\prime }-\lfloor \frac{a}{b}\rfloor b{y}^{\prime }$；

$\Rightarrow ax+by=a{y}^{\prime }+b(x^{\prime }-\lfloor \frac{a}{b}\rfloor y^{\prime })$；

$\Rightarrow x=y^{\prime },y=x^{\prime }-\lfloor \frac{a}{b}\rfloor y^{\prime }$。

当欧几里得算法运行到最后一步时，$b=0$，$gcd(a_n,b_n)=a_n$，

有 $x=1,y=0$。

对递推过程使用数学归纳法，可知扩欧成立。

证毕。

定理二（贝祖定理）

对于不定方程 $ax+by=c$，当且仅当 $gcd(a,b)|c$ 时，方程有整数解。

证明：

已知 $gcd(a,b)|a$，$gcd(a,b)|b$，由整除的性质，$\forall a,b\in Z$，则$gcd(a,b)|ax+by$，即
$gcd(a,b)|c$。

3.求最小整数解

给出不定方程 $ax+by=gcd(a,b)$，已知一组特解 $x_0,y_0$，求最小整数解 $x$。

推导：

因为 $ax+by=gcd(a,b),a{x}_0+b{y}_0=gcd(a,b)$，所以 $ax+by=a{x}_0+b{y}_0$。

上式转化为 $a(x_0-x)=-b(y_0-y)$。

两边同除以 $gcd(a,b)$，得
$\frac{a}{gcd(a,b)}(x_0-x)=-\frac{b}{gcd(a,b)}(y_0-y)$

易得$\frac{a}{gcd(a,b)}$ 与 $\frac{b}{gcd(a,b)}$ 互质，则$\frac{b}{gcd(a,b)}|(x_0-x)$，$\frac{a}{gcd(a,b)}|(y_0-y)$。

则$x=x_0-k\frac{b}{gcd(a,b)},k\in Z$，即 $x=x_0\mathrm{m}\mathrm{o}\mathrm{d}\frac{b}{gcd(a,b)}$。

而此时的 $x$ 有可能是负数，所以还要再加上 $\frac{b}{gcd(a,b)}$ 再 $\mathrm{m}\mathrm{o}\mathrm{d}\frac{b}{gcd(a,b)}$。

则 $x=(x_0\mathrm{m}\mathrm{o}\mathrm{d}\frac{b}{gcd(a,b)}+\frac{b}{gcd(a,b)})\mathrm{m}\mathrm{o}\mathrm{d}\frac{b}{gcd(a,b)}$

4.代码

void exgcd(int a, int b, int &x, int &y)
{
	if(b == 0)
	{
		x = 1; y = 0;
	}
	else
	{
		int x1 = x, y1 = y;
		exgcd(b, a % b, x1, y1);
		x = y1;
		y = x1 - (a / b) * y1;
	}
} 
int gcd(int a, int b)
{
	if(b == 0)	return a;
	else	return gcd(b, a % b);
}
int main()
{
	cin >> a >> b;
	exgcd(a, b, x, y);
	d = gcd(a, b);
	cout << (x % (b / d) + b / d) % (b / d);
}

5.题目

Luogu P4549 【模板】裴蜀定理

Luogu P1082 [NOIP2012 提高组] 同余方程

Luogu P1516 青蛙的约会

三、线性同余方程

定义一

给定 $a,b,c\in Z$，求 $x\in Z$，使得 $ax\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}b)$，或者给出无解。

因为未知数的系数为1，所以我们称之为一次同余方程，也称线性同余方程。

原方程等价于 $ax-c\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}b)$，即 $\exists y\in Z$，使得 $-by=(ax-c)$，即 $ax+by=c$。

因此当 $gcd(a,b)|c$ 时，使用扩欧求解即可。

换而言之，有解的线性同余方程皆可用扩欧求解。

定义二

特别地，当 $gcd(a,b)=1$（即 $a,b$ 互质时），则称 $ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$ 的一个解为 $a$ 模 $b$ 的逆。

四、乘法逆元

给定 $a,b\in Z$ 且 $a,b$ 互质，则 $\exists x\in Z$，使得 $ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$，称 $x$ 为 $a$ 的模 $b$ 乘法逆元，记作 $\frac{1}{a}(\mathrm{m}\mathrm{o}\mathrm{d}b)$ 或 $a^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}b)$。

1.扩欧求逆元

将 $ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$ 转化为 $ax+by=1$ ，因为 $a,b$ 互质，所以 $gcd(a,b)|1$，使用扩欧求解即可。

#include <iostream>
using namespace std;
int m, n, x, y;
void exgcd(int a, int b, int &x,int &y)
{
	if(b == 0)
		x = 1, y = 0;
	else
	{
		int x2, y2;
		exgcd(b, a % b, x2, y2);
		x = y2;
		y = x2- (a / b) * y2;
	}
}

2.快速幂求逆元

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

当b为质数时，由费小得，$a^{b-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$，则 $a\cdot a^{b-2}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$，即 $a^{-1}\equiv a^{b-2}(\mathrm{m}\mathrm{o}\mathrm{d}b)$，使用快速幂求出 $a^{b-2}$ 即可。

int quickpow(int a, int b)
{
	long long p = b - 2, ans = 1;
	while(p != 0)
	{
		if(p & 1)
			ans = ans * a % b;
		p >>= 1;
		a = a * a % b; 
	}
	return ans % b;
} 

3.递推求逆元

求整数 $i=1...M-1$ 模 $M$ 的逆元 $inv[i]$，则$$inv[i]=(M-\lfloor M/i\rfloor )\cdot inv[M\mathrm{m}\mathrm{o}\mathrm{d}i]\mathrm{m}\mathrm{o}\mathrm{d}M$$

推导过程：

设 $t=\lfloor M/i\rfloor$， $k=M\mathrm{m}\mathrm{o}\mathrm{d}i$.

$k\equiv M-ti(\mathrm{m}\mathrm{o}\mathrm{d}M)$

$\Rightarrow k+ti\equiv M(\mathrm{m}\mathrm{o}\mathrm{d}M)$

$\Rightarrow k+ti\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}M)$

$\Rightarrow -ti\equiv k(\mathrm{m}\mathrm{o}\mathrm{d}M)$

两边同时除以 $ik$，得

$-\frac{t}{k}\equiv \frac{1}{i}(\mathrm{m}\mathrm{o}\mathrm{d}M)$

$\Rightarrow -t{k}^{-1}\equiv i^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}M)$

$\Rightarrow -t\cdot inv[k]\equiv inv[i](\mathrm{m}\mathrm{o}\mathrm{d}M)$

把 $t$ 和 $k$ 替换掉，得

$inv[i]=-\lfloor M/i\rfloor \cdot inv[M\mathrm{m}\mathrm{o}\mathrm{d}i]\mathrm{m}\mathrm{o}\mathrm{d}M$

$\Rightarrow inv[i]=(M-\lfloor M/i\rfloor )\cdot inv[M\mathrm{m}\mathrm{o}\mathrm{d}i]\mathrm{m}\mathrm{o}\mathrm{d}M$

void invmod(int m)
{
	inv[1] = 1;
	for(int i = 2; i < m; i++)
		inv[i] = (m - m / i) * inv[m % i] % m;
} 

4.阶乘求逆元

	fact[0] = 1;
	for(re ll i = 1; i <= n; i++)	fact[i] = fact[i - 1] * i % mod;
	inv[n] = exgcd(n, mod, x, y);
	for(re ll i = n - 1; i >= 1; i--)	inv[i] = inv[i + 1] * (i + 1) % mod;

5.题目

Luogu P3811【模板】乘法逆元

五、中国剩余定理（CRT）

1.定理

今有物不知其数，三三数之剩二，五五数之剩三，七七数之剩二，问物几何？
答曰：‘二十三’。
术曰：三三数之剩二，置一百四十；五五数之剩三，置六十三，七七数之剩二，置三十，并之。得二百三十三，以二百一十减之，即得。凡三三数之剩一，则置七十；五五数之剩一，则置二十一；七七数之剩一，则置十五；一百六以上以一百五减之即得。
——《孙子算经》

找出三个数：从3和5的公倍数中找出被7除余1的最小数15，从3和7的公倍数中找出被5除余1 的最小数21，最后从5和7的公倍数中找出除3余1的最小数70。

用15乘以2（2为最终结果除以7的余数），用21乘以3（3为最终结果除以5的余数），同理，用70乘以2（2为最终结果除以3的余数），然后把三个乘积相加15∗2+21∗3+70∗2得到和233。

用233除以3，5，7三个数的最小公倍数105，得到余数23，即233%105=23。这个余数23就是符合条件的最小数。

以上就是孙子定理，国际上称为中国剩余定理（CRT），定理如下：

若 $m_1,m_2\dots m_n$ 是两两互质的正整数，$M={\prod }_{i=1}^n{m}_i$，$M_i=M/m_i$，$t\equiv M_i^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$。对于$\forall a_1,a_2\dots a_n\in Z$，则同余方程组
$$\{\begin{array}{cc}& x\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_1)\\ & x\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{m}_2)\\ & \dots \dots \\ & x\equiv a_n(\mathrm{m}\mathrm{o}\mathrm{d}{m}_n)\end{array}$$
有整数解，方程组的解为 $$x=\sum _{i=1}^n{a}_i{M}_i{t}_i$$并且在模 $M$ 意义下有唯一解。

推导： 设 $x_i\equiv a_i(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$。

由 $a\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}b)\iff a+kb\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}b)$ 可知，要使 $x\equiv a_i(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$，则 $x_i$ 是除
$m_i$ 之外所有模数的倍数。

又因为模数两两互质，所以 $x_i$ 是 $M_i=M/m_i$ 的倍数，$M_i$ 也是除 $m_i$ 以外的模数的最小公倍数。

设 $n_i=\frac{x_i}{M_i}$，则 $n_i{M}_i\equiv a_i(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$，两边同时除以 $a_i$，则
$\frac{n_i}{a_i}Mi\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$，则 $t_i=\frac{x_i}{a_i{M}_i}$。

则 $x_i=a_i{M}_i{t}_i$， $x={\sum }_{i=1}^n{a}_i{M}_i{t}_i$。

2.代码

int Crt()
{
	int ans = 0, M = 1, Mi, x, y;
	for(int i = 1; i <= n; i++)
		M *= m[i];
	for(int i = 1; i <= n; i++)
	{
		Mi = M / m[i];
		exgcd(Mi, m[i], x, y);
		ans = (ans + Mi * a[i] * x) % M;
	}
	return (ans + M) % M;
}

3.题目

Luogu P1495 【模板】中国剩余定理(CRT)/曹冲养猪

六、扩展中国剩余定理 （EXCRT）

1.定理

给定 $n$ 组非负整数 $a_i,m_i$，求解关于 $x$ 的方程组的最小非负整数解。

$$\{\begin{array}{cc}& x\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_1)\\ & x\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{m}_2)\\ & \dots \dots \\ & x\equiv a_n(\mathrm{m}\mathrm{o}\mathrm{d}{m}_n)\end{array}$$

推导：

设 $x_0\equiv a_{i-1}(\mathrm{m}\mathrm{o}\mathrm{d}{m}_{i-1})$，$x_0\equiv a_i(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i),i=2\dots n$,

则 $k{m}_1+a_1=-k^{\prime }{m}_2+a_2,k,k^{\prime }\in Z$，

上式可化为 $k{m}_1+k^{\prime }{m}_2=a_2-a_1$，求得最小正整数解 $k,k^{\prime }$。

将 $k,k^{\prime }$ 代入原式，可求得 $x_0$。

则此时 $x\equiv x_0(\mathrm{m}\mathrm{o}\mathrm{d}lcm(m_{i-1},m_i))$。

将上式与下一式联立，用上述方法持续求解，直到 $i=n$时即可。

2.代码

int excrt()
{
	for(int i = 2; i <= n; i++)
	{
		a = M[i - 1]; b = M[i]; c = A[i] - A[i - 1];
		d = exgcd(a, b, x, y);
		x = (x * (c / d) % (b / d) + (b / d)) % (b / d);
		A[i] = a * b / d;
		B[i] = x * a + B[i - 1];
	}
	return B[n];
}

3.题目

Luogu P4777 【模板】扩展中国剩余定理（EXCRT）