扩展欧几里德 中国剩余定理 合并模线性方程组

1.1.1 扩展欧几里得

要说扩展必须先从它的非扩展版本说起，对于求两个数的最大公约数，我们有辗转相除法，其核心就是gcd(a,b)=gcd(b,a%b) (a>=b) (1)为什么呢，我们来证明一下令a=k*b+t 则a%b=t ,若设d是a，b的一个公约数，a%d==0 k*b%d==0 又因为(k*b+t)%d==0 所以t%d==0，这个d包含了a和b的最大公约数，于(1)得证。有了这个作为基础我们来看下扩展欧几里得到底是个什么东东，对于这样一个不定方程ax+by=c 方程有解当且仅当 c%gcd(a,b)==0,为什么呢，因为(ax+by)%gcd(a,b)==0,得证。于是对于ax1+by1=gcd(a,b) 该方程的一组解(x1,y1) 只需都乘上c/gcd(a,b)，就可以得到原方程的一组(x,y)解。我们现在只需要考虑ax1+by1=gcd(a,b)上.由于gcd(a,b)=gcd(b,a%b),故a*x1+b*y1=gcd(a,b)=gcd(b,a%b)=b*x2+(a%b)*y2,令k=a/b r=a%b a=k*b+r,r=a-k*b将其带入不定方程得a*x1+b*y1=b*x2+(a-(a/b)*b)y2=a*y2+b*(x2-(a/b)*y2) (x1=y2,y1=x2-(a/b)*y2),我们还可得到通解

xx=x1+b/gcd(a,b)*t

yy=y1- a/gcd(a,b)*t (t为任意正整数，x1,y1为一组基解)，扩展欧几里得一个很大的应用就是求乘法逆元，逆元是什么？离散数学全还给老师了-_-||,定义群G，x∈G，若存在x'∈G，e是单位元，x'· x=e则称x'是x的左逆元 x·x'=e,则称x'是x的右逆元,若两者一样，则x'是x的逆元，注意在G中元素可能没有逆元。在模运算中a×b ≡1 (mod m)，则b是a模m的乘法逆元。乘法逆元1/b(mod m)存在当且仅当gcd(b,m)=1，且满足1/b (mod m)存在的b的个数为phi(m)个。下面来看下如何求逆元，对于1/a (mod n) = x 若存在，则可以写成ax %n =1 即 ax - (ax/n)*n=1 等价于 ax-ny=1,求解即可。

1.1.2 中国剩余定理

对于所有ai%mi=x方程组 (1<=i<=n),m值两两互素，且都>1,则存在解x满足上面方程组，对于另一个x' != x 满足 x' % M = x 其中M=m1*m2....*mn。这个证明参看《计算数论》p 112 - p 113 。其解得形式可以写成 x sigma( ai Mi Mi' % m) (1<=i<=n) 其中m =m1*m2...mn, Mi=m/mi ,Mi* Mi' %mi =1。

1.1.3 合并模线性方程组

先简单介绍下线性同余方程 ax%n=b 其有解需满足(a,n)|b,且其解得个数有(a,n)个，令d=gcd(a,n) 设t为(a/d)x % (n/d) = (b/d)的唯一解，则ax%n=b的d个解为t,t+n/d,t+2n/d....t+(d-1)n/d。证明其有解过程如下：若有解，则(a,n)|b满足，设s为一个解，则as%n=b，则n|(as-b),故as-b=kn(k为整数) 因为 (a,n)|a,(a,n)|kn,所以(a,n)|b证毕。

问题模型：给定模线性同余方程组 r[i]%a[i]=m (1<i<=n,r[i],a[i]均为正整数)，判断是否存在最小非负整数解，为一般方程不考虑中国剩余定理的情形。先看前两个方程a1*x+r1=m a2*y+r2=m, 可以写成a*x-b*y=r2-r1 方程有解需保证(a,b)|(r2-r1),问题等价于寻求ax+by=n的解，我们利用扩展欧几里得先求得ax+by=(a,b)的一组解(x0,y0)，于是(x1,y1)=(x0,y0)*n/(a,b)为ax+by=n的解。其他解 x_1,2...≡x1+(b/(a,b))*t,y_1,2...≡y1-(a/(a,b))*t (0 <=t<=(a,b)-1)。设d=gcd(a,b) x=x+(a2/d)*t,则x的最小非负整数解为 x=(x%(a2/d)+(a2/d))%(a2/d); 令f=a2/d;则x=(x%f+f)%f;新的方程变为a1*(x+(a2/d)*t)+r1=m 我们化简得(a1*a2/d)*t+(a1*x+r1)=m r1=a1*x+r1; a1=lcm(a1,a2)=(a1*a2)/d,故m=r1(mod a1)两个线性同余方程合并成一个等效方程，并继续合并下去,最后r1就是方程组的解了，该方法复杂度为O(n*extend_gcd)。

我们来看一个应用 http://poj.org/problem?id=3708
f(j) = a[j] (0<j<d)

f(d*n+j) = d*f(n) + b[j] (0<=j<d) 有点d进制的感觉，这是一个很好的思维突破口

对于f(N) 把N表示成d进制为N=x[m]*d^m+x[m-1]*d^(m-1)+....+x[1]*d+x[0]

则f(N)= d*(f(x[m]*d^(m-1)))+...d*f(x[1])+b[x[0]]

......

f(N)=d^m*a[x[m]]+d^(m-1)*b[x[m-1]]+....d*b[x[1]]+b[x[0]]

目标是经过若干次f变换，得到M

(x[0],x[1],...x[m])是一个群，而b数组的作用就是对任意一个x[i]起到一个置换的作用，把他变成了b[x[i]]，从而让其进入到一个循环中

我们可以获得每个x[i]的循环节大小w1[i]以及到达目标的步长w2[i] （关于群，置换，循环节的知识将会在以后的文章中出现）

设经历了p次变换，于是对于任何一个i必有 w1[i]*ki+w2[i]=p (注意a部分只有1个，我们可以单独计算其循环节及步长，最终获得一个类似的方程)

接下来的事情就可以利用上面所说得合并方程组来搞定了，这个问题的关键在于识破d进制的本质，通过置换群的概念 抽化出线性方程组，最终通过合并来解决