51Nod 1079 中国剩余定理 Label:数论
生活随笔
收集整理的這篇文章主要介紹了
51Nod 1079 中国剩余定理 Label:数论
小編覺(jué)得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
一個(gè)正整數(shù)K,給出K Mod 一些質(zhì)數(shù)的結(jié)果,求符合條件的最小的K。例如,K % 2 = 1, K % 3 = 2, K % 5 = 3。符合條件的最小的K = 23。 Input 第1行:1個(gè)數(shù)N表示后面輸入的質(zhì)數(shù)及模的數(shù)量。(2?<=?N?<=?10)
第2?-?N?+?1行,每行2個(gè)數(shù)P和M,中間用空格分隔,P是質(zhì)數(shù),M是K?%?P的結(jié)果。(2?<=?P?<=?100,?0?<=?K?<?P) Output 輸出符合條件的最小的K。數(shù)據(jù)中所有K均小于10^9。 Input示例 3
2?1
3?2
5?3 Output示例 23
代碼
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<algorithm> 5 #define ll long long 6 using namespace std; 7 8 ll x,y,n,a[15],m[15]; 9 ll exgcd(ll a,ll b){ 10 if(b==0){ 11 x=1;y=0;return a; 12 } 13 ll g=exgcd(b,a%b); 14 ll t=x; 15 x=y; 16 y=t-(a/b)*y; 17 return g; 18 } 19 20 21 ll CRT() 22 { 23 ll M = 1; 24 ll ans = 0; 25 for(ll i=1; i<=n; i++) 26 M *= m[i]; 27 for(ll i=1; i<=n; i++) 28 { 29 x=0, y=0; 30 ll Mi = M / m[i]; 31 exgcd(Mi, m[i]); 32 ans = (ans + Mi * x * a[i]) % M; 33 } 34 if(ans < 0) ans += M; 35 return ans; 36 } 37 38 int main(){ 39 // freopen("01.in","r",stdin); 40 41 scanf("%d",&n); 42 for(ll i=1;i<=n;i++){ 43 scanf("%d%d",&m[i],&a[i]); 44 } 45 cout<<CRT()<<endl; 46 return 0; 47 } 轉(zhuǎn)載代碼 不予展示請(qǐng)?jiān)谂蛘弑诚?擴(kuò)展歐幾里得算法,同余?的概念或代碼后來(lái)看這道題
建議先去百度百科逛一圈(我已經(jīng)放棄了wiki。。)
聽(tīng)大神講CRT http://blog.csdn.net/acdreamers/article/details/8050018
看大神證明CRT?http://blog.csdn.net/hard_man/article/details/7732795
提示,好好看~
轉(zhuǎn)載于:https://www.cnblogs.com/radiumlrb/p/5930831.html
總結(jié)
以上是生活随笔為你收集整理的51Nod 1079 中国剩余定理 Label:数论的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 几段表单处理的JQuery代码
- 下一篇: ignite通过注解配置查询