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

[51nod1079]中国剩余定理

时间:2017-05-27 19:15:31      阅读:123      评论:0      收藏:0      [点我收藏+]

标签:for   cin   code   剩余定理   ace   typedef   nod   clu   its   

解题关键:注意爆long long

$x \equiv {M_1}M_1^{ - 1}{a_1} + ... + {M_k}M_k^{ - 1}{a_k}(\bmod m)$

其中,$m = \prod\limits_{j = 1}^k {{m_j}}$,$\forall 1 \le j \le k$,${M_j} = \frac{m}{{{m_j}}}$,$M_j^{ - 1}$是满足${M_j}M_j^{ - 1} \equiv 1(\bmod m)$的一个整数

复杂度$O(n\log n)$

 

 

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 typedef long long ll;
 4 ll a[100],b[100];
 5 ll x,y;
 6 ll extgcd(ll a,ll b,ll &x,ll &y){
 7     int d=a;
 8     if(b){
 9         d=extgcd(b,a%b,y,x);
10         y-=a/b*x;
11     }else{
12         x=1,y=0;
13     }
14     return d;
15 }
16 int main(){
17     ll n,m=1;
18     ll ans=0;
19     cin>>n;
20     for(int i=0;i<n;i++){ cin>>b[i]>>a[i];m*=b[i];}
21     for(int i=0;i<n;i++){
22         ll mi=m/b[i];
23         extgcd(mi,b[i],x,y);
24         x=(x+b[i])%b[i];
25         ans=(ans+mi*x*a[i]+m)%m;
26     }
27     cout<<ans<<endl;
28     return 0;
29 }

 

[51nod1079]中国剩余定理

标签:for   cin   code   剩余定理   ace   typedef   nod   clu   its   

原文地址:http://www.cnblogs.com/elpsycongroo/p/6914070.html

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