中国剩余定理及其扩展学习笔记

最新推荐文章于 2023-08-20 10:47:47 发布
最新推荐文章于 2023-08-20 10:47:47 发布
阅读量110 收藏
点赞数
原文链接：http://www.cnblogs.com/VictoryCzt/p/10053429.html
版权

Chinese Remainder Theorem(CRT)

又称孙子定理，是用来解决一类同余方程组问题的方法。

这个要用到欧几里得算法相关知识，不会走这里GO

我们从一个简单的例子入手：

现有如下式子（p1,p2p1,p2为质数）：

①{x≡a1(mod p1)x≡a2(mod p2)①{x≡a1(mod p1)x≡a2(mod p2)

按照同余式，上式可以写成：

②{x=a1+k1p1x=a2+k2p2②{x=a1+k1p1x=a2+k2p2其中k1,k2∈Zk1,k2∈Z

通过联立起来，将其写成a1+k1p1=a2+k2p2a1+k1p1=a2+k2p2,k1p1−k2p2=a2−a1k1p1−k2p2=a2−a1

又因为gcd(p1,p2)=1gcd(p1,p2)=1,所以此方程肯定有解。

对应ax+by=cax+by=c式子，那么⎧⎩⎨a=p1b=−p2c=a2−a1{a=p1b=−p2c=a2−a1，所以根据这种形式，我们用exgcdexgcd可以求出一组解(k^1,k^2)(k^1,k^2),然后我们将其带回原方程组②②可得x0x0，然后通过exgcdexgcd的知识我们可以由这组(k^1,k^2)(k^1,k^2)特解，推导出其它的解：{k1=p2t+k^1k2=p1t+k^2,k1,k2∈Z{k1=p2t+k^1k2=p1t+k^2,k1,k2∈Z。

带入原方程求解得：x=a1+k1p1=a1+p1p2t+k^1=x0+p1p2tx=a1+k1p1=a1+p1p2t+k^1=x0+p1p2t等价于x≡x0(mod p1p2)x≡x0(mod p1p2)。
我们就得到了方程的答案，同时将模数合并了起来。

所以对于多个方程构成的方程组，我们每次选择两个合并出一个新的，然后再用新的去和另一个合并（注意：模数相乘后小心爆范围long longlong long）。

那么对于如下的模方程组：

⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪x≡a1 (mod m1)x≡a2 (mod m2)x≡a3 (mod m3)⋯x≡an (mod mn){x≡a1 (mod m1)x≡a2 (mod m2)x≡a3 (mod m3)⋯x≡an (mod mn)

所有的mimi均互质，所以用exgcdexgcd求出特解即可解决。（LRJ的蓝书上也有讲解）。

模板代码：

\\[TJOI2009猜数字]
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int M=21;
int n;
ll a[M],b[M];
void exgcd(ll a,ll b,ll &x,ll &y){
  if(!b){x=1;y=0;}
  else {exgcd(b,a%b,y,x);y-=x*(a/b);}
}
ll mul(ll a,ll b,ll mod){
  ll c=0;
  for(;b>0;b>>=1,a=(a+a)%mod)if(b&1) c=(c+a)%mod;
  return c;
}
ll CRT(ll *a,ll *b,int n){
  ll M=1,d=0,x,y;
  for(int i=1;i<=n;i++) M*=b[i];
    for(int i=1;i<=n;i++){
      ll w=M/b[i];
      exgcd(b[i],w,x,y);y=(y%b[i]+b[i])%b[i];//exgcd求解
      d=(d+mul(mul(w,y,M),(a[i]%b[i]+b[i])%b[i],M))%M;
    }
    return (d+M)%M;
}
int main(){
  scanf("%d",&n);
  for(int i=1;i<=n;i++)scanf("%lld",&a[i]);
  for(int i=1;i<=n;i++)scanf("%lld",&b[i]);
  ll ans=CRT(a,b,n);
  printf("%lld\n",ans);
  return 0;
}

扩展中国剩余定理

我们现在解决的一类模方程的求法，但是这是在模数互质的条件下。大多数时候模数都不会互质，那么应该如何求解呢？
我们仍然从简单的入手

看如下一个方程组：

{x≡a1 (mod m1)x≡a2 (mod m2){x≡a1 (mod m1)x≡a2 (mod m2)

同样的我们将其变形得到 a1+k1m1=a2+k2m2a1+k1m1=a2+k2m2，如果该方程要有整数解，那么根据裴蜀定理，要满足 gcd(m1,m2)|(a2−a1)gcd(m1,m2)|(a2−a1)。

这里我们先假设用exgcdexgcd求出了一组解特解(k1,k2)(k1,k2)，那么通解有是什么呢？
接下来令g=gcd(m1,m2)g=gcd(m1,m2)，又因为g|(a2−a1)g|(a2−a1)（因为我们先假设它有解），所以原方程可以写成

k1m1g−k2m2g=a2−a1gk1m1g−k2m2g=a2−a1g
，这样我们就可得知 gcd(m1g,m2g)=1gcd(m1g,m2g)=1，系数就转换为之前的互质的情况了，互质后我们再通过之前的方法得到一组特解 (k^1,k^2)(k^1,k^2)，通解便为：
{k1=m2gt+k^1k2=m2gt+k^2t∈Z{k1=m2gt+k^1k2=m2gt+k^2t∈Z

回带入原方程，我们就可以得到:
x=a1+k1m1x=a1+k1m1

=x0+m1m2gt=x0+m1m2gt

=x0+lcm(m1,m2)t=x0+lcm(m1,m2)t

lcm为最小公倍数。

实质上就等于解：x≡x0 (mod lcm(m1,m2))x≡x0 (mod lcm(m1,m2))这个方程。

举个例子，解如下方程：

⎧⎩⎨①x≡2 (mod 3)②x≡3 (mod 5)③x≡2 (mod 7){①x≡2 (mod 3)②x≡3 (mod 5)③x≡2 (mod 7)

先自己试着解一下再看下面解答。

先联立①②①②得④x≡ 8(mod 15)④x≡ 8(mod 15)，再联立③④③④得到最终答案x≡23 (mod 105)x≡23 (mod 105)，那么答案就为x=23x=23，反带回去发现确实成立。

【模板IN】

模板代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int M=201;
ll xx(ll a,ll b,ll mod){
    if(a>b) swap(a,b);ll c=0;
    for(ll i=0;(1ll<<i)<=b;i++){if(b&(1ll<<i)){c=(c+(1ll<<i)*a)%mod;}}
    return c;
}
void exgcd(ll a,ll b,ll &d,ll &x,ll &y){
    if(!b){d=a;x=1;y=0;}
    else {exgcd(b,a%b,d,y,x);y-=x*(a/b);}
}
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
ll inv(ll a,ll b){
    ll g,x,y;
    exgcd(a,b,g,x,y);
    if(g!=1) return -1;
    return (x%b+b)%b;
}
bool merge(ll a1,ll b1,ll a2,ll b2,ll &a,ll &b){
    ll g=gcd(b1,b2);
    ll c=a2-a1;
    if(c%g) return 0;//判断有无解
    c=(c%b2+b2)%b2;
    b1/=g;b2/=g;
    c/=g;c=xx(c,inv(b1,b2),b2);
    b=b1*b2*g;
    c=xx(c,xx(b1,g,b),b);
    c+=a1;
    a=(c%b+b)%b;
    return 1;
}
ll CRT(ll *a,ll *b,int n){
    ll aa=a[1],bb=b[1];
    for(int i=2;i<=n;i++){
        ll da=a[i],db=b[i];
        ll nexa,nexb;
        //合并
        if(!merge(aa,bb,da,db,nexa,nexb)) return -1;
        aa=nexa;bb=nexb;
    }
    return (aa%bb+bb)%bb;
}
ll a[M],b[M];
int n;
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){scanf("%lld%lld",&b[i],&a[i]);}
    ll ans=CRT(a,b,n);
    printf("%lld\n",ans);
    return 0;
}
//模板2
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define LL long long
using namespace std;
const int N=210;
int n;
LL x,y,lcm,equ[N][2];
LL multi(LL a,LL b,LL p)
{
    a=(a%p+p)%p;b=(b%p+p)%p;LL ans=0;
    for(;a;a>>=1,b=(b*2)%p)if(a&1)ans=(ans+b)%p;
    return ans;
}
LL exgcd(LL a,LL b,LL &x,LL &y)
{
    if(!b)
    {
        x=1,y=0;
        return a;
    }
    LL val=exgcd(b,a%b,x,y);
    LL t=x;x=y;y=t-a/b*y;return val;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
     scanf("%lld%lld",&equ[i][1],&equ[i][0]);
    for(int i=1;i<n;i++)
    {
        LL val=exgcd(equ[i][1],equ[i+1][1],x,y);
        lcm=equ[i][1]/val*equ[i+1][1];
        if((equ[i+1][0]-equ[i][0])%val)
        {
            printf("No Answer\n");
            return 0;
        }
        val=multi(y,(equ[i+1][0]-equ[i][0])/val,lcm);
        equ[i+1][0]=(multi(equ[i+1][1],-val,lcm)+equ[i+1][0]+lcm)%lcm;
        equ[i+1][1]=lcm;
    }
    printf("%lld\n",(equ[n][0]%equ[n][1]+equ[n][1])%equ[n][1]);
    return 0;
}

一个额外的小技巧：
例如NOI2018的屠龙勇士。
当同余方程组有系数，如下式子：

⎧⎩⎨⎪⎪⎪⎪⎪⎪a1x≡b1 (mod m1)a2x≡b2 (mod m2)⋯anx≡bn(mod mn){a1x≡b1 (mod m1)a2x≡b2 (mod m2)⋯anx≡bn(mod mn)

如果a,ba,b都互质，我们用扩展欧几里得定理就可以求得逆元，然后把系数除过去就可以计算了。

但是当a,ba,b不一定互质怎么办？

我们观察式子，ax≡b(mod m)ax≡b(mod m),可以将其写成ax=km+bax=km+b,我们记g=gcd(a,m)g=gcd(a,m)(gcd表示求最大公约数)，那么左右可以写成(gw1)x=k(gw2)+b(gw1)x=k(gw2)+b,由于左边为gg的倍数，右边那么同理也应该为gg的倍数，所以显然bb也应该为gg的倍数，若不为则该方程组无整数解，有的话将a,b,ma,b,m同时除以gg，就得到了系数互质的方程，用原来的方法解就可以了。

End

Thanks For Reading!

转载于:https://www.cnblogs.com/VictoryCzt/p/10053429.html

weixin_30360497
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
[数论]扩展中国剩余定理（EX-CRT）
zero_orez6的博客
05-06 165
前言 CRT和EX-CRT目的都是为了解决同余的问题，但他们在方法上有着本质的区别，这篇文章主要用数学归纳法来证明EX-CRT。 问题 给你nnn个同余式子 {x≡a1(mod m1)x≡a2(mod m2)......x≡an(mod mn)\begin{cases}x\equiv a_1(mod\ m_1)\\x\equiv a_2(mod\ m_2)\\......\\x\equiv a_n(mod\ m_n)\end{cases}⎩⎪⎪⎪⎨⎪⎪⎪⎧​x≡a1​(mod&
中国剩余定理（孙子定理）详解
热门推荐
彼岸花苏陌的博客
09-09 2万+
题目引出 在《孙子算经》中有这样一个问题：“今有物不知其数，三三数之剩二（除以3余2），五五数之剩三（除以5余3），七七数之剩二（除以7余2），问物几何？”这个问题称为“孙子问题”，该问题的一般解法国际上称为“中国剩余定理”。 解决方法 1.找出三个数：从3和5的公倍数中找出被7除余1的最小数15，从3和7的公倍数中找出被5除余1 的最小数21，最后从5和7的公倍数中找出除3余1的最小数70。 2.用15乘以2（2为最终结果除以7的余数），用21乘以3（3为最终结果除以5的余数），同理，用70乘以2（2为
中国剩余定理及扩展
最新发布
weixin_74088105的博客
08-20 558
所以，孙子问题解法的本质是从5和7的公倍数中找一个除以3余2的数n1，从3和7的公倍数中找一个除以5余3的数n2，从3和5的公倍数中找一个除以7余2的数n3，再将三个数相加得到解。这就牵涉到一个最基本数学定理，如果有a%b=c，则有(a+k∗b)%b=c(k为非零整数)，换句话说，如果一个除法运算的余数为c，那么被除数与kk倍的除数相加（或相减）的和（差）再与除数相除，余数不变。首先，我们假设n1是满足除以3余2的一个数，比如2，5，8等等，也就是满足3∗k+2（k>=0）的一个任意数。
中国剩余定理(孙子定理)详解
aixiaqin3479的博客
02-11 1206
问题：今有物不知其数，三三数之剩二，五五数之剩三，七七数之剩二。问物几何？ 简单点说就是，存在一个数x，除以3余2，除以5余三，除以7余二，然后求这个数。上面给出了解法。再明白这个解法的原理之前，需要先知道一下两个定理。 定理1：两个数相加，如果存在一个加数，不能被整数a整除，那么它们的和，就不能被整数a整除。 定理2：两数不能整除，若除数扩大（或缩小）了几倍，而被除数不变，则其商...
中国剩余定理详解
jk13171217的专栏
07-29 4133
中国剩余定理 孙子算经里有这样一个问题：“今有物不知其数，三三数之剩二，五五数之剩三，七七数之剩二，问物几何？” 翻译成现在的数学问题就是 x%3 == 2，x%5 == 3，x%7 ==  2,求x 的值； 遇到这这样一个问题很多C语言初学者不禁会想到用暴力可以算出来，还要这样一个定理干嘛？ 如果数据相当大呢？计算机就会计算相当困难。然而这个问题早早的就被孙子解决了。  求出3,5,7
中国剩余定理算法详解（互质与不互质情况）
codeswarrior的博客
07-15 1万+
参考 中国剩余定理： 参考博客 中国剩余定理介绍 在《孙子算经》中有这样一个问题：“今有物不知其数，三三数之剩二（除以3 余2），五五数之剩三（除以5 余3），七七数之剩二（除以7 余2），问物几何？”这个问题称为“孙子问题”，该问题的一般解法国际上称为“中国剩余定理”。 具体解法分三步： 找出三个数，从3和5的公倍数中找出被7除余1（≡≡\equiv1 mod 7）的最小...
秘密共享方案_中国剩余定理
11-15
#某电信安数学基础实验（3）基于中国剩余定理的秘密共享方案（更新）： #涉及位数的均为二进制位 #python #使用中国剩余定理 #使用（t,n）门限来控制秘密，即：n个子秘密中任取t个或以上即可计算得到秘密，而任意t-1...
C语言实现中国剩余定理
01-21
计算一次同余方程的解，并且调用miracl库函数，实现大数解的计算。
Python实现的中国剩余定理算法示例
09-21
主要介绍了Python实现的中国剩余定理算法,结合实例形式分析了中国剩余定理的概念、原理及具体算法实现技巧,需要的朋友可以参考下
2_中国剩余定理_信息安全数学基础_
09-30
中国剩余定理（Chinese Remainder Theorem，CRT）是数论中的一个重要概念，它在密码学，特别是信息安全...通过对这些代码的分析和学习，我们可以更好地理解和应用中国剩余定理，提升在信息安全领域的理论和技术能力。
CQOJ 扩展中国剩余定理 【数论】
baicang2022的博客
07-11 572
原题： 　　求在小于等于 n 的正整数中有多少个X满足： 　　X = b[1] (mod a[1]) 　　X = b[2] (mod a[2]) 　　..... 题解： 　　对于中国剩余定理的理解： 　　　　对于前两式 　　　　　　X=B[1] (mod A[1])　 => A[1]*a+B[1]=X 　　　　　　X=B[2] (mod A[2])　 =...
【JZOJ5796】【2018提高组】模拟A组&省选 划分（容斥+扩展中国剩余定理+扩展欧几里得）
qq_36551189的博客
08-10 223
Problem   有一个未知的序列x，长度为n。它的K-划分序列y指的是每连续K个数的和得到划分序列，y[1]=x[1]+x[2]+….+x[K]，y[2]=x[K+1]+x[K+2]+….+x[K+K]….。若n不被K整除，则y[n/K+1]可以由少于K个数加起来。比如n=13，K=5，则y[1]=x[1]+…+x[5]，y[2]=x[6]+….+x[10]，y[3]=x[11]+x[12]...
中国剩余定理 && 扩展中国剩余定理 （模板）
辞树
10-19 578
中国剩余定理 &amp;&amp; 扩展中国剩余定理 一个整数除以三余二，除以五余三，除以七余二，求这个整数。 例题： 一个正整数K，给出K Mod 一些质数的结果，求符合条件的最小的K。例如，K % 2 = 1, K % 3 = 2, K % 5 = 3。符合条件的最小的K = 23。 #include&lt;bits/stdc++.h&gt; using namespace std;...
中国剩余定理总结
YouthDance
11-10 2320
中国剩余定理总结   当你遇到x == c (mod p) 要你求解x的时候，是不是很容易的想到了这样转换---> x - py = c运用extgcd得到答案。但是现在如果有很多个x呢？如：                                           x == c1(mod p1)                                      
中国剩余定理与扩展模板
KobeDuu
12-12 353
中国剩余定理： 解决模数互质的线性同余方程组： 中国剩余定理给出了这类方程的解法： （1）先求 X1 使得 X1 整除5、7并且 X1%3==2，再求 X2 使得 X2 整除3、7并且 X2%5==3，最后求 X3 使得 X3 整除3、5 并且 X3%7==2 （2）那么 X = X1 + X2 + X3 必然是上述同余方程组的一个解，验证第一个方程：X2%3==0，X3%3==0，...
中国剩余定理 扩展中国剩余定理
henulmh的博客
08-05 256
一、中国剩余定理 对于 x ≡ a1(mod m1) x ≡ a2(mod m2) x ≡ a3(mod m3) ………… 且 m1 m2 m3…两两互质 此时可运用中国剩余定理 首先求出 s=m1m2m3 M[i]=s/m[i] t[i]为M[i]关于m[i]的逆元 则最终求得的x =∑(a[i]*t[i]*M[i])%s 代码如下 #include<bits/stdc++.h> #...
扩展中国剩余定理讲解
weixin_30621711的博客
10-28 105
扩展中国剩余定理讲解 1.运用领域 　　扩展中国剩余定理是解决向下面列出的一元线性同余方程组的一种数论知识，可以求出下面方程组中最下的正整数$x$。但是扩展中国剩余定理和中国剩余定理有什么区别呢？中国剩余定理对于$mod$是有限制的，他对于$mod$要求为两两互质，然而扩展中国剩余定理对于$mod$没有要求。 　　$\begin{cases} x\equiv x_1 （mod_1）\\ x\...
中国剩余定理与扩展中国剩余定理
路人黑的纸巾
08-06 427
中国剩余定理（CRT） 我好蔡啊 不学这个东东我连任意模数NTTNTTNTT都学不了 问题 中国剩余定理用于求解同余方程组 {x≡a1(mod&ThinSpace;&ThinSpace;m1)x≡a2(mod&ThinSpace;&ThinSpace;m2)......x≡ak(mod&ThinSpace;&ThinSpace;mk) \left\{...
中国剩余定理及其应用解析
"这篇资源是加州州立大学圣伯纳迪诺分校的Jacquelyn Ha Lac的一篇论文，探讨了中国剩余定理及其应用。该论文属于代数领域，可以在CSUSB ScholarWorks平台上找到。作者通过这篇论文研究了中国剩余定理在数学，特别是...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值