1079 中國剩余定理
阿新 • • 發佈:2017-08-14
aps 關註 logs 數據 pre 結果 技術 fine 空間限制
基準時間限制:1 秒 空間限制:131072 KB 分值: 0 難度:基礎題
收藏
關註
一個正整數K,給出K Mod 一些質數的結果,求符合條件的最小的K。例如,K % 2 = 1, K % 3 = 2, K % 5 = 3。符合條件的最小的K = 23。
Input
Input示例
第1行:1個數N表示後面輸入的質數及模的數量。(2 <= N <= 10) 第2 - N + 1行,每行2個數P和M,中間用空格分隔,P是質數,M是K % P的結果。(2 <= P <= 100, 0 <= K < P)Output
輸出符合條件的最小的K。數據中所有K均小於10^9。
3 2 1 3 2 5 3Output示例
23
chuansong門
#include<iostream> #include<cstdio> using namespace std; #define ll long long int #define maxn 1000000 ll n,a[maxn],m[maxn]; ll exgcd(ll a,ll b,ll &x,ll &y) { if(b==0){ x=1; y=0; return a; } else { ll rView Code=exgcd(b,a%b,x,y); ll tmp=x; x=y; y=tmp-(a/b)*y; return r; } } ll crt() { ll M=1,mi=0,ans=0; for(ll i=1;i<=n;i++) M*=m[i]; for(ll i=1;i<=n;i++) { ll x=0,y=0; mi=M/m[i]; exgcd(mi,m[i],x,y); ans=(ans+a[i]*x*mi)%M; }if(ans<0) ans+=M; return ans; } int main() { scanf("%lld",&n); for(ll i=1;i<=n;i++) scanf("%lld%lld",&m[i],&a[i]); cout<<crt(); return 0; }
#include<iostream> #include<cstdio> using namespace std; #define ll long long int #define maxn 1000000 ll n,a[maxn],m[maxn]; ll exgcd(ll a,ll b,ll &x,ll &y) { if(b==0){ x=1; y=0; return a; } else { ll r=exgcd(b,a%b,x,y); ll tmp=x; x=y; y=tmp-(a/b)*y; return r; } } ll crt() { ll a1=a[1],a2,m2,c,d,m1=m[1]; for(int i=2;i<=n;++i) { ll x=0,y=0; a2=a[i],m2=m[i]; c=a2-a1; d=exgcd(m1,m2,x,y); ll mod=m2/d; if(c%d) return -1; x=x*c/d; x=(x%mod+mod)%mod; a1+=m1*x; m1*=mod; } if(a1==0) a1+=m1; return a1; } int main() { scanf("%lld",&n); for(ll i=1;i<=n;i++) scanf("%lld%lld",&m[i],&a[i]); cout<<crt(); return 0; }View Code
1079 中國剩余定理