【原创】 中国剩余定理 和 拓展中国剩余定理

孙子

Preface

数论学习Part 7。
每天进步一点点，退役不会太丢脸。

---

联赛只剩十五天，隔靴搔痒智何添？
剩余定理四天一篇，动归图论又几何，数据结构似云烟。
莫再等闲，莫忘时间。

CRT

我们小时候都做过这种问题，一个数$x$它满足$x\equiv \{\begin{array}{ll}{a}_1& \mathrm{m}\mathrm{o}\mathrm{d}{p}_1\\ a_2& \mathrm{m}\mathrm{o}\mathrm{d}{p}_2\\ & ⋮\\ a_n& \mathrm{m}\mathrm{o}\mathrm{d}{p}_n\end{array}$，然后让你求出一个最小的x出来。要求所有的$p$互质。

遇到这种问题以前的我一般会各写几项出来，然后找相同。
这个时候孙子站了出来，他说出：记$P=\prod p$，则$x={\sum }_{k=1}^{k=n}{a}_k\frac{P}{p_k}{t}_k$，$t_k是\frac{P}{p_k}关于p_k的逆元$。这是x的一个解，如果你逆元是最小正逆元，那这个x就是最小正解。

为什么呢？
证明啊……。

$先看x\equiv a_1\mathrm{m}\mathrm{o}\mathrm{d}{p}_1。$
$因为\frac{P}{p_1}与p_1互质，所以肯定找得到数t_1使得\frac{P}{p_1}{t}_1\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}{p}_1,这个t_1被称为\frac{P}{p_1}关于p_1的逆元。$
$那么就有x\equiv a_1\equiv a_1\times 1\equiv a_1\frac{P}{p_1}{p}_1\mathrm{m}\mathrm{o}\mathrm{d}{p}_1。$
$因为\frac{P}{p_1}=p_2\times p_3\times \cdots \times p_n，所以有a_1\frac{P}{p_1}{p}_1\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}{p}_2),(\mathrm{m}\mathrm{o}\mathrm{d}{p}_3),\cdots ,(\mathrm{m}\mathrm{o}\mathrm{d}{p}_n)。$←不标准的表达，应当扣分
$以此类推，a_2\frac{P}{p_2}{t}_2就满足模p_2时与x同余，模其它p的时候得到0。$
$那么这些东西的和——s={\sum }_{k=1}^{k=n}{a}_k\frac{P}{p_k}{t}_k，一定满足s\equiv a_1+0\times (n-1)\mathrm{m}\mathrm{o}\mathrm{d}{p}_1，$$s\equiv a_2+0\times (n-1)\mathrm{m}\mathrm{o}\mathrm{d}{p}_2，$$s\equiv a_3+0\times (n-1)\mathrm{m}\mathrm{o}\mathrm{d}{p}_3，\cdots ，$$s\equiv a_n+0\times (n-1)\mathrm{m}\mathrm{o}\mathrm{d}{p}_n$，所以可以认为$s$是原方程的一个解。那么通解就是$s+zP,z\in \mathbb{Z}$。

Q.E.D
就这？

代码

洛谷板子题，双倍经验。

曹冲养猪

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

typedef long long ll;
void Read(ll &p)
{
	p=0;
	int f=1;
	char c=getchar();
	while(c<'0' || c>'9') 
	{
		if(c=='-') f=-1;
		c=getchar();
	}
	while(c>='0' && c<='9')
		p=p*10+c-'0',c=getchar();
	p*=f;
}

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

ll n,ans,d,x,y,P=1,a[1024],p[1024];

int main()
{
	Read(n);
	for(ll i=1;i<=n;i++) Read(p[i]),Read(a[i]),P*=p[i];
	for(ll i=1;i<=n;i++) exgcd(P/p[i],p[i],d,x,y),ans+=a[i]*(P/p[i])%p*x,((ans%=P)+=P)%=P;
	

[TJOI2009]猜数字。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

typedef long long ll;
void Read(ll &p)
{
	p=0;
	int f=1;
	char c=getchar();
	while(c<'0' || c>'9') 
	{
		if(c=='-') f=-1;
		c=getchar();
	}
	while(c>='0' && c<='9')
		p=p*10+c-'0',c=getchar();
	p*=f;
}

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

ll mul(ll a,ll b,ll p)
{
	return ((a*b-(ll)((long double)a/p*b+1e-8)*p)%p+p)%p;
}

ll n,ans,d,x,y,P,a[1024],p[1024];

int main()
{
	Read(n),P=1;
	for(ll i=1;i<=n;i++) Read(a[i]);
	for(ll i=1;i<=n;i++) Read(p[i]),((a[i]%=p[i])+=p[i])%=p[i],P*=p[i];
	for(ll i=1;i<=n;i++) exgcd(P/p[i],p[i],d,x,y),ans+=mul(mul(a[i],P/p[i],P),x,P),((ans%=P)+=P)%=P;
	printf("%lld\n",ans);
}

EXCRT

现在这些p不一定互质了。

我们先看头两个方程吼，$\{\begin{array}{l}x\equiv a_1\mathrm{m}\mathrm{o}\mathrm{d}{p}_1\\ x\equiv a_2\mathrm{m}\mathrm{o}\mathrm{d}{p}_2\end{array}$。
那么就有$x=k_1\times p_1+a_1=k_2\times p_2+a_2\Rightarrow k_1\times p_1-k_2\times p_2=a_2-a_1$。
令$gcd(p_1,p_2)=g$，则同时除就有$k_1\frac{p_1}{g}-k_2\frac{p_2}{g}=\frac{a_2-a_1}{g}⟺k_1\frac{p_1}{g}\equiv \frac{a_2-a_1}{g}\mathrm{m}\mathrm{o}\mathrm{d}\frac{p_2}{g}$，其中$gcd(\frac{p_1}{g},\frac{p_2}{g})=1$。
注意，如果$g\nmid a_2-a_1$，则无解。
因为这两个数互质，所以可以求出逆元来，就有：$k_1\equiv \frac{a_2-a_1}{g}\times (\frac{p_1}{g}{)}^{-1}\mathrm{m}\mathrm{o}\mathrm{d}\frac{p_2}{g}$。
所以设$k_1=q\times \frac{p_2}{g}+\frac{a_2-a_1}{g}\times (\frac{p_1}{g}{)}^{-1}$。

知道$k_1$了，我们就知道了$x$。
$x=k_1\times p_1+a_1=\left(q\times \frac{p_2}{g}+\frac{a_2-a_1}{g}\times (\frac{p_1}{g}{)}^{-1}\right)\times p_1+a_1=q\times \frac{p_1{p}_2}{g}+\frac{a_2-a_1}{g}\times (\frac{p_1}{g}{)}^{-1}\times p_1+a_1$
$⟺x\equiv \frac{a_2-a_1}{g}\times (\frac{p_1}{g}{)}^{-1}\times p_1+a_1\mathrm{m}\mathrm{o}\mathrm{d}\frac{p_1{p}_2}{g}$。

我们把两个同余式合成了一个，太棒了。
看来我们可以做很多次，最后得到一个方程$x\equiv A\mathrm{m}\mathrm{o}\mathrm{d}B$。
$x$肯定可以随便解了。

代码

洛谷板子题，【模板】扩展中国剩余定理（EXCRT）。

居然一次过编译过样例，还直接A掉了。
我佛了。

太丑了。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

typedef long long ll;
void Read(ll &p)
{
	p=0;
	int f=1;
	char c=getchar();
	while(c<'0' || c>'9') 
	{
		if(c=='-') f=-1;
		c=getchar();
	}
	while(c>='0' && c<='9')
		p=p*10+c-'0',c=getchar();
	p*=f;
}

ll gcd(ll a,ll b){return !b?a:gcd(b,a%b);}
void exgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
	if(!b){d=a,x=1,y=0; return;}
	exgcd(b,a%b,d,y,x),y-=(a/b)*x;
}

ll mul(ll a,ll b,ll p)
{
	return ((a*b-(ll)((long double)a/p*b+1e-8)*p)%p+p)%p;
}

ll n,A,B,C,D,E,d,p,x,y;

int main()
{
	Read(n),n--,Read(B),Read(A);
	while(n--) Read(D),Read(C),d=gcd(B,D),exgcd(B/d,D/d,p,x,y),E=B/d*D,(A+=mul(mul((C-A)/d,x,E),B,E))%=E,B=E;
	printf("%lld\n",(A%B+B)%B);
}

CanKaoWenXian

CRT：其实是在我看了好几篇举个例子就说自己证好了的博客以后，再看到这篇博客的时候突然懂了，所以就以这篇博客为参考文献了。不过这篇文章确实给了我参考。

---

EXCRT
贞德真的好看。（不是贞德当我没说）

说在后面

问了问数竞的同学，他们的中国剩余定理是证明其有解，然后数学归纳法就整出来了。
好吧，他们确实算不出来逆元。