首页 > 其他 > 详细

51nod1079中国剩余定理

时间:2016-03-30 22:08:26      阅读:159      评论:0      收藏:0      [点我收藏+]
/**
 *中国剩余定理
 */
#include<iostream>
#include<cstdio>
#include<map>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<vector>
#include<stack>
#include<cstdlib>
#include<cctype>
#include<cstring>
#include<cmath>
#define  LL __int64
using namespace std;
/**
 *gcd(a,b)=d;则存在x,y,使d=ax+by
 *extended_euclid(a,b)=ax+by
 */
LL extended_euclid(LL a,LL b,LL &x,LL &y){//扩张欧几里的算法
    int d;
    if(b==0){
        x=1; y=0;
        return a;
    }
    d=extended_euclid(b,a%b,y,x);
    y=y-a/b*x;
    return d;
}
/**
 *x=b[i](modw[i]) o<i<len
 *w[i]>0,且w[]中任意两个数互质
 */
LL chinese_remainder(int b[],int w[],int len){
    LL res,i,d,x,y,n,m;
    res=0; n=1;
    for(i=0;i<len;i++) n*=w[i];
    for(i=0;i<len;i++){
        m=n/w[i];
        extended_euclid(w[i],m,x,y);
        res=(res+y*m*b[i])%n;
    }
    return (n+res%n)%n;
}

int main()
{
    int len,b[12],w[12];
    while(cin>>len){
        for(int i=0;i<len;i++){
            cin>>w[i]>>b[i];
        }
        cout<<chinese_remainder(b,w,len)<<endl;

    }
    return 0;
}

 

51nod1079中国剩余定理

原文:http://www.cnblogs.com/OMG-By/p/5338838.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!