51nod 1079中國剩餘定理
阿新 • • 發佈:2019-02-20
基準時間限制:1 秒 空間限制:131072 KB 分值: 0 難度:基礎題 收藏 關注一個正整數K,給出K Mod 一些質數的結果,求符合條件的最小的K。例如,K % 2 = 1, K % 3 = 2, K % 5 = 3。符合條件的最小的K = 23。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。Input示例
3 2 1 3 2 5 3Output示例
23
中國剩餘定理是用來解決模線性方程組問題的
求出 x 滿足這個等式
第一種方法:
直接寫
#include<cstdio> #include<cstring> #include<algorithm> #include<cmath> using namespace std; typedef long long ll; const int maxn=11; int main() { int n; while(~scanf("%d",&n)) { ll a[maxn],m[maxn]; ll lcm=1; ll ans=0; for(int i=0;i<n;i++) { scanf("%lld %lld",&a[i],&m[i]); lcm*=a[i]; } for(int i=0;i<n;i++) { int j=1; ll M=lcm/a[i]; if(M%a[i]!=0) for(j=1;j*M%a[i]!=1;j++); ans+=M*m[i]*j; } printf("%lld\n",ans%lcm); } return 0; }
第二種,利用exgcd
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<cmath> using namespace std; typedef long long ll; const int maxn=105; int n; ll m[maxn],p[maxn]; ll exgcd(ll a,ll b,ll &x,ll &y) { if(!b) { x=1; y=0; return a; } int r=exgcd(b,a%b,y,x); y-=x*(a/b); return r; } ll CRT(ll *m,ll *a) { ll M=1,ans=0; for(int i=0;i<n;i++) M*=m[i]; for(int i=0;i<n;i++) { ll mi=M/m[i],x,y; exgcd(mi,m[i],x,y); //要注意為止不能放反。。。 ans=(ans+x*mi*a[i])%M; } return (ans+M)%M; } int main() { while(~scanf("%d",&n)) { for(int i=0;i<n;i++) scanf("%lld %lld",m+i,p+i); printf("%lld\n",CRT(m,p)); } return 0; }