码迷,mamicode.com
首页 > 其他好文 > 详细

【Vijos】【1164】曹冲养猪

时间:2015-02-03 00:40:41      阅读:392      评论:0      收藏:0      [点我收藏+]

标签:

中国剩余定理

  没啥重要的……模板题,中国剩余定理就是解出模线性方程组的一个可行解(好像也是唯一解?)

  这是一种神奇的构造方法……明白了为什么这样构造是对的就行了=。=至于怎么想到这种构造方法的……去问孙子去→_→

技术分享
 1 //Vijos 1164
 2 #include<cstdio>
 3 #include<cstdlib>
 4 #include<cstring>
 5 #include<iostream>
 6 #include<algorithm>
 7 #define rep(i,n) for(int i=0;i<n;++i)
 8 #define F(i,j,n) for(int i=j;i<=n;++i)
 9 #define D(i,j,n) for(int i=j;i>=n;--i)
10 using namespace std;
11 
12 int getint(){
13     int v=0,sign=1; char ch=getchar();
14     while(ch<0||ch>9) {if (ch==-) sign=-1; ch=getchar();}
15     while(ch>=0&&ch<=9) {v=v*10+ch-0; ch=getchar();}
16     return v*=sign;
17 }
18 typedef long long LL;
19 /*******************tamplate********************/
20 
21 void exgcd(LL a,LL b,LL &d,LL &x,LL &y){
22     if(!b){ d=a; x=1; y=0; }
23     else { exgcd(b,a%b,d,y,x); y-=x*(a/b); }
24 }
25 LL china(int n,int *a,int *m){
26     LL M=1,d,y,x=0;
27     F(i,1,n) M*=m[i];
28     F(i,1,n){
29         LL w=M/m[i];
30         exgcd(m[i],w,d,d,y);
31         x=(x+y*w*a[i])%M;
32     }
33     return (x+M)%M;
34 }
35 int a[11],b[11];
36 int main(){
37     int n=getint();
38     F(i,1,n) a[i]=getint(),b[i]=getint();
39     printf("%lld\n",china(n,b,a));
40     return 0;
41 }
View Code

 

【Vijos】【1164】曹冲养猪

标签:

原文地址:http://www.cnblogs.com/Tunix/p/4269016.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!