数论常用内容——中国剩余定理

xiaoxiao2021-02-28  74

中国剩余定理

中国剩余定理是中国古代求解一次同余式组(见同余)的方法。是数论中一个重要定理。又称孙子定理。

用现代数学的语言来说明的话,中国剩余定理给出了下图所示的一元线性同余方程组有解的判定条件,并用构造法给出了在有解情况下解的具体形式:

说明:假设整数m1,m2, … ,mn两两互质,则对任意的整数:a1,a2, … ,an,方程组 有解,并且通解可以用如下方式构造得到:

设 M = m1*m2*…*mn 是整数m1,m2, … ,mn的乘积,并设 Mi = M/mi 是除了mi以外的n-1个整数的乘积。

设 ti = Mi^(-1) 为 Mi 模 mi 的数论倒数(ti 为 Mi 模 mi 意义下的逆元) Mi*ti ≡ 1(mod mi),i∈{1,2,3,…,n} 方程组 (S) 的通解形式为 x = k*M + ∑(ai*ti*Mi),i∈{1,2,3,…,n},k∈Z 在模 M 的意义下,方程组(S)只有一个解如下图所示

代码实现

互质

//求M%A=a,M%B=b,...中的M,其中A,B,C...互质 int CRT(int a[],int m[],int n){ int M = 1; int ans = 0; for(int i=1; i<=n; i++) M *= m[i]; for(int i=1; i<=n; i++){ int x, y; int Mi = M / m[i]; ex_gcd(Mi, m[i], x, y); ans = (ans + Mi * x * a[i]) % M; } if(ans < 0) ans += M; return ans; }

非互质

一般的中国剩余定理要求mi两两互质,但是保证互质条件太苛刻了,若mi并不满足两两互质时,就要采用两两合并的思想,假设要合并如下两个方程 x=a1+m1*x1 x=a2+m2*x2

得到 a1+m1*x1 = a2+m2*x2 → m1*x1+m2*x2 = a2-a1

再通过扩展欧几里得算法解出x1的最小正整数解,代入 x=a1+m1*x1

得到x后合并为一个方程的结果为 y ≡ x(mod lcm(m1,m2))

这样一直合并下去,最终可以求得同余方程组的解。

下面来看代码:

bool merge(LL a1, LL m1, LL a2, LL m2, LL &a3, LL &m3) { LL d = gcd(m1, m2); LL c = a2 - a1; if(c % d) return false; c = (c % m2 + m2) % m2; m1 /= d; m2 /= d; c /= d; c *= Inv(m1, m2);//Inv为乘法逆元,数论常用内容——欧几里得算法与扩展欧几里得算法 c %= m2; c *= m1 * d; c += a1; m3 = m1 * m2 * d; a3 = (c % m3 + m3) % m3; return true; } LL CRT(LL a[], LL m[], int n) { LL a1 = a[1]; LL m1 = m[1]; for(int i=2; i<=n; i++) { LL a2 = a[i]; LL m2 = m[i]; LL m3, a3; if(!merge(a1, m1, a2, m2, a3, m3)) return -1; a1 = a3; m1 = m3; } return (a1 % m1 + m1) % m1; }
转载请注明原文地址: https://www.6miu.com/read-71132.html

最新回复(0)