算法——中国剩余定理
中国剩余定理
中国剩余定理又称中国余数定理,是数论中一个关于一元线性同余方程组的定理,说明了一元线性同余方程组有解的准则以及求解方法。
求解方法
中国剩余定理原理
主要用到两个定理:
a mod b = c 等价于 (a + nb) mod b = c;
a mod b = c 等价于 2a mod 2b = 2c;
中国剩余定理代码表示
void Gcd(int a, int b, int &x, int &y){
//根据欧几里德定理
if(b == 0){
//任意数与0的最大公约数为其本身。
x = 1;
y = 0;
}else{
int x1, y1;
Gcd(b, a % b, x1, y1);
if(a*b < 0){
//异号取反
x = - y1;
y = a/b*y1 - x1;
}else{
//同号
x = y1;
y = x1 - a/b* y1;
}
}
}
// a为被除数组,m为余数,k为被除数组个数1
int China_res_con(int a[],int m[],int k)
{
int N[k];
int min_com_mult = 1; //最小公倍数
int result = 0;
for(int i = 0;i < k;i++)
{
min_com_mult *= m[i];
}
for(int j = 0;j < k; j++)
{
int L, J;
Gcd(min_com_mult/m[j],-m[j],L,J);
N[j] = m[j] * J + 1;//1
N[j] = mm/m[j] * L;//2 【注】1和2这两个值应该是相等的。
result += N[j]*a[j];
}
return result % min_com_mult;
}
参考文献
[1] 算法导论(第三版)
[2] hard_man.中国剩余定理(孙子定理)的证明和c++求解.2012.07.10