扩展欧几里得、中国剩余定理和扩展中国剩余定理

1 扩展欧几里得算法

扩展欧几里得算法主要可以解决：
·求乘法逆元
·求解同余方程/二元一次方程
其中,求乘法逆元可以转化成求解同余方程。

1.1求解方程$ax+by=gcd(a,b)$的一组特解

先考虑以下二元一次方程：

$ax+by=gcd(a,b)(1)$

$b{x}^{\prime }+(a\%b)y^{\prime }=gcd(b,a\%b)=gcd(a,b)(2)$

$...$

$gcd(a,b)x^{\prime \prime }+0\ast y^{\prime \prime }=gcd(a,b)(3)$

对于方程$(3)$，显然它的一个解是$x^{\prime \prime }=1,y^{\prime \prime }=0;$

想要求解方程$(1)$,我们就要找到方程$(1)(2)$中$x,y$与$x^{\prime },y^{\prime }$之间的关系。

将方程$(2)$变形:

$b{x}^{\prime }+(a\%b)y^{\prime }=gcd(a,b)$

$b{x}^{\prime }+(a-\lfloor a/b\rfloor b)y^{\prime }=gcd(a,b)$

将含$a$和$b$的项合并，得到

$a{y}^{\prime }+b[(x^{\prime }-\lfloor a/b\rfloor )y^{\prime }]=gcd(a,b)(4)$

将其与方程$(1)$比较：

$ax+by=gcd(a,b)(1)$

容易看出，可以取

$x=y^{\prime }$

$y=x^{\prime }-\lfloor a/b\rfloor y^{\prime }$

作为方程$(1)$的解。

因此，可以在递归调用欧几里得算法时不断更新$x$和$y$，得到方程的一组解。

$Code:$

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

1.2求出二元一次方程$ax+by=gcd(a,b)$的通解和最小正整数解

$<$结论$>$设方程$ax+by=gcd(a,b)$有特解$x=x_0,y=y_0$,

则方程$ax+by=gcd(a,b)$通解形如

$x=x_0+t\ast (b/gcd(a,b))$

$y=y_0-t\ast (a/gcd(a,b))(t\in \mathbb{Z})$,

最小正整数解$x_{min}=(x_0+p)\%p.$

$<$证明$>$在$(1.1)$中，我们解出了$ax+by=gcd(a,b)$的一组通解$(x_0,y_0)$，设该方程的通解形如

$x=x_0+m$

$y=y_0+n$

代入方程得

$a(x_0+m)+b(y_0+n)=gcd(a,b)=a{x}_0+b{y}_0$

$a{x}_0+am+b{y}_0+bn=a{x}_0+b{y}_0$

$am+bn=0$

$a/b=-(n/m)=(a/gcd(a,b))/(b/(gcd(a,b)))$

由于$a/gcd(a,b)$和$b/gcd(a,b)$一定互质(分数为最简分数)，那么$m$和$n$一定是在这个分数基础上分子分母同乘一个整数$t$得到的，因此一定有

$m=t\ast (b/gcd(a,b))$

$n=(-t)\ast (a/gcd(a,b))$

因此方程通解

$x=x_0+t\ast (b/gcd(a,b))$

$y=y_0-t\ast (a/gcd(a,b))$

$t\in \mathbb{Z}.$

若要求出最小正整数的$x_{min}$,记$b/gcd(a,b)=p$,

当$x\ge 0$时，$x_{min}=(x_0+tp)\%p=x\%p$;

当$x<0$时，$x_{min}=(x_0+p)\%p.$

1.3求出二元一次方程$ax+by=c$的特解

$<$结论$>$若方程$ax+by=gcd(a,b)$有特解$x=x_0,y=y_0$,

则方程$ax+by=c$当$c$能被$gcd(a,b)$整除时有特解

$x=k{x}_0$

$y=k{y}_0,$

其中$k=c/gcd(a,b)$.

$<$证明$>$通过上面$(1.1)$的算法可以求解方程$ax+by=gcd(a,b)$的一组解$(x_0,y_0)$，即

$a{x}_0+b{y}_0=gcd(a,b)$

设$c/gcd(a,b)=k$(仅当$c$被$gcd(a,b)$整除，方程有解)，于是

$k(a{x}_0+b{y}_0)=k\cdot gcd(a,b)$

$ka{x}_0+kb{y}_0=c$

$a(k{x}_0)+b(k{y}_0)=c$

因此，$x=k{x}_0,y=k{y}_0$是方程$ax+by=c$的一组特解。

1.4求出二元一次方程$ax+by=c$的通解

$<$结论$>$若方程$ax+by=c$有特解$x=x_1,y=y_1$,

则方程$ax+by=c$有通解

$x=x_1+m\ast (b/gcd(a,b))$

$y=y_1-m\ast (a/gcd(a,b)),(m\in \mathbb{Z}).$

$<$证明$>$记$ax+by=gcd(a,b)$的一组特解为$x_1,y_1$,那么由$(1.3)$中的推导可知，$x=k{x}_1,y=k{y}_1$是方程$ax+by=c$的解。

假设方程有特解形如$x^{\prime }=k{x}_1+\Delta x,y^{\prime }=k{y}_1+\Delta y$,

就有

$a{x}^{\prime }+b{y}^{\prime }=a(k{x}_1+\Delta x)+b(k{y}_1+\Delta y)=c$

$ak{x}_1+bk{y}_1+a\Delta x+b\Delta y=c=ak{x}_1+bk{y}_1$

即$a\Delta x+b\Delta y=0.$

类似于$(1.2)$中的过程，可得

$a/b=-(\Delta y/\Delta x)=(a/gcd(a,b))/(b/gcd(a,b))$

因此

$\Delta x=m\ast (b/gcd(a,b))$

$\Delta y=-m\ast (a/gcd(a,b))$

所以$ax+by=c$的通解

$x=k{x}_1+\Delta x=k{x}_1+m\ast (b/gcd(a,b))$

$y=k{y}_1+\Delta y=k{x}_2-m\ast (a/gcd(a,b))$

$m\in \mathbb{Z}.$

1.5求乘法逆元

称一个数$a$在模$p(gcd(a,p)=1)$意义下的乘法逆元$x$为满足条件
$ax\equiv 1(modp)$的一个$x.$

由同余的意义，可以将这个条件转化为

$ax=kp+1(k\in \mathbb{Z})$

即$ax-kp=1(gcd(a,p)=1)$

这时就转化为$(1.1)$中方程$ax+by=gcd(a,b)$当$gcd(a,b)=1$时的情况，就可以用扩展欧几里得算法求解了.然而求出的一个解$x_0$有可能为负值，如求解$5$在模$13$意义下的乘法逆元：

$5x-13k=1$

这时解$x=-5,k=-2$也成立。

而正确的乘法逆元为$8$.此时对$x$加上模数$13$即可。

$<$证明$>$若$a{x}_0-k_{0}p=1,-p<x_0<0(a,x_0,k_0,p\in \mathbb{Z})$则$x^{\prime }=x_0+p$是一个乘法逆元。

要证明$x^{\prime }$是一个$a$的乘法逆元，只需要证明

$a{x}^{\prime }\equiv 1(modp)$

而$a{x}^{\prime }=a(x_0+p)=a{x}_0+ap=k_{0}p+1+ap=(k_0+a)p+1$

即$a{x}^{\prime }=(k_0+a)p+1$,两边同时取模得

$a{x}^{\prime }\equiv 1(modp).\square$

为了与$x\ge 0$的情况统一，最终对$x$进行$x=(x+p)\%p$可确保得到一个$[0,p-1]$的逆元。

$Code:$

int inv(int a,int mod)
{
	int x,y;
	exgcd(a,mod,x,y);
	return (x+mod)%mod;
}

2 中国剩余定理(CRT)

《孙子算经》中记载了一个问题：“有物不知其数，三三数之剩二，五五数之剩三，七七数之剩二。问物几何？”转化为数学语言为：

$x\equiv 2(mod3)(1)$

$x\equiv 3(mod5)(2)$

$x\equiv 2(mod7)(3)$

求解$x$。利用中国剩余定理可以解决这种问题模数之间互质的情况。(不互质可以用扩展中国剩余定理求解)

设上面三个方程的解分别为$x_1,x_2,x_3$,先考虑使$X=x_1+x_2+x_3$对于方程$(1)$成立的条件：

若$x_2,x_3$均是$3$的倍数，那么

$X\%3=(x_1+x_2+x_3)\%3=(x_1+3k_2+3k_3)\%3=x_1\%3=2.$

这样，我们就找到了一个使$x_1+x_2+x_3$对方程$(1)$成立的条件。

同时用方程$(2)(3)$来限制$X=x_1+x_2+x_3$，可以得到:

$x_1$是$5$和$7$的倍数;

$x_2$是$3$和$7$的倍数;

$x_3$是$3$和$5$的倍数.

那么，如果我们分别找到满足这些条件的一组$x_1,x_2,x_3$，由上面的推导可知，$X=x_1+x_2+x_3$就是方程组$(1)(2)(3)$的一个解。

这时利用乘法逆元，有：

$x\cdot inv(x,p)\equiv 1(modp)$

其中$inv(x,p)$为$x$在模$p$意义下的乘法逆元。

代入$x=5\ast 7=35,$得：

$35\ast inv(35,3)\equiv 1(mod3)$

等式两边同时乘以方程$(1)$的余数$2$,得：

$35\ast inv(35,3)\ast 2\equiv 2(mod3)$.

取$x_1=35\ast inv(35,3)\ast 2$即可，$x_2,x_3$可用类似方法求出。

我们就取得了一个$X=x_1+x_2+x_3$满足方程组$(1)(2)(3).$

然而这个$X$不一定是满足方程组最小的$x$，如果要找出最小的$x$只需要把$X$对$3\ast 5\ast 7=105$取余即可。

经过这个例子，我们可以总结出更一般的求解算法：

对于同余方程组

$x\equiv b_1(mod{a}_1)$

$x\equiv b_2(mod{a}_2)$

$...$

$x\equiv b_n(mod{a}_n)$

其中$a_1,a_2,...,a_n$互质。

记$PI={\prod }_{i=1}^n{a}_i=a_1{a}_2...a_n$,

那么可以取

$x_1=(PI/a_1)\ast inv(PI/a_1,b_1)\ast b_1$

$x_2=(PI/a_2)\ast inv(PI/a_2,b_2)\ast b_2$

$...$

$x_n=(PI/a_n)\ast inv(PI/a_n,b_n)\ast b_n$

$X=x_1+x_2+...+x_n$

满足要求的最小$x_{min}=X\%PI$.

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

$Code:$

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 1E5+1;
typedef long long LL;
int n;
int a[MAXN],b[MAXN];
void exgcd(LL a,LL b,LL &x,LL &y)
{
	if(b == 0)
	{
		x = 1;
		y = 0;
		return;
	}
	exgcd(b,a%b,x,y);
	LL _x = x;
	x = y;
	y = _x-a/b*y;
}
LL inv(LL a,LL mod)
{
	LL x,y;
	exgcd(a,mod,x,y);
	return (x%mod+mod)%mod;
}
LL add(LL a,LL b,LL MOD)
{
	return (a%MOD+b%MOD)%MOD;
}
int main()
{
	scanf("%d",&n);
	LL PI = 1;
	for(int i=1;i<=n;i++)
	{
		scanf("%d%d",a+i,b+i);
		PI *= a[i];
	}
	LL sum = 0;
	for(int i=1;i<=n;i++)
	{
		sum = add(sum , b[i]*inv(PI/a[i],a[i])*(PI/a[i]) , PI);
	}
	printf("%lld",sum);
	return 0;
} 

3 扩展中国剩余定理

当方程组

$x\equiv b_1(mod{a}_1)$

$x\equiv b_2(mod{a}_2)$

$...$

$x\equiv b_n(mod{a}_n)$

中的$a_1,a_2,...,a_n$不两两互质时，需要采用扩展中国剩余定理算法。

取出其中的两个方程：

$x\equiv b_1(mod{a}_1)$

$x\equiv b_2(mod{a}_2)$

由同余的意义，有

$x=k_1{a}_1+b_1$

$x=k_2{a}_2+b_2$

消去$x$，得

$k_1{a}_1+b_1=k_2{a}_2+b_2$

$a_1{k}_1-k_2{a}_2=b2-b1$

为一个二元一次方程，可以用扩展欧几里得算法解出一个特解（$(1.3)$内容）$(k_1,k_2)$

由$(1.4)$结论，方程组有通解

$x=x_1+m\ast (b/gcd(a,b))$

$y=y_1-m\ast (a/gcd(a,b)),(m\in \mathbb{Z}).$

又

$x=(k_1+mr)a_1+b_1,$其中$r=a_2/gcd(a_1,a_2)$

所以

$x=k_1{a}_1+mr{a}_1+b_1$

$x=k_1{a}_1+mlcm(a_1,a_2)+b_1$

方程两边同模$lcm(a_1,a_2)$得

$x\equiv k_1{a}_1+b_1(modlcm(a_1,a_2))$

这样，两个方程就合并成了一个方程，求解最后一个方程的最小正整数解即可。

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

$Code:$

#include<cstdio>
#include<cctype>
typedef __int128 LL;
const int MAXN = 1E5+1;
int n;
LL ai[MAXN],bi[MAXN];
inline void read(__int128 &X)
{
	X = 0;
	int w=0; char ch=0;
	while(!isdigit(ch)) {w|=ch=='-';ch=getchar();}
	while(isdigit(ch)) X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
	if (w) X = -X;
}
void print(LL x)
{
	if (!x) return ;
	if (x < 0) putchar('-'),x = -x;
	print(x / 10);
	putchar(x % 10 + '0');
}
LL gcd(LL a,LL b)
{
	return !b?a:gcd(b,a%b);
}
LL exgcd(LL a,LL b,LL &x,LL &y)
{
	if(b == 0)
	{
		x = 1;
		y = 0;
		return a;
	}
	LL r = exgcd(b,a%b,x,y);
	LL x_ = x;
	x = y;
	y = x_-a/b*y;
	return r;
}
LL lcm(LL a,LL b)
{
	return a*b/gcd(a,b);
}
LL solve(LL a,LL b,LL c)
{
	LL x,y;
	exgcd(a,b,x,y);
	LL x1 = c/gcd(a,b)*x;
	LL mod = b/gcd(a,b);
	return (x1%mod+mod)%mod;
}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		read(ai[i]),read(bi[i]);
	for(int i=1;i<n;i++)
	{
		LL a1 = ai[i],a2 = ai[i+1];
		LL b1 = bi[i],b2 = bi[i+1];
		LL k1 = solve(a1,a2,b2-b1);
		ai[i+1] = lcm(a1,a2);
		bi[i+1] = k1*a1+b1;
	}
	print(bi[n]);
	return 0;
}

本博客为自己学习中的理解，如有错误还请指出。