【ACM】中国剩余定理(CRT)

中国剩余定理是求解一次线性同余方程组的方法。

中国剩余定理:
假设整数m1,m2, … ,mn两两互素,则同余方程组
这里写图片描述
有整数解。
这里写图片描述 是m1,m2,m3…mn的乘积,并设 这里写图片描述是除了mi以外的n-1个整数的乘积。
这里写图片描述是Mi模mi的逆元(数论中的倒数)。

那么在模M下的解是唯一的。

这里写图片描述

模板代码如下:

void Extend_Euclid(int a,int b,int& x,int& y)
{
    if(b == 0)
    {
        x = 1;
        y = 0;
        return; //表示程序退出,但是没有返回值
    }
    Extend_Euclid(b,a%b,x,y);
    int temp = x;
    x = y;
    y = temp - (a/b) * y;
}

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];
        Extend_Euclid(Mi,m[i],x,y); // 求逆元
        ans = (ans + a[i]*Mi*x) % M;
    }
    if(ans<0)
        ans += M;
    return ans;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值