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

扩展中国剩余定理

时间:2023-05-03 10:13:05浏览次数:46  
标签:剩余 frac gcd pmod 定理 扩展 cases inv equiv

前置知识

1.乘法逆元
2.朴素欧几里得

问题

已知\(\begin{cases}x\equiv c_1\pmod{m_1}\\x\equiv c_2\pmod{m_2}\\x\equiv c_3\pmod{m_3}\\...\\x\equiv c_n\pmod{m_n}\end{cases}\),求未知数\(x\)

求解

若能将这两个式子\(\begin{cases}x\equiv c_1\pmod{m_1}\\x\equiv c_2\pmod{m_2}\end{cases}\)合并为一个,则可以递推合并整个式子,故问题转化为求解

\[\begin{cases}x\equiv c_1\pmod{m_1}\\x\equiv c_2\pmod{m_2}\end{cases} \]

可以转化为$$c_1+m_1k_1=c_2+m_2k_2$$
两边同时除以\(\gcd(m_1,m_2)\)得

\[\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)}\pmod{\frac{m_2}{\gcd(m_1,m_2)}} \]

两边同时除以\(\frac{m_1}{\gcd(m_1,m_2)}\)得

\[k_1=\frac{c_2-c_1}{\gcd(m_1,m_2)}*inv(\frac{m_1}{\gcd(m_1,m_2)})\pmod{\frac{m_2}{\gcd(m_1,m_2)}} \]

\[k_1=\frac{c_2-c_1}{\gcd(m_1,m_2)}*inv(\frac{m_1}{\gcd(m_1,m_2)})+{\frac{m_2}{\gcd(m_1,m_2)}}*y \]

代入\(x=k_1m_1+c_1\)得

\[x=\frac{c_2-c_1}{\gcd(m_1,m_2)}*inv(\frac{m_1}{\gcd(m_1,m_2)})*m_1+c_1+{\frac{m_1m_2}{\gcd(m_1,m_2)}}*y \]

\[x\equiv{{\frac{c_2-c_1}{\gcd(m_1,m_2)}*inv(\frac{m_1}{\gcd(m_1,m_2)})*m_1+c_1}\pmod{{\frac{m_1m_2}{\gcd(m_1,m_2)}}}} \]

\[c_1'={\frac{c_2-c_1}{\gcd(m_1,m_2)}*inv(\frac{m_1}{\gcd(m_1,m_2)})*m_1+c_1} \]

\[m_1'={\frac{m_1m_2}{\gcd(m_1,m_2)}} \]

不断递推直到只剩下一个式子\(x\equiv{c'\pmod{m’}}\),得\(x=c'+km'\)

标签:剩余,frac,gcd,pmod,定理,扩展,cases,inv,equiv
From: https://www.cnblogs.com/cdx1221/p/17368713.html

相关文章

  • 浅聊Java核心技术之高可扩展利器SPI
    SPI的概念JAVASPI =基于接口的编程+策略模式+配置文件的动态加载机制SPI的使用场景Java是一种面向对象语言,虽然Java8开始支持函数式编程和Stream,但是总体来说,还是面向对象的语言。在使用Java进行面向对象开发时,一般会推荐使用基于接口的编程,程序的模块与模块之前不会直接进行......
  • m基于matlab的AODV,leach自组网网络平台仿真,对比吞吐量,端到端时延,丢包率,剩余节点
    1.算法仿真效果matlab2022a仿真结果如下:     2.算法涉及理论知识概要       AODV是一种应用于无线网状网络的路由协议。它源节点需要发送数据时才进行路由发现。当没有数据发送请求时并不执行。在路由发现过程中首先检查路由表中是否存在从源节点到目的......
  • Dockers下php容器中安装redis扩展
    首先进入php容器dockerexec-it容器ID或名称查看php安装位置  whichphp查看php已安装扩展  php-m1、下载redis扩展包   redis扩展下载地址【https://pecl.php.net/package/redis 】下载相应版本的扩展2、解压扩展包   tar-zxvfredis-5.1.1.tg......
  • 中国剩余定理
    中国剩余定理:代码实现://互质版中国剩余定理(CRT)#include<iostream>usingnamespacestd;typedeflonglongLL;constintN=20;LLa[N],b[N];intn;voidexgcd(LLa,LLb,LL&x,LL&y){if(!b){x=1,y=0;return;}exgcd(b,......
  • 浅谈裴蜀定理
    前置知识扩展欧几里得问题给定$a,b,$设$s=ax+by$,求当$s>$0时,求s的最小值定理$\min(s)=\gcd(a,b)$证明见扩展欧几里得引理给定n个数,分别为$A_1$,$A_2$,$A_3$...$A_n$任取n个数,分别为$X_1$,$X_2$,$X_3$...$X_n$设$s=\sum_{i=1}^NA_i*X_i$使$s>0且使s$......
  • 浅谈扩展欧几里得
    前置知识朴素欧几里得问题已知$a$$,$$b$,求一组$(x,y)$满足$ax+by=c$定理无解:$c\mid\gcd(a,b)$不成立有解a,b中一个为负则对其加另一个直至其为正,两个为负则翻转正负(包括答案) voidex_gcd(inta,intb,int&x,int&y){if(!b)x=1,y=0;elseex_gcd(b,a%b,......
  • 采样定理
    信号\(x(t)\)的频谱为\(X(\mathrm{j}\omega)\)。对信号使用周期单位冲激串采样得到采样信号\(x_p(t)\):\[x_p(t)=x(t)p(t)\]其中,\(p(t)\)为采样函数,是周期为\(T\)的周期单位冲激串,并且\(p(t)\)的基波频率\(\omega_s=\frac{2\pi}{T}\)。采样函数如下:\[p(t)=\sum_{......
  • MFC-GetExtendedStyle获取扩展样式
     DWORDExStyles=mylist4.GetExtendedStyle();//获取扩展样式DWORDoldstyle=mylist4.SetExtendedStyle(ExStyles|LVS_EX_FULLROWSELECT);//设置扩展样式/*指定的扩展样式LVS_EX_GRIDLINES//绘制表格LVS_EX_SUBITEMIMAGES//......
  • 中国剩余定理(CRT)学习笔记
    约定\(A\perpB\)表示\(\gcd(A,B)=1\)。\(A\midB\)表示\(B\equiv0\pmod{A}(A\neq0)\)。引入考虑以下这道题:有物不知其數,三三數之剩二,五五數之剩三,七七數之剩二。問物幾何?——《孫子算經》也就是说,求出下列关于\(x\)方程组的最小整数解:\[\begin{cases}x\equi......
  • (Edge,Chrome)编写扩展应用,替代IE ActiveX插件
    资料来源#这次以Edge作为例子,Chrome其实也差不多Edge扩展应用资料:https://docs.microsoft.com/zh-cn/microsoft-edge/extensions-chromium用到的浏览器Api资料:https://developer.mozilla.org/zh-CN/docs/Mozilla/Add-ons/WebExtensions/API/runtime/sendMessagehttps://developer......