同余与二元一次方程求解(上)

同余方程

同余：

设m是正整数，若$a$和$b$是整数，且$m|(a-b)$（意思是$a-b$能被$m$整除），则$a$和$b$模$m$同余$，a\equiv b(\mathrm{mod}m)$ ，也就是说$a$除以$m$的余数等于$b$除以$m$的余数

现在我们考虑一元线性同余方程$ax\equiv b(\mathrm{mod}m)$ $,a>0,n>0$
它表示$ax-b$是$m$的倍数，设为$-y$倍，则有二元线性丢番图方程$ax+my=b$，即求解一元线性同余方程等价于求解二元线性丢番图方程

如何求二元线性丢番图方程呢？

先给出一个推论：当且仅当$gcd(a,m)|b$的时候，方程对于未知量$x$有$gcd(a,m)$个模$m$不同余的解，另外，如果$(x_0,y_0)$是方程的一个特解,那么所有的通解可以表示为$x=x_0+n(b/gcd(a,m)),y=y_0-n(a/gcd(a,m))$，其中$n$为任意正整数

证明：
令$a=gcd(a,m)a^{\prime },m=gcd(a,m)m^{\prime }$,有$ax+my=gcd(a,m)(a^{\prime }x+m^{\prime }y)=b$,所以如果$a$和$m$都是整数，那$b$也必须是$gcd(a,m)$的整数倍才有整数解

所以问题转化为怎么求出这个二元方程$ax+my=b$的一个特解$(x_0,y_0)$,既然$gcd(a,m)|b$，不如先求出$a{x}^{\prime }+m{y}^{\prime }=gcd(a,m)$，最后左右两边都乘以$b/gcd(a,m)$,

$ax+my=gcd(a,m)$ 扩展欧几里得算法

• 当$m=0$时，$x=1,y=0$时等式成立
• 当$m>0$时，已知：$gcd(a,m)=gcd(m,a-m)=gcd(m,a$%$m),a$%$m=a-m\ast a/m$
  则有$mx+(a$%$m)y=gcd(m,a$%$m)=gcd(a,m)$
  $mx+(a-m\ast a/m)y=gcd(m,a$%$m)=gcd(a,m)$
  $ay+mx-(a/m)\ast my=gcd(a,m)$
  $ay+m(x-a/m\ast y)=gcd(a,m)$

对比$mx+a$%$m\ast y=gcd(a,m)$
所以$x=y,$$y=(x-a/m\ast y),$$a=m,$$m=a$%$m$

typedef long long ll;
ll extend_gcd(ll a,ll m,ll &x,ll &y){//返回特解x0,y0
	if(m==0){
		x=1;y=0;
		return a;
	}else{
		ll d=extend_gcd(m,a%m,y,x);//swap(x,y)
		y-=(a/m)*x;//x，y交换了，这句容易出错
		return d;//返回gcd(a,m)
	}
}

得到了$a{x}^{\prime }+m{y}^{\prime }=gcd(a,m)$的特解$(x_0^{\prime },y_0^{\prime })$
则$ax+my=b$的特解为$x_0=x_0^{\prime }\ast b/gcd(a,m),y_0=y_0^{\prime }\ast b/gcd(a,m)$
通解为$x=x_0+n\ast (b/gcd(a,m)),y=y_0-n\ast (a/gcd(a,m))$

x=x0 + n*(b/__gcd(a,m));
y=y0 - n*(b/__gcd(a,m));

对于同余方程$ax\equiv b(\mathrm{mod}m)$ $,a>0,n>0$
调用$extend$_$gcd(a,m,x,y)$，得出的$x_0$就是一个解，通解$x=x_0+(m/gcd(a,m))\ast n$，$0\le n\le d-1$.,当$n从0到d-1$,就取遍了模d的完全剩余系

x=x0+(m/__gcd(a,m))*n;

P1082

题目描述

求关于 $x$ 的同余方程 $ax\equiv 1(\mathrm{mod}b)$ 的最小正整数解。

输入格式

一行，包含两个整数 $a,b$，用一个空格隔开。

输出格式

一个整数 $x_0$，即最小正整数解。输入数据保证一定有解。

样例 #1

样例输入 #1

3 10

样例输出 #1

7

数据规模与约定

• 对于 $100\%$ 的数据，$2\le a,b\le 2,000,000,000$。

$ax\equiv 1(\mathrm{mod}b)$即方程$ax+my=1$，先用$extend$_$gcd()$求出一个特解，然后用取模运算来得出最小正整数解

x=(x%m+m)%m;//保证返回正整数解

---

补充

$ax\equiv 1(\mathrm{mod}b)$的解，为$a$模$m$的逆，记为$a^{-1}$,所以扩展欧几里得算法就可以用来求逆，下面再介绍一种求逆的方法
费马小定理：设n是素数，a是正整数且与n互素，有 $a^{n-1}\equiv 1(\mathrm{mod}n)$

$a\ast a^{n-2}\equiv 1(\mathrm{mod}n)$，那按照$a^{n-2}(\mathrm{mod}n)$就是$a$模$n$的逆！！！这个部分计算用快速幂取模来操作
但是！！！仅用费马小定理来测试素数是很危险的，原因是存在少量的伪素数也满足费马小定理，被称为Carmichael数,一亿个数里面大约有200多个，
欧几里得算法（辗转相除法）

int gcd(int a,int b){
	return b==0?a:gcd(b,a%b);
}

Thanks♪(･ω･)ﾉ