码迷,mamicode.com
首页 > 其他好文 > 详细

搞搞中国剩余定理和它的扩展

时间:2019-10-01 10:16:39      阅读:113      评论:0      收藏:0      [点我收藏+]

标签:扩展   crt   同余   gcd   exgcd   ++   情况   相关   ret   

中国剩余定理

解决同余方程组的相关问题,然而貌似扩展一下发生质变?

同余方程组问题

\[x\equiv a_1 \mod m_1 \]

\[x\equiv a_2 \mod m_2 \]

\[\dots\dots\]

\[x\equiv a_n \mod m_n \]

求解满足上述同余方程组的\(x\)

中国剩余定理(CRT)

  • 要求:模数\(m_1,m_2,\dots,m_n\) 必须互质(为什么呢?)
  • CRT告诉我们,这个方程组在模\(M=m_1*m_2*\dots*m_n\) 意义下的解是唯一的。
  • CRT还告诉我们这个解为\[x\equiv a_1 M_1 M_1^{-1}+a_2 M_2M_2^{-1}+\dots+a_n M_nM_n^{-1} \mod M\]其中\(M_i=M/m_i\) ,而\(M_i^{-1}\) 是模\(m_i\) 意义下的逆元。
  • \[x=\sum a_i\frac{M}{m_i}inv(\frac{M}{m_i},\,m_i)\]
  • 这个解是构造性的解,可以验证其是满足上述方程组的。
  • 为什么两两互质?逆元的要求!
const int maxn=1000;
int a[maxn];
int m[maxn];

int CRT(int n){
    int M=1;
    int ans=0;
    for(int i=1;i<=n;i++){
        M*=m[i];
    }
    for(int i=1;i<=n;i++){
        int x,y;
        int Mi=M/m[i];
        exgcd(Mi,m[i],x,y);
        ans=(ans+a[i]*Mi*x)%M;
    }
    if(ans<0)ans+=M;
    return ans;
}

扩展中国剩余定理(exCRT)

  • 模数\(m_1,m_2,\dots,m_n\) 不互质的时候也可以有解啦!
  • 考虑一下同余方程组仅有两个方程的情况。
  • \[x\equiv c_1 \mod m_1 \]\[x\equiv c_2 \mod m_2\]
  • \[x=c_1 +m_1k_1$ ,$x=c_2 +m_2k_2\]
  • \[c_1 +m_1k_1=c_2 +m_2k_2\]
  • \[m_1k_1 =c_2-c_1+m_2k_2\]
  • 这个方程有解的充要条件是\[\gcd(m_1,m_2)\,|\,(c_2-c_1)\]
  • \[\frac{m_1k_1}{\gcd(m_1,m_2)} =\frac{c_2-c_1}{\gcd(m_1,m_2)}+\frac{m_2k_2}{\gcd(m_1,m_2)}\]
  • \[\frac{m_1k_1}{\gcd(m_1,m_2)} \equiv\frac{c_2-c_1}{\gcd(m_1,m_2)}\mod \frac{m_2}{\gcd(m_1,m_2)}\]
  • \[k_1\equiv inv(\frac{m_1}{\gcd(m_1,m_2)},\frac{m_2}{\gcd(m_1,m_2)})*\frac{c_2-c_1}{\gcd(m_1,m_2)}\mod \frac{m_2}{\gcd(m_1,m_2)}\]
  • \(inv(a,b)\) 表示\(a\) 在模\(b\) 意义下的逆元
  • \[k_1=inv(\frac{m_1}{\gcd(m_1,m_2)},\frac{m_2}{\gcd(m_1,m_2)})*\frac{c_2-c_1}{\gcd(m_1,m_2)}+K* \frac{m_2}{\gcd(m_1,m_2)}\]
  • 带回\(x=c_1 +m_1k_1\)
  • \[x=c_1 +m_1*(inv(\frac{m_1}{\gcd(m_1,m_2)},\frac{m_2}{\gcd(m_1,m_2)})*\frac{c_2-c_1}{\gcd(m_1,m_2)}+K* \frac{m_2}{\gcd(m_1,m_2)})\]
  • \[x\equiv c_1 +m_1*inv(\frac{m_1}{\gcd(m_1,m_2)},\frac{m_2}{\gcd(m_1,m_2)})*\frac{c_2-c_1}{\gcd(m_1,m_2)}\mod \frac{m_1m_2}{\gcd(m_1,m_2)}\]
  • 到这里,我们已经完成了两个同余方程的整合!有形式\(x\equiv a \mod m\)
  • 其中\[a=c_1 +m_1*inv(\frac{m_1}{\gcd(m_1,m_2)},\frac{m_2}{\gcd(m_1,m_2)})*\frac{c_2-c_1}{\gcd(m_1,m_2)}\] \[m=\frac{m_1m_2}{\gcd(m_1,m_2)}\]
  • 反复重复以上操作!(二合一,加一,合一,······)
const int maxn=1000;
int a[maxn];
int m[maxn];

int gcd(int a,int b){
    return b?gcd(b,a%b):a;
}

int exgcd(int a,int b, int &x,int &y){
    if(b==0){
        x=1;
        y=0;
        return a;
    }
    int d=exgcd(b,a%b,x,y);
    int tmp=x;
    x=y;
    y=tmp-a/b*y;
    return d;
}

int inv(int a,int mod){
    int x,y;
    int d=exgcd(a,mod,x,y);
    return d==1?(x%mod+mod)%mod:-1;
}

int exCRT(int n){
    int m1, m2, c1, c2, d;
    for(int i=2;i<=n;i++){
        m1=m[i-1];//这里是与前一个进行合并
        m2=m[i];
        c1=a[i-1];
        c2=a[i];
        d=gcd(m1,m2);
        if((c2-c1)%d!=0){
            return -1;//无法合并
        }
        m[i]=m[i-1]*m[i]/d;//公式(倒数第二个)
        a[i]=c1+m1*inv(m1/d,m2/d)%(m2/d)*(c2-c1)/d;//公式(倒数第二个)
        a[i]=(a[i]%m[i]+m[i])%m[i];
    }
    return a[n];
}

搞搞中国剩余定理和它的扩展

标签:扩展   crt   同余   gcd   exgcd   ++   情况   相关   ret   

原文地址:https://www.cnblogs.com/tongseli/p/11614623.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!