关闭

(ACM数论)中国剩余定理(孙子定理)

标签: acm
244人阅读 评论(0) 收藏 举报
分类:

中国剩余定理是一种能求解一次同余方程组的姿势~

首先我们来看看什么是一次同余方程组:

X % a[1] = b[1]

X % a[2] = b[2]

….

X % a[n] = b[n]

这里的 % 是取余数的意思,使用这个姿势的前提是a[1]、a[2]、a[3]、… 、a[n]两两互质。

令 X = X[1] + X[2] + … + X[n],此时只要找到一组 X[i] 满足 X[i] % a[i] = b[i] 且 X[i] % a[1…n(除了i)]= 0 成立即可。

要使 X[i] % a[1…n(除了i)]= 0 成立,只需要让 X[i] 等于他们乘积的倍数即可。令M为他们乘积,于是有 X[i] = k * M(其中k为任意整数)。

将 X[i] = k * M 代入 X[i] % a[i] = b[i] 得到 (k * M) % a[i] = b[i]

即 k * M + o * a[i] = b[i] (其中o为任意整数)

两边同时除以 b[i] 可得到 k’ * M + o’ * a[i] = 1

即 k’ * M % a[i] = 1

此时 k’ 为 M 对 a[i](质数)的 乘法逆元 (关于乘法逆元可以参考我的这篇文章:http://blog.csdn.net/kasumimasami/article/details/52354715

所以对式子 k’ * M + o’ * a[i] = 1 两边同时乘以 b[i] 可得到

k’ * b[i] * M + o’ * b[i] * a[i] = b[i] 即 k’ * b[i] * M % a[i] = b[i]

所以 X[i] = k’ * b[i] * M

所以 X = ni[1] * b[1] * M[1] + ni[2] * b[2] * M[2] + … + ni[n] * b[n] * M[n],其中 ni[i] 为 M[i] 对 a[i](质数)的 乘法逆元

代码如下:

int fastpow(int a,int b,int mod)// a ^ b % mod
{
    int r = 1;
    while(b)
    {
        if(b&1)
            r = r * a % mod;
        a = a * a % mod;
        b = b >> 1;
    }
    return r;
}

int CRT(int a[], int b[], int len)//数组a,b如上定义,len为方程的个数
{
    int m = 1 , X = 0;
    for(int i = 1; i <= len; i++)
    {
        m = m * a[i];
    }
    for(int i = 1; i <= len; i++)
    {
        int M = m / a[i];//除了a[i]的乘积M[i]
        int ni = fastpow(M,(a[i]-2),a[i]);//求M[i]关于a[i]的乘法逆元ni[i]
        X = X + ni * M * b[i];
    }
    //返回满足的最小正整数
    while(X < 0)
        X = X + m;
    return X % m;
}
1
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:3716次
    • 积分:121
    • 等级:
    • 排名:千里之外
    • 原创:8篇
    • 转载:0篇
    • 译文:0篇
    • 评论:2条
    文章分类
    文章存档
    最新评论