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

中国剩余定理

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

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

说明:假设整数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;  
}  
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值