扩展欧几里德 中国剩余定理 合并模线性方程组
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
我们来看一个应用 http://poj.org/problem?id=3708
f(j) = a[j] (0
f(d*n+j) = d*f(n) + b[j] (0<=j
对于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进制的本质，通过置换群的概念 抽化出线性方程组，最终通过合并来解决