数论——同余

定义：

若 $a,b$ 为两个整数，满足 $a\mathrm{mod}m=b\mathrm{mod}m$ 或 $(a-b)|m$ 则称 $a$ 和 $b$ 关于模 $m$ 同余，记为 $a\equiv b(\mathrm{mod}m)$ 。

性质：

1. $a\equiv a(\mathrm{mod}m)$

2. 若 $a\equiv b(\mathrm{mod}m)$ ， 则 $b\equiv a(\mathrm{mod}m)$ 。

3. 若 $a\equiv b(\mathrm{mod}m)\text{ , }b\equiv a(\mathrm{mod}m)$ 。

4. 若 $a\equiv b(\mathrm{mod}m)$ ， 则 $a+c\equiv b+c(\mathrm{mod}m)$ 。

5. 若 $a\equiv b(\mathrm{mod}m)$ ， 则 $a\times c\equiv b\times c(\mathrm{mod}m)$ 。

6. $a\times b\mathrm{mod}k=(a\mathrm{mod}k\times b\mathrm{mod}k)\mathrm{mod}k$ 。

7. $(a+b)\mathrm{mod}k=(a\mathrm{mod}k+b\mathrm{mod}k)\mathrm{mod}k$ 。

8. 若 $a\equiv b(\mathrm{mod}m)$ ， 则 $a^n\equiv b^n(\mathrm{mod}m)$ 。

9. 若 $a\equiv b(\mathrm{mod}m)\text{ , }c\equiv d(\mathrm{mod}m)$ ， 则 $a\times c\equiv b\times d(\mathrm{mod}m)$ 。

证明：
有：

$$a\times c\mathrm{mod}m=(a\mathrm{mod}m\times c\mathrm{mod}m)\mathrm{mod}m$$

$$b\times d\mathrm{mod}m=(b\mathrm{mod}m\times d\mathrm{mod}m)\mathrm{mod}m$$

$$a\mathrm{mod}m=b\mathrm{mod}m\text{ , }c\mathrm{mod}m=d\mathrm{mod}m$$

则：

$$a\times c\mathrm{mod}m=b\times d\mathrm{mod}m$$

所以：

$$a\times c\equiv b\times d(\mathrm{mod}m)$$

10. 若 $a\times p\equiv x(\mathrm{mod}m)\text{ , }a\times q\equiv x(\mathrm{mod}m)$ ， $p,q$ 互质 ， 则 $a\mathrm{mod}p\times q=x$ 。

证明：

有 $s,t$ 使

$$a=s\times p+x\text{ , }a=t\times q+x$$

则：

$$s\times p=t\times q$$

因为 $p,q$ 互质，所以 $p|t\text{ , }q|s$ 。

那么一定有 $r$ 使

$$s=q\times r$$

则：

$$a=r\times p\times q+x$$

所以：

$$a\mathrm{mod}p\times q=x$$

前置芝士gcd

拓展欧几里得算法

用途： 求 $a\times x+b\times y=\gcd (a,b)$

解法：

有：

$$\gcd (a,b)=\gcd (b,a\mathrm{mod}b)$$

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

所以：

$$b\times x^{\prime }+(a\mathrm{mod}b)\times y^{\prime }=a\times x+b\times y$$

因为 $a\mathrm{mod}b=a-\lfloor a/b\rfloor \times b$

所以右式整理得：

$$b\times x^{\prime }+(a\mathrm{mod}b)\times y^{\prime }$$

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

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

所以

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

就可以递归求解：

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

定理1：

若同余方程 $a\times x+b\times y=c$ 有整数解，则有 $\gcd (a,b)|c$

定理2：

若 $\gcd (a,b)=1$ 且 $x_0,y_0$ 为 $a\times x+b\times y=c$ 的一组解，那么方程的任意一组解为 $x=x_0+b\times t,y=y_0+a\times t$

特解：

求 $x$ 的最小正整数解，即 $t=b/\gcd (a,b)\text{ , }x=(x\mathrm{mod}t+t)\mathrm{mod}t$

例题：青蛙的约会

题目传送门

思路

因为两个青蛙同时跳到同一个点上才算碰面，设 $t$ 为跳的次数， $p$ 为两个青蛙跳的圈数之差，有如下式子：

$$(x+m\times t)-(y+n\times t)=p\times L$$

整理得：

$$(n-m)\times t+L\times p=x-y$$

首先，要判断 $\gcd (n-m,L)\nmid x-y$ 的情况。

然后，直接套拓展欧几里得算法，求出 $(n-m)\times t+L\times p=\gcd (n-m,L)$ 一组解 $t_0\text{ , }{p}_0\text{ , }gcd=\gcd (n-m,L)$ 。

带入特解（本题求的是 $t$ 的最小解），求出 $(n-m)\times t+L\times p=\gcd (n-m,L)$ 中 $t_0$ 的最小解 $t^{\prime }$

$$t^{\prime }=(t_0\mathrm{mod}\frac{L}{gcd}+\frac{L}{gcd})\mathrm{mod}\frac{L}{gcd}$$

有

$$t=t^{\prime }\times \frac{x-y}{gcd}$$

所以

$$t=(\frac{x-y}{gcd}\times t_0\mathrm{mod}\frac{L}{gcd}+\frac{L}{gcd})\mathrm{mod}\frac{L}{gcd}$$

$Code$

#define int long long
signed main()
{
	int x,y,m,n,l;
	int X,Y,gcd;
	cin>>x>>y>>m>>n>>l;
	int A=x-y,B=n-m;
	if(B<0)//注意取反！！！
	{
		B=-B;
		A=-A;
	}
	if(A%(gcd=exgcd(B,l,X,Y))!=0)
		cout<<"Impossible"<<endl;
	else
		cout<<(  A/gcd * X % (l/gcd) + (l/gcd)  ) % (l/gcd)<<endl;
}

中国剩余定理（CRT）

思想

用于求：

$$\{\begin{array}{l}x\equiv a_1(\mathrm{mod}mo{d}_1)\\ x\equiv a_2(\mathrm{mod}mo{d}_2)\\ x\equiv a_3(\mathrm{mod}mo{d}_3)\\ ...\\ x\equiv a_{n-2}(\mathrm{mod}mo{d}_{n-2})\\ x\equiv a_{n-1}(\mathrm{mod}mo{d}_{n-1})\\ x\equiv a_n(\mathrm{mod}mo{d}_n)\end{array}$$

其中 $mo{d}_1,mo{d}_2,mo{d}_3,...,mo{d}_n$ 互质。

接下来以求如下同余方程组为例：

$$\{\begin{array}{l}x\equiv 2(\mathrm{mod}3)\\ x\equiv 3(\mathrm{mod}5)\\ x\equiv 2(\mathrm{mod}7)\end{array}$$

先求如下三个方程组：

$$(1)\{\begin{array}{l}x\equiv 1(\mathrm{mod}3)\\ x\equiv 0(\mathrm{mod}5)\\ x\equiv 0(\mathrm{mod}7)\end{array}(2)\{\begin{array}{l}x\equiv 0(\mathrm{mod}3)\\ x\equiv 1(\mathrm{mod}5)\\ x\equiv 0(\mathrm{mod}7)\end{array}(3)\{\begin{array}{l}x\equiv 0(\mathrm{mod}3)\\ x\equiv 0(\mathrm{mod}5)\\ x\equiv 1(\mathrm{mod}7)\end{array}$$

以 $(1)$ 为例：

求

$$\{\begin{array}{l}x\equiv 1(\mathrm{mod}3)\\ x\equiv 0(\mathrm{mod}5)\\ x\equiv 0(\mathrm{mod}7)\end{array}$$

因为 $x\equiv 0(\mathrm{mod}5)x\equiv 0(\mathrm{mod}7)$ ，所以 $x$ 是 $35$ 的倍数，等同于求

$$35y\equiv 1(\mathrm{mod}3)(x=35y)$$

等同于求

$$2y\equiv 1(\mathrm{mod}3)(x=35y)$$

所以

$$x\equiv (y=2)\times 35\equiv 70(\mathrm{mod}(105=3\times 5\times 7))$$

所以解为：

$$(1)=70(2)=21(3)=15$$

同理，有：

$$(1)\{\begin{array}{l}x\equiv 2(\mathrm{mod}3)\\ x\equiv 0(\mathrm{mod}5)\\ x\equiv 0(\mathrm{mod}7)\end{array}=35(2)\{\begin{array}{l}x\equiv 0(\mathrm{mod}3)\\ x\equiv 3(\mathrm{mod}5)\\ x\equiv 0(\mathrm{mod}7)\end{array}=63(3)\{\begin{array}{l}x\equiv 0(\mathrm{mod}3)\\ x\equiv 0(\mathrm{mod}5)\\ x\equiv 2(\mathrm{mod}7)\end{array}=30$$

有：

$$\{\begin{array}{l}x\equiv 2(\mathrm{mod}3)\\ x\equiv 3(\mathrm{mod}5)\\ x\equiv 2(\mathrm{mod}7)\end{array}=((1)+(2)+(3))\mathrm{mod}105=128\mathrm{mod}105=23$$

其中 $\mathrm{mod}105$ 的原因是 $x$ 可以同时去掉 $3,5,7$ 使 $x$ 更小，并且不会影响正确性。

$Code$

int China()
{
	int x,y,ans=0,M,N=1;
	for(int i=1; i<=n; i++)
	N*=mod[i];
	for(int i=1; i<=n; i++)
	{
		M=N/mod[i];
		exgcd(mod[i],M,x,y);//直接套模板。
		ans=(ans+(y*M*a[i])%N+N)%N;
	}
	if(a>0)return ans;
	else
	return ans+N;
}

逆元

思想

用于求 $a\times x\equiv 1(\mathrm{mod}b)$ （ $a,b$ 互质） 的 $x$ ， 称 $x$ 为 $a$ 的逆元，记为 $a^{-1}$ 。

可用拓展欧几里得算法计算。

也可用费马定理计算。

费马定理： 若 $a,p$ 互质且 $p$ 为素数， $a^{p-1}\equiv 1(\mathrm{mod}p)$

还可用递推计算。

首先：

$$1^{-1}\equiv 1(\mathrm{mod}p)$$

设 $p=k\times i+r,r<i,1<i<p$

有

$$k\times i+r\equiv 0(\mathrm{mod}p)$$

两边同乘上 $i^{-1}$ , $r^{-1}$ 。

$$k\times r^{-1}+i^{-1}\equiv 0(\mathrm{mod}p)$$

$$i^{-1}\equiv -k\times r^{-1}(\mathrm{mod}p)$$

$$i^{-1}\equiv -\frac{p}{i}\times (p\mathrm{mod}i{)}^{-1}(\mathrm{mod}p)$$

$Code$

long long exgcd_inv(long long a,long long mod)
{
	long long d,x,y;
	d=exgcd(a,n,x,y);
	return d==1?(x+n)%n:-1;
}

long long fermat_inv(long long a,long long mod)
{
	return pow_mod(a,mod-2,mod);
}

inv[i]=-(p/i)*inv[p%i];

威尔逊定理

若 $p$ 为素数，则 $(p-1)!\equiv -1(\mathrm{mod}p)$ 。 若 $(p-1)!\equiv -1(\mathrm{mod}p)$ 为素数，则 $p$ #

费马定理

若 $a,p$ 互质且 $p$ 为数， $a^{p-1}\equiv 1(\mathrm{mod}p)$

费马小定理

若 $p$ 为素数， $a^p\equiv a(\mathrm{mod}p)$

未完待续…