关于中国剩余定理的理解及模板

对于中国剩余定理的表述: 
设m1,m2,m3….mk两两互素,则同余方程组 
这里写图片描述

有整数解。并且在模M=m1*m2*……*mk下的解是唯一的,解为: 
这里写图片描述

这个方程十分关键,用来解题十分方便。
其中Mi=M/m[i],而Mi^-1为Mi模m[i]的逆元。用扩展欧几里得求出,就是求出方程x-k*m[i]=a[i]中的k。

故可以写程序求出x。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long LL;
LL p[20], r[20];
void ex_gcd(LL a, LL b, LL &x, LL &y){
    if(!b){
        x = 1;
        y = 0;
        return;
    }
    ex_gcd(b, a%b, x, y);
    LL temp = x;
    x = y;
    y = temp - a/b * y;
}
int main(){
    int T, n, kase = 0;
    scanf("%d", &T);
    while(T--){
        scanf("%d", &n);
        LL MOD = 1;
        for(int i = 0; i < n; i++){
            scanf("%lld%lld", &p[i], &r[i]);
            MOD *= p[i];
        }
        LL x, y;
        LL ans = 0;
        for(int i = 0; i < n; i++){
            ex_gcd(MOD/p[i], p[i], x, y);
            ans = (ans + MOD/p[i]*x*r[i] + MOD) % MOD;
        }
        printf("Case %d: %lld\n",++kase, (ans + MOD) % MOD);
    }
    return 0;
}
#include<bits/stdc++.h>
# define ll long long
using namespace std;
int extend_gcd(int a,int b,int &x,int &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    int gcd=extend_gcd(b,a%b,x,y);
    int temp=x;
    x=y;
    y=temp-a/b*y;
    return gcd;
}
int china(int w[],int b[],int k)//b[i]就相当于那个余数
{
    int x,y,a=0,m,n=1,i;
    for(i=0;i<k;i++)
      n*=w[i];
    for(i=0;i<k;i++)
        {
            m=n/w[i];
            extend_gcd(w[i],m,x,y);
            a=(a+y*m*b[i])%n;
        }
    if(a>0)return a;
    else return a+n;
}
int main()
{

}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值