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

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

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

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;
}
  • 3
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值