首页 > 其他分享 >扩展中国剩余定理(EXCRT)

扩展中国剩余定理(EXCRT)

时间:2023-06-15 17:04:10浏览次数:45  
标签:剩余 gcd dfrac 定理 km EXCRT mod 1m 式子

中国剩余定理(CRT)不能解决模数不互质情况的模线性同余方程组。这是中国剩余定理的原理所决定的。

但当我们的模数不互质时,这个方式显然就寄掉了,因此我们要打破原有的思路,去找一个新的方式解不定方程组,这时我们的扩展中国剩余定理(EXCRT)就出现了

假设我们现在有如下不定方程组

\[\begin{cases} x \equiv r_1(mod \ m_1) \\ x \equiv r_2(mod \ m_2) \\ x \equiv r_3 (mod \ m_3) \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot\end{cases} \]

显然我们可以得到 $ x = m_1k_1 + r_1 = m_2k_2 + r_2$ ---------- \((1)\)

则我们移项可得 $ m_1k_1 + m_2k_2 = r_2 - r_1$ ---------- \((2)\)

这里我们 \(m_2k_2\) 的系数变为正,是因为 \(k_2\) 是任意整数,所以前面的系数不会影响整个式子)

我们已知式子 \(m_1p_1 + m_2p_2 = \gcd(m_1,m_2)\) 将 \(d\) 设为 $ \gcd(m_1,m_2) $ 。

根据裴蜀定理,此时我们要进行判断,判断 \(r_2-r_1\) 是否为 \(\gcd(m_1,m_2)\) 的倍数,若是,那么有解,若不是,那么无解。

判断完是否有解后我们可以的到一个显然的式子 \(m_1p_1+m_2p_2+km_1m_2-km_1m_2=m_1p_1+m_2p_2=d=\gcd(m_1,m_2)\)

此时最左边的式子我们可以化简为 \(m_1(p_1+km_2)+m_2(p_2-km_1)=d\)

由于我们要找到原式子的通解,所以此时我们将等式两边同乘 $ \dfrac {r_2-r_1}{\gcd(m_1,m_2)} $

则可以得到式子 \(m_1(\dfrac {p_1(r_2-r_1)}{\gcd(m_1,m_2)} + \dfrac {km_2(r_2-r_1)}{\gcd(m_1,m_2)}) + m_2(\dfrac {p_2(r_2-r_1)}{\gcd(m_1,m_2)}- \dfrac {km_1(r_2-r_1)}{\gcd(m_1,m_2)}) = r_2-r_1\) ---------- \((3)\)

由上面的等式 \((2)\) 得可将 \((3)\) 式右面 \(r_2-r_1\) 的式子代换成 $ m_1k_1 + m_2k_2 $

从而得到的等式 $m_1(\dfrac {p_1(r_2-r_1)}{\gcd(m_1,m_2)} + \dfrac {km_2(r_2-r_1)}{\gcd(m_1,m_2)}) + m_2(\dfrac {p_2(r_2-r_1)}{\gcd(m_1,m_2)}- \dfrac {km_1(r_2-r_1)}{\gcd(m_1,m_2)}) = m_1k_1 + m_2k_2 $

由此我们显然可以得到

\[m_1k_1 = m_1(\dfrac {p_1(r_2-r_1)}{\gcd(m_1,m_2)} + \dfrac {km_2(r_2-r_1)}{\gcd(m_1,m_2)}) \]

\[m_2k_2 = m_2(\dfrac {p_2(r_2-r_1)}{\gcd(m_1,m_2)} - \dfrac {km_1(r_2-r_1)}{\gcd(m_1,m_2)}) \]

消掉第一个式子中的 \(m_1\) 。则我们求得 \(k_1=p_1\dfrac {(r_2-r_1)}{\gcd(m_1,m_2)} + km_2 \dfrac {(r_2-r_1)}{\gcd(m_1,m_2)}\) ---------- \((4)\)

我们将现在得到的等式 \((4)\) 带入等式 \((1)\) 则显然可以得到

\[x = m_1(p_1\dfrac {(r_2-r_1)}{\gcd(m_1,m_2)} + km_2 \dfrac {(r_2-r_1)}{\gcd(m_1,m_2)}) + r_1 \]

化简得

\[x = p_1m_1\dfrac {(r_2-r_1)}{\gcd(m_1,m_2)} + km_1m_2 \dfrac {(r_2-r_1)}{\gcd(m_1,m_2)} + r_1 \]

由于 \(k\) 是任意的整数,因此我们可以将 \(k\) 的系数化成 \(m_1m_2 \dfrac {1}{\gcd(m_1,m_2)}\)

从而得到最终的式子

\[x = p_1m_1\dfrac {(r_2-r_1)}{\gcd(m_1,m_2)} + km_1m_2 \dfrac {1}{\gcd(m_1,m_2)} + r_1 \]

我们此时显然是知道 \(m_1,m_2,r_1,r_2,p_1,\gcd(m_1,m_2)\) 的值,但我们不知道 \(k\) 的值,那么这时我们只要模上 \(k\) 的系数即可消掉 \(k\) 。即把模数变为 \(m_1m_2 \dfrac {(r_2-r_1)}{\gcd(m_1,m_2)}\)

此时我们就将两个不定方程组合并完成了

\(x \equiv p_1m_1\dfrac {(r_2-r_1)}{\gcd(m_1,m_2)} + r_1 (mod \ \dfrac {m_1m_2}{\gcd(m_1,m_2)})\)

则此时新的不定方程的 \(r\) 为 \(\dfrac {(r_2-r_1)}{\gcd(m_1,m_2)} + r_1\) , \(m\) 为 \(\dfrac {m_1m_2}{\gcd(m_1,m_2)}\)

再继续往下合并

最后我们可以合并成唯一一个式子 \(x \equiv r(mod \ m)\)

则 \(x\) 的最小整数解为 $ x = ( r \ mod \ m + m ) \ mod \ m $

( \(r \ mod \ m+m\) 中加 \(m\) 是为了防止出现 \(r\) 为负数的情况)

标签:剩余,gcd,dfrac,定理,km,EXCRT,mod,1m,式子
From: https://www.cnblogs.com/jueqingfeng/p/17483395.html

相关文章

  • 隐函数定理的几何应用
    隐函数定理的几何应用一、平面曲线的切线与法线设平面曲线由方程\[F(x,y)=0\tag{1}\]确定,它在\(P_0(x_0,y_0)\)的某领域上满足隐函数定理的条件,于是在点\(P_0\)附近所确定的连续可微隐函数\(y=f(x)\)(或\(x=g(y)\))和方程\((1)\)在\(P_0\)附近表示同一曲线,从而该曲......
  • Luogu P1495 【模板】中国剩余定理(CRT)/ 曹冲养猪
    【模板】中国剩余定理(CRT)/曹冲养猪题目描述自从曹冲搞定了大象以后,曹操就开始捉摸让儿子干些事业,于是派他到中原养猪场养猪,可是曹冲满不高兴,于是在工作中马马虎虎,有一次曹操想知道母猪的数量,于是曹冲想狠狠耍曹操一把。举个例子,假如有\(16\)头母猪,如果建了\(3\)个猪圈,剩下......
  • CAP原则(CAP定理)、BASE理论
    一、CAP原则 CAP原则又称CAP定理,指的是在一个分布式系统中,Consistency(一致性)、Availability(可用性)、Partitiontolerance(分区容错性),三者不可得兼。CAP原则是NOSQL数据库的基石。分布式系统的CAP理论:理论首先把分布式系统中的三个特性进行了如下归纳:一致性(C):在分......
  • 调题时出现的问题 in 『中国剩余定理』
    1(焯冲养pig/板子)【模板】中国剩余定理(CRT)/曹冲养猪要注意这东西不能用费马小定理,只能用扩欧.因为费马小定理的适用条件是模数为质数.......
  • 数学分析复习:Weierstrass 逼近定理, Müntz–Szász 定理
    本学期的“数学分析(不是实验班)”讲了一堆Approximationtheory,这是怎么绘事呢?定理1(Weierstrass).连续函数\(f\in\mathrmC[0,1]\)可被多项式一致逼近.对任意\(\varepsilon>0\)和\(x\in[0,1]\),设随机变量\(X\)服从二项分布\(\mathrmB(n,x)\),由Chebysh......
  • Lucas(卢卡斯定理)
    \(C^m_n\equivC^{m/p}_{n/p}*C^{m\mod\p}_{n\mod\p}\)首先,我们可以知道如下定理我们令\(n=ap+b\),\(m=cp+d\)则由二项式定理得\((1+x)^n\equiv\Sigma_{i=0}^nC^i_nx^i(mod\p)\)---------(1)由\(n=ap+b\)可知\((1+x)^n\equiv(1+x)^{......
  • 奈奎斯特采样定理中的奈奎斯特到底是谁?
    当用手机和家人通话、视频的时候,你有没有想过你的声音、影像为什么能传送到千里之外的地方? 这个问题还要从模拟信号说起。模拟信号是指用连续变化的物理量表示的信息。我们所在的世界中充满了各种各样的模拟信号,如声音、温度、湿度、压力、转速等等。可是这些连续的信号对于只知道......
  • 隐函数存在唯一性定理
    ......
  • k 倍区间(同余定理,组合数)
    题目描述给定一个长度为 N 的数列,1,2,⋯A1​,A2​,⋯AN​,如果其中一段连续的子序列 ,+1,⋯(≤)Ai​,Ai+1​,⋯Aj​(i≤j) 之和是 K 的倍数,我们就称这个区间 [,][i,j] 是 K 倍区间。你能求出数列中总共有多少个 K 倍区间吗?输入格式第一行包含两个整数 N 和 K......
  • 信道容量与香农定理、信源编码、信道编码总结
    1信道容量定义1.1信道容量:信道中平均每个符号所能传递的最大互信息量$I(X;Y)$$C=\mathop{max}\limits_{p(x)}{I(X;Y)}$单位:bit/符号1.2单位时间t内信道容量:$C_t=\frac{C}{t}$单位:bit/s1.3最佳输入概率$p(x)$分布时,传输的信息能达到信道容量1.4信道容量反映信道特性,表示信......