中国剩余定理模板(互质的情况)

在《孙子算经》中有这样一个问题:“今有物不知其数,三三数之剩二(除以3余2),五五数之剩三(除以5余3),七七数之剩二(除以7余2),问物几何?”这个问题称为“孙子问题”,该问题的一般解法国际上称为“中国剩余定理”。具体解法分三步:

    1. 找出三个数:从3和5的公倍数中找出被7除余1的最小数15,从3和7的公倍数中找出被5除余1 的最小数21,最后从5和7的公倍数中找出除3余1的最小数70。
    2. 用15乘以2(2为最终结果除以7的余数),用21乘以3(3为最终结果除以5的余数),同理,用70乘以2(2为最终结果除以3的余数),然后把三个乘积相加(15*2+21*3+70*2)得到和233。
    3. 用233除以3,5,7三个数的最小公倍数105,得到余数23,即233%105=23。这个余数23就是符合条件的最小数。


在理解这个计算过程是非常简单的,但真正在代码中实现却不然,刚开始学知道模板对就套,现在有时间来详解一下中国剩余定理的模板
(感觉网上的大神们没有解释代码的(对他们来说太简单了!^^)),接下来我就详解一下中国剩余定理的模板,拿15*2来举例吧,代码中是这样呈现的
d=Extended_Euclid(w[i],m,x,y);
        ret=(ret+y*m*a[i])%n;
15是3和5的最小公倍数,现在我们需要找3和5的公倍数中除7余2的那个,即(15*x)%7==2,我们知道(15*x)%7,=gcd(15,7)中x为15的逆元
所以我们所求的x其实是15 的逆元,有因为所有数互质所以gcd(15,7)=1,所以15*x+7*y=gcd(15,7)所求除的值是15*x+7*y=1的解,
所以最后还需要乘以w[i](这里就是2了)

模板:
#include <iostream>
#include <cstdio>
using namespace std;


int Extended_Euclid(int a,int b,int &x,int &y)    //扩展欧几里得算法
{
    int d;
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    d=Extended_Euclid(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
int Chinese_Remainder(int a[],int w[],int len)    //中国剩余定理  a[]存放余数  w[]存放两两互质的数
{
    int i,d,x,y,m,n,ret;
    ret=0;
    n=1;
    for (i=0; i<len; i++)
        n*=w[i];
    for (i=0; i<len; i++)
    {
        m=n/w[i];
        d=Extended_Euclid(w[i],m,x,y);
        ret=(ret+y*m*a[i])%n;
    }
    return (n+ret%n)%n;
}
int main()
{
    int n,i;
    int w[15],b[15];
    while (scanf("%d",&n),n)
    {
        for (i=0; i<n; i++)
        {
            scanf("%d%d",&w[i],&b[i]);
        }
        printf("%d\n",Chinese_Remainder(b,w,n));
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值