中国剩余定理

最新推荐文章于 2020-09-09 16:24:59 发布

G-Howe

最新推荐文章于 2020-09-09 16:24:59 发布

阅读量378

点赞数

分类专栏：中国剩余定理文章标签：中国剩余定理孙子定理

本文链接：https://blog.csdn.net/jk211766/article/details/81134685

版权

中国剩余定理专栏收录该内容

2 篇文章 0 订阅

订阅专栏

中国剩余定理详解见：https://www.cnblogs.com/linyujun/p/5199415.html

我是没怎麽理解这个定理，此题套用模板

一个正整数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

Sample Output

代码：

#include<algorithm>
using namespace std;
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<stdlib.h>
#include<stdio.h>
typedef long long int LL;
void ex_gcd(LL a,LL b,LL &x,LL &y)
{
    if(!b)
    {
        x=1;
        y=0;
    }
    else
    {
        ex_gcd(b,a%b,y,x);
        y-=x*(a/b);
    }
}
LL inv(LL m,LL n)
{
    LL x,y;
    ex_gcd(m,n,x,y);
    return x%n+n;
}
LL CRT(LL *n,LL *m,LL t)
{
    LL s=0,k=1,i;
    for(i=0;i<t;i++)
        k*=m[i];
    for(i=0;i<t;i++)
    {
        LL temp=k/m[i];
        s=(s+temp*inv(temp,m[i])*n[i])%k;
    }
    return (s+k)%k;
}
int main()
{
    LL m[20],n[20],k,i,s=0;
    while(~scanf("%lld",&k))
    {
        for(i=0;i<k;i++)
            scanf("%lld%lld",&m[i],&n[i]);
        s=CRT(n,m,k);
        printf("%lld\n",s);
    }
    return 0;
}