扩展中国剩余定理 exCRT 学习笔记
前言
由于 \(\{\mathrm{CRT}\}\subseteq\{\mathrm{exCRT}\}\),而且 CRT 又太抽象了,所以直接学 exCRT 了。
这么抽象的东西我怎么可能会写
前置技能
- gcd/lcm
- exgcd
- 快速乘
参考资料
- 一篇未通过的洛谷日报 by AH_ljq
- 比较直观的 exCRT 学习笔记 by Milky Way
- 我之前写过的 exgcd 学习笔记
- huyufeifei 对 CRT 的劝退
用途
用于求一个关于 \(x\) 的同余方程组
\begin{matrix}
x\equiv a_1\pmod{b_1}\\\
x\equiv a_2\pmod{b_2}\\\
\cdots\\\
x\equiv a_n\pmod{b_n}
\end{matrix}\right.
\]
的解。
解法推导
考虑到如果 \(x\) 对原方程组成立,那么 \(x\) 对其中任意几个方程也都成立。那么如果要满足 \(x\) 对前 \(i\) 个方程都成立,一个必要条件是对前 \(i-1\) 个方程都成立。
采用一种类似“增量法”的思路来合并每个式子。
假定前 \(i-1\) 个式子已经合并完了,得到 \(x\equiv a_{i-1}\pmod{b_{i-1}}\)。此时合并
\begin{aligned}
&x\equiv a_{i-1}&\pmod{b_{i-1}}\\\
&x\equiv a_i&\pmod{b_i}
\end{aligned}\right.
\]
它的本质是存在 \(k_1,k_2\),满足(上下颠倒了一下,方便下面推导)
\begin{aligned}
&x+k_1b_i=a_i\\\
&x+k_2b_{i-1}=a_{i-1}
\end{aligned}\right.
\]
两式相减,得到
\]
此时,仅有 \(k_1,k_2\) 是变量,剩下的都已知。但是我们用 exgcd 解方程时,要求等式右边是 \(\gcd(b_i,b_{i-1})\)。此时我们列出一个新方程
\]
解出 \(k_1'\) 的一个特解 \(k_0\)。
那么为了满足新方程和原方程都成立,那么
\]
因此 \(k_1\) 有特解
\]
则 \(k_1\) 的通解是
\]
不定方程的通解
对于 \(ax+by=c\) 这个方程,假定我们已经解出来一组解为
\[\left\{
\begin{matrix}
x=x_0\\\
y=y_0
\end{matrix}
\right.
\]我们带入后对左侧式子变形,得到 \(ax_0+S+by_0-S=c\)。
我们要从 \(ax_0+S\) 中提一个 \(a\) 出来,从 \(by_0-S\) 中提一个 \(b\) 出来,得到 \(a(x_0+\frac Sa)+b(y_0-\frac Sb)=c\)。
此时最小的满足 \(x_0+\frac Sa,\ y_0-\frac Sb\) 都是整数的 \(S\) 就是 \(a,b\) 的最小公倍数了。
那么每个相邻的 \(x\) 之间相差 \(\frac Sa=\frac{\operatorname{lcm}(a,b)}a=\frac b{\gcd(a,b)}\),即 \(x\) 的通解为 \(x_0+t\times \frac b{\gcd(a,b)},\ t\in \mathbb Z\)。
后面的式子仅与 \(a,b\) 有关。
又因为
\]
我们带入 \(k_1\) 的通解,得
\]
化简
t\times\operatorname{lcm}(b_i,b_{i-1})=a_i-x-b_ik_{1_0},\ t\in \mathbb Z
\]
由于 \(t\in \mathbb Z\),所以
a_i-x-b_ik_{1_0}\equiv 0\pmod{\operatorname{lcm}(b_i,b_{i-1})}\\\
x\equiv a_i-b_ik_{1_0}\pmod{\operatorname{lcm}(b_i,b_{i-1})}
\]
\(\equiv\) 符号右边的都是已知量,这时我们就把两个方程结合到一块了。
我们依次合并 \((2)=[(1),(2)],(3)=[(2),(3)],\cdots,(n)=[(n-1),(n)]\),最终得到第 \(n\) 个式子
\]
(\(a_n\) 是合并 \(n-1\) 个式子后的新 \(a_n\),上述过程中每一步都会更新 \(a_i\))
此时 \(a_n\) 就是解了,取模取正数就是最小正整数解。
(上下两部分实质相同)
解法简述
对于
\begin{aligned}
&x\equiv a_{i-1}&\pmod{b_{i-1}}\\\
&x\equiv a_i&\pmod{b_i}
\end{aligned}\right.
\]
有
\begin{aligned}
&x+k_1b_i=a_i\\\
&x+k_2b_{i-1}=a_{i-1}
\end{aligned}\right.
\]
解出
\]
的特解 \(k_1'\),两边同乘 \(\frac{a_i-a_{i-1}}{\gcd(b_i,b_{i-1})}\),得到
\]
那么 \(k_ 1'\) 也被乘了 \(\frac{a_i-a_{i-1}}{\gcd(b_i,b_{i-1})}\),因此 \(k_1\) 有特解
\]
而 \(k_1\) 的通解是
\]
代回去就是
\]
得到
\]
作为第 \(i\) 个方程即可。
注意事项
在 \(k_1,k_{1_0},a_i\) 的计算中是有可能爆 long long
的,而这三个计算恰好又都是模意义下,所以使用快速乘。
一般情况下不考虑 \(\operatorname{lcm}(b_1,b_2,\cdots,b_n)\) 爆 long long
的情况。
代码
#include<cstdio>
#define ll long long
ll exgcd(ll a,ll b,ll &x,ll &y)
{
if(!b)
{
x=1,y=0;
return a;
}
ll g=exgcd(b,a%b,y,x);
y-=a/b*x;
return g;
}
ll qmul(ll x,ll y,ll p)
{
ll ans=0,f=1;
if(x<0)
{
x=-x;
f=-f;
}
if(y<0)
{
y=-y;
f=-f;
}
while(y)
{
if(y&1)
ans=(ans+x)%p;
x=(x+x)%p;
y>>=1;
}
return ans*f;
}
ll a[100100],b[100100];
int main()
{
int n;
ll x,y;
scanf("%d",&n);
for(int i=1;i<=n;++i)
scanf("%lld%lld",&b[i],&a[i]);
for(int i=2;i<=n;++i)
{
ll g=exgcd(b[i],b[i-1],x,y);
ll t=b[i-1]/g;//k1 的最小波动 Δ
x=qmul(x,(a[i]-a[i-1])/g,t);
x=(x%t+t)%t;
a[i]-=qmul(x,b[i],b[i]/g*b[i-1]);
b[i]=b[i]/g*b[i-1];
a[i]=(a[i]%b[i]+b[i])%b[i];
}
printf("%lld\n",(a[n]%b[n]+b[n])%b[n]);
return 0;
}
扩展中国剩余定理 exCRT 学习笔记的更多相关文章
- 扩展中国剩余定理 (ExCRT)
扩展中国剩余定理 (ExCRT) 学习笔记 预姿势: 扩展中国剩余定理和中国剩余定理半毛钱关系都没有 问题: 求解线性同余方程组: \[ f(n)=\begin{cases} x\equiv a_1\ ...
- 扩展中国剩余定理 (exCRT) 的证明与练习
原文链接https://www.cnblogs.com/zhouzhendong/p/exCRT.html 扩展中国剩余定理 (exCRT) 的证明与练习 问题模型 给定同余方程组 $$\begin{ ...
- 中国剩余定理(CRT) & 扩展中国剩余定理(ExCRT)总结
中国剩余定理(CRT) & 扩展中国剩余定理(ExCRT)总结 标签:数学方法--数论 阅读体验:https://zybuluo.com/Junlier/note/1300035 前置浅讲 前 ...
- 扩展中国剩余定理(EXCRT)快速入门
问题 传送门 看到这个问题感觉很难??? 不用怕,往下看就好啦 假如你不会CRT也没关系 EXCRT大致思路 先考虑将方程组两两联立解开,如先解第一个与第二个,再用第一个与第二个的通解来解第三个... ...
- CRT(中国剩余定理)学习笔记
先扔个模板题.链接. 简化题意:他让我求 \(x \equiv a_i \pmod{m_i}\) 的解. 例如,\( \begin{cases} x \equiv 1 \pmod{3} \\ x \e ...
- 扩展中国剩余定理(EXCRT)学习笔记
扩展中国剩余定理(EXCRT)学习笔记 用途 求解同余方程组 \(\begin{cases}x\equiv c_{1}\left( mod\ m_{1}\right) \\ x\equiv c_{2} ...
- P4777 【模板】扩展中国剩余定理(EXCRT)/ poj2891 Strange Way to Express Integers
P4777 [模板]扩展中国剩余定理(EXCRT) excrt模板 我们知道,crt无法处理模数不两两互质的情况 然鹅excrt可以 设当前解到第 i 个方程 设$M=\prod_{j=1}^{i-1 ...
- P4777 【模板】扩展中国剩余定理(EXCRT)
思路 中国剩余定理解决的是这样的问题 求x满足 \[ \begin{matrix}x \equiv a_1(mod\ m_1)\\x\equiv a_2(mod\ m_2)\\ \dots\\x\eq ...
- P4777 【模板】扩展中国剩余定理(EXCRT)&& EXCRT
EXCRT 不保证模数互质 \[\begin{cases} x \equiv b_1\ ({\rm mod}\ a_1) \\ x\equiv b_2\ ({\rm mod}\ a_2) \\ ... ...
随机推荐
- [JAVA] 冻结Excel的第一行或第一列
可以按照如下设置创建冻结窗口. sheet.createFreezePane( 3, 2, 3, 2 ); 前两个参数是你要用来拆分的列数和行数.后两个参数是下面窗口的可见象限,其中第三个参数是右边区 ...
- linux c MQTT客户端实现
linux c MQTT客户端实现 摘自:https://www.jianshu.com/p/d309de966379 一.前言:mqtt协议是轻量级的消息订阅和发布(publish/subscrib ...
- python 的输入和输出
内置函数:raw_inpurt('place input') print getpass 隐藏输入密码 import getpass pwd = getpass.getpass("> ...
- Transferring Data Between ASP.NET Web Pages
14 July 2012 20:24 http://www.mikesdotnetting.com/article/192/transferring-data-between-asp-net-web- ...
- CodeForces 339D Xenia and Bit Operations (线段树)
题意:给定 2的 n 次方个数,对这些数两个两个的进行或运算,然后会减少一半的数,然后再进行异或运算,又少了一半,然后再进行或运算,再进行异或,不断重复,到最后只剩下一个数,要输出这个数,然后有 m ...
- OpenCV源码解析
OpenCV K-means源码解析 OpenCV 图片读取源码解析 OpenCV 视频播放源码解析 OpenCV 追踪算法源码解析 OpenCV SIFT算法源码解析 OpenCV 滤波源码分析:b ...
- vmware中安装centos 6.7
centos 6.7 软件下载地址:http://b.mirrors.lanunion.org/CentOS/6.7/isos/i386/ 引用:http://www.cnblogs.com/sees ...
- 不用SQL给打印记录编号
以QUICKREPORT为例 页面设置如下: 其中ID为编号. 设置为表的ID字段. QUICKREPORT所在的FORM添加一个变量: var FprnT6: TFprnT6; Vxh:intege ...
- CSS基础知识:常见选择器示例
CSS(Cascading Style Sheet),中文译为层叠样式表,可以让设计者方便灵活地控制Web页面的外观表现.CSS是1996年由W3C审核通过并且推荐使用的.CSS的引入,就是为了使HT ...
- C++裁剪文件,截断文件,_chsize()
errno_t _chsize_s( int fd, __int64 size ); 详见msdn知识库 _chsize将文件裁剪为指定大小,大小的度量方法与 long ftell(FILE * fp ...