1. 程式人生 > >HDU——1788 Chinese remainder theorem again

HDU——1788 Chinese remainder theorem again

ron %d pri strong 結束 str 最小 include 問題

        Chinese remainder theorem again

        Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
            Total Submission(s): 3174 Accepted Submission(s): 1371

Problem Description 我知道部分同學最近在看中國剩余定理,就這個定理本身,還是比較簡單的:
假設m1,m2,…,mk兩兩互素,則下面同余方程組:
x≡a1(mod m1)
x≡a2(mod m2)

x≡ak(mod mk)
在0<=<m1m2…mk內有唯一解。
記Mi=M/mi(1<=i<=k),因為(Mi,mi)=1,故有二個整數pi,qi滿足Mipi+miqi=1,如果記ei=Mi/pi,那麽會有:
ei≡0(mod mj),j!=i
ei≡1(mod mj),j=i
很顯然,e1a1+e2a2+…+ekak就是方程組的一個解,這個解加減M的整數倍後就可以得到最小非負整數解。
這就是中國剩余定理及其求解過程。
現在有一個問題是這樣的:
一個正整數N除以M1余(M1 - a),除以M2余(M2-a), 除以M3余(M3-a),總之, 除以MI余(MI-a),其中(a<Mi<100 i=1,2,…I),求滿足條件的最小的數。 Input 輸入數據包含多組測試實例,每個實例的第一行是兩個整數I(1<I<10)和a,其中,I表示M的個數,a的含義如上所述,緊接著的一行是I個整數M1,M1...MI,I=0 並且a=0結束輸入,不處理。 Output 對於每個測試實例,請在一行內輸出滿足條件的最小的數。每個實例的輸出占一行。 Sample Input 2 1 2 3 0 0 Sample Output 5 Author lcy Source 2007省賽集訓隊練習賽(10)_以此感謝DOOMIII Recommend lcy | We have carefully selected several similar problems for you: 1695 1452 1060 1370 1905 題目一樣的裸題、、、、(中國剩余定理裸題) 代碼:
#include<cstdio>
#include
<cstring> #include<cstdlib> #include<iostream> #include<algorithm> #define N 20 #define ll long long using namespace std; ll n,m[N],a[N],m1,e; ll read() { ll x=0,f=1; char ch=getchar(); while(ch<0||ch>9){if(ch==-)f=-1; ch=getchar();} while(ch>=
0&&ch<=9){x=x*10+ch-0; ch=getchar();} return x*f; } ll exgcd(ll a,ll b,ll &x,ll &y) { if(b==0) { x=1,y=0; return a; } ll r=exgcd(b,a%b,x,y),tmp; tmp=x,x=y,y=tmp-a/b*y; return r; } ll crt() { ll a1=a[1],a2,m2,d,c;m1=m[1]; for(ll i=2;i<=n;++i) { a2=a[i],m2=m[i]; c=a2-a1;ll x=0,y=0; d=exgcd(m1,m2,x,y); if(c%d) return -1; x=x*c/d; int mod=m2/d; x=(mod+x%mod)%mod; a1+=m1*x;m1*=mod; } return a1; } int main() { while(1) { n=read(),e=read(); if(n==0&&e==0) break; for(int i=1;i<=n;i++) m[i]=read(),a[i]=m[i]-e; printf("%lld\n",crt()); } return 0; }

HDU——1788 Chinese remainder theorem again