中国剩余定理

模板

#include<stdio.h>
#include <iostream>
using namespace std;
//扩展欧几里得算法
int exgcd(int a,int b,int &x,int &y)
{
    int d;
    if(b==0)
    {
        x=1;y=0;
        return a;
    }
    d=exgcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
 //中国剩余定理 ,r[]存放余数 ,prime[]存放两两互质的数
int Chinese_Remainder(int r[],int prime[],int len)
{
    int i,d,x,y,m,n=1,sum=0;
    //计算所以除数的积n，也是所以除数的最小公倍数
    for(i=0;i<len;i++)
        n*=prime[i];
    //计算符合所以条件的数
    for(i=0;i<len;i++)
    {
        m=n/prime[i];//计算除去本身的所有除数的积m
        d=exgcd(prime[i],m,x,y);//计算w[i]*x+m*y=gcd(w[i],m)的一个解y
        //累加整数解y的同并不断对n取余，其利用公式:(a+b)%c=(a%c+b%c)%c
        sum=(sum+y*m*r[i])%n;
    }
    return (n+sum%n)%n;//满足所以方程的最小解
}
int main()
{
    int n,i;
    int prime[15],r[15];
    while (printf("请输入组数n：\n"),scanf("%d",&n)!=EOF)
    {
        printf("请依次输入每组的除数和余数：\n");
        for (i=0;i<n;i++)
        {
            scanf("%d%d",&prime[i],&r[i]);
        }
        //printf("%d\n",Chinese_Remainder(b,w,n));
        printf("符合条件的最小整数：%d\n\n",Chinese_Remainder(r,prime,n));
    }
    return 0;
}

核心代码：

int Chinese_Remainder(int r[],int prime[],int len)
{
    int i,d,x,y,m,n=1,sum=0;
    for(i=0;i<len;i++)
        n*=prime[i];
    for(i=0;i<len;i++)
    {
        m=n/prime[i];
        d=exgcd(prime[i],m,x,y);
        sum=(sum+y*m*r[i])%n;
    }
    return (n+sum%n)%n;
}

一个正整数K，给出K Mod 一些质数的结果，求符合条件的最小的K。例如，K % 2 = 1, K % 3 = 2, K % 5 = 3。符合条件的最小的K = 23。

Input

第1行：1个数N表示后面输入的质数及模的数量。（2 <= N <= 10)
第2 - N + 1行，每行2个数P和M，中间用空格分隔，P是质数，M是K % P的结果。（2 <= P <= 100, 0 <= K < P)

Output

输出符合条件的最小的K。数据中所有K均小于10^9。

Sample Input

3
2 1
3 2
5 3

Sample Output

23

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
using namespace std;
typedef long long ll;
ll extend_Euclid(ll a,ll b,ll &x,ll &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return x;
    }
    extend_Euclid(b,a%b,x,y);
    ll t;
    t=x;
    x=y;
    y=t-a/b*y;
    return x;
}
ll Chinese_Remainder(ll r[],ll prime[],ll len)
{
    ll i,x,y,m,n=1,sum=0;
    for(i=0;i<len;i++)
        n*=prime[i];
    for(i=0;i<len;i++)
    {
        m=n/prime[i];
        extend_Euclid(m,prime[i],x,y);
        sum=(sum+x*m*r[i])%n;
    }
    return (n+sum%n)%n;
}
int main()
{
    ll prime[15],r[15],i,n;
    while(~scanf("%lld",&n))
    {
        for(i=0;i<n;i++)
        {
            scanf("%lld %lld",&prime[i],&r[i]);
        }
        printf("%lld\n\n",Chinese_Remainder(r,prime,n));
    }
    return 0;
}