1. 程式人生 > >1079 中國剩余定理

1079 中國剩余定理

aps 關註 logs 數據 pre 結果 技術 fine 空間限制

基準時間限制: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 3
Output示例
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 r
=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; }
View Code

技術分享
#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 中國剩余定理