中国剩余定理CRT (互质)

中国剩余定理给出了以下的一元线性同余方程组:
中国剩余定理CRT (互质)_第1张图片
有解的判定条件,并用构造法给出了在有解情况下解的具体形式。
中国剩余定理说明:假设整数m1,m2, … ,mn两两互质,则对任意的整数:a1,a2, … ,an,方程组 有解,并且通解可以用如下方式构造得到:
设M是整数m1,m2, … ,mn的乘积,并设Mi是除了mi以外的n- 1个整数的乘积。
设为Mi^-1模mi的数论倒数 :
方程组的通解形式为 这里写图片描述:在模M的意义下,方程组 这里写图片描述只有一个解:
例题 POJ 1006
51nod 1079

#include 
#include 
#include 
using namespace std;
typedef long long ll ;
ll exgcd(ll a ,ll b , ll &x ,ll &y ) {
    if(!b) {
        x = 1 , y = 0 ;
        return a ;
    }
    else {
        ll d = exgcd(b,a%b,x,y) ;
        ll t = x ;
        x = y ;
        y = t - a/b*y ;
        return d ;
    }
}
int main() {
    ll n ;
    while(~scanf("%lld",&n)){
            vector p,m,a;
            p.clear() ; m.clear() ; a.clear() ;
            long long M = 1 , ans = 0 ;
        for(ll pi,ai,i=0 ; iscanf("%lld%lld",&pi,&ai) ;
            p.push_back(pi) ;
            a.push_back(ai) ;
            M *= pi ;
        }
        for(ll i=0 ; iif(x<0)x+=p[i] ;
           // cout << x << " " << y <
            ans = (ans+a[i]*x*m[i])%M;
            //cout << ans << endl;
        }
//        ans = (ans+M)%M ;
        std::cout<std::endl;
    }
    return 0 ;
}

你可能感兴趣的:(ACM-数学)