扩展欧几里得(exgcd)与同余详解
exgcd入门以及同余基础
gcd,欧几里得的智慧结晶,信息竞赛的重要算法,数论的...(编不下去了
讲exgcd之前,我们先普及一下同余的性质:
- 若,那么
- 若,,且p1,p2互质,
有了这三个式子,就不用怕在计算时溢出了。
下面我会用与分别表示a与b的最大公约数与最小公倍数。
首先会来学扩欧的同学肯定都会欧几里得算法(即辗转相除法)了吧
而通过观察发现:,先除后乘防溢出。
所以与的代码如下:
inline int gcd(int a,int b)
{return (b==)?a:gcd(b,a%b);}
inline int lcm(int a,int b)
{return a/gcd(a,b)*b;
讲exgcd之前先引入一种方程——不定方程
所谓不定方程,是指未知数的个数多于方程个数,且未知数受到某些限制(如要求是有理数、整数或正整数等等)的方程或方程组。
————百度百科
就是形如的方程,其中a,b,c已知。
1.判断是否有解
如果,那么方程无解。
2.转化
方程可转化为,
其中,,
3.求一组特解
接着就用到了exgcd。
我们知道gcd有一个性质
如果,一直循环下去,b将等于0,那么x将等于c/a,y=0。
inline void exgcd(int a,int b,int c)
{
if(!b)
{x=c/a;y=;return;}
exgcd(b,a%b,c);
x=y;
y=(c-a*x)/b;
return;
}
这就求出了一组特解。
exgcd的模板我也在这摆出来
inline void exgcd(int a,int b)
{
if(!b)
{x=;y=;return;}
exgcd(b,a%b);
k=x;x=y;
y=k-a/b*y;
return;
}
这是求时用的,后面讲同余方程会讲。
4.构造通解
我们假设x1,y1是我们求出的一组特解,那么
同余类问题
1.单个同余方程
求的是关于x的解
转化一下,就成了,就可以直接套exgcd模板。
2.同余方程组
1.有解的充要条件
2.
下式减上式得
再用exgcd求出y1和y2,
3.关于通解
所有的x mod lcm(p1,p2)有唯一解,这样就可以通过特解,求通解了。
4.至于式子更多的同余方程组,就先联立两个,就可以得出新的方程
再联立下一个。
exgcd用处及题目讲解
1.求同余方程的解
例如这道题P1082
这是一道裸的扩欧模板题,变形之后就是求。
套模板即可。
inline void exgcd(int a,int b)
{
if(b==)
{x=;y=;return;}
exgcd(b,a%b);
k=x;x=y;
y=k-a/b*y;
return;
}
int main()
{
int n,m;
read(n),read(m);
exgcd(n,m);
printf("%d",(x+m)%m);
}
还有一道模板P1516
仔细观察,推一下后我们发现,这在就是在求的解。
进而可以推出
合并同类项后
把一些东西移到左边来后
把(x-y),(n-m)各看成一个整体后,问题就成了解一个不定方程。
inline int exgcd(long long a,long long b)
{
if(b==)
{x=;y=;return a;}
ans=exgcd(b,a%b);
k=x;x=y;
y=k-a/b*y;
return ans;
}
int main()
{
long long x1,y1,m,n,l;
read(x1),read(y1),read(m),read(n),read(l);
if(n-m<)swap(x1,y1);
exgcd(std::abs(n-m),l);
if((x1-y1)%ans!=)
printf("Impossible");
else
printf("%lld",((x*((x1-y1)/ans))%(l/ans)+(l/ans))%(l/ans));
}
还有一道也是模板P4777,涉及同余方程组求解,上面已详细的讲了,近期我也会发一篇中国剩余定理的博客
inline long long mul(long long a,long long b,long long mod)
{
long long res=;
while(b>)
{
if(b&) res=(res+a)%mod;
a=(a+a)%mod;
b>>=;
}
return res;
}
long long exgcd(long long a,long long b,long long &x,long long &y)
{
if(!b)
{x=;y=;return a;}
long long gcd=exgcd(b,a%b,x,y);
k1=x;x=y;
y=k1-a/b*y;
return gcd;
}
int main()
{
io::begin();
io::read(n);
for(register int i=;i<=n;i++)
io::read(b1[i]),io::read(a1[i]);
long long x,y,k;
long long m=b1[],ans=a1[];
for(int i=;i<=n;i++)
{
long long a=m,b=b1[i],c=(a1[i]-ans%b+b)%b;
long long gcd=exgcd(a,b,x,y);
long long p=b/gcd;
x=mul(x,c/gcd,p);
ans+=x*m;
m*=p;
ans=(ans%m+m)%m;
}
printf("%lld",(ans%m+m)%m);
}
2.扩欧求逆元
这是一种很重要的算法,至于逆元怎么跟扩欧扯上关系,大家可以点这里乘法逆元及两道模板题详解
这里就不多赘述了,大家可以用扩欧a一下P3811,P2613。
我要讲的讲完了,如果觉得讲的还好,请关注我的blog,谢谢
扩展欧几里得(exgcd)与同余详解的更多相关文章
- 同余问题(一)——扩展欧几里得exgcd
前言 扩展欧几里得算法是一个很好的解决同余问题的算法,非常实用. 欧几里得算法 简介 欧几里得算法,又称辗转相除法. 主要用途 求最大公因数\(gcd\). 公式 \(gcd(a,b)=gcd(b,a ...
- 浅谈扩展欧几里得[exgcd] By cellur925
关于扩展欧几里得从寒假时就很迷,抄题解过了同余方程,但是原理并不理解. 今天终于把坑填上了qwq. 由于本人太菜,不会用markdown,所以这篇总结是手写的(什么).(字丑不要嫌弃嘛) ****** ...
- 扩展欧几里得(exgcd)-求解不定方程/求逆元
贝祖定理:即如果a.b是整数,那么一定存在整数x.y使得ax+by=gcd(a,b).换句话说,如果ax+by=m有解,那么m一定是gcd(a,b)的若干倍.(可以来判断一个这样的式子有没有解)有一个 ...
- 扩展欧几里得 exGCD
Elementary Number Theory - Extended Euclid Algorithm Time Limit : 1 sec, Memory Limit : 65536 KB Jap ...
- 数论--扩展欧几里得exgcd
算法思想 我们想求得一组\(x,y\)使得 \(ax+by = \gcd(a,b)\) 根据 \(\gcd(a,b) = \gcd(b,a\bmod b)\) 如果我们现在有\(x',y'\) 使得 ...
- J - 青蛙的约会(扩展欧几里得)
https://vjudge.net/contest/218366#problem/J 第一步追及公式要写对:y+nk-(x+mk)=pL => (n-m)k+lp=x-y 可以看出扩展欧几里得 ...
- EXGCD 扩展欧几里得
推荐:https://www.zybuluo.com/samzhang/note/541890 扩展欧几里得,就是求出来ax+by=gcd(x,y)的x,y 为什么有解? 根据裴蜀定理,存在u,v使得 ...
- exgcd扩展欧几里得求解的个数
知识储备 扩展欧几里得定理 欧几里得定理 (未掌握的话请移步[扩展欧几里得]) 正题 设存在ax+by=gcd(a,b),求x,y.我们已经知道了用扩欧求解的方法是递归,终止条件是x==1,y==0: ...
- UVA 12169 Disgruntled Judge 枚举+扩展欧几里得
题目大意:有3个整数 x[1], a, b 满足递推式x[i]=(a*x[i-1]+b)mod 10001.由这个递推式计算出了长度为2T的数列,现在要求输入x[1],x[3],......x[2T- ...
随机推荐
- 关掉那些windows上因权限无法关闭的服务
ProcessExplorer 下载地址 :https://docs.microsoft.com/en-us/sysinternals/downloads/process-explorer 以管理员身 ...
- 关于nginx unit服务非正常关闭后,无法重新启动问题的处理
昨天在前领导技术大牛吕哥的帮忙下,python服务管理从nginx+supervisor+uwsgi+python3改为了轻便结构nginx + unit + python3,部署和配置起来顿时轻松起 ...
- mongoDB概述
mongoDB介绍 是一个开源的文档数据库,基于分布式文件存储的数据库.由 C++ 语言编写,并领先的 NoSQL 数据库(非关系型数据库)MongoDB 是由 c++语言编写. MongoDB 是一 ...
- 【spring源码分析】IOC容器初始化(六)
前言:经过前几篇文章的讲解,我们已经得到了BeanDefinition,接下来将分析Bean的加载. 获取Bean的入口:AbstractApplicationContext#getBean publ ...
- kvm虚拟化介绍
一.虚拟化分类 1.虚拟化,是指通过虚拟化技术将一台计算机虚拟为多台逻辑计算机.在一台计算机上同时运行多个逻辑计算机,每个逻辑计算机可运行不同的操作系统,并且应用程序都可以在相互独立的空间内运行而互相 ...
- centos7之zabbix简单检查之端口监控
1.创建模板 模板名字可以随便起 2.重点在监控项和触发器上. 3.比如我们创建一个监控25端口号的监控项 net.tcp.service[tcp,,25] 这个选项不是手动写上去的,是需要点击选择. ...
- [转]html5监听任何App自带返回键javascript事件
1.前言 如今h5新特性.新标签.新规范等有很多,而且正在不断完善中,各大浏览器商对它们的支持,也是相当给力.作为前端程序员,我觉得我们还是有必要积极关注并勇敢地加以实践.接下来我将和各位分享一个特别 ...
- 【MySQL 读书笔记】普通索引和唯一索引应该怎么选择
通常我们在做这个选择的时候,考虑得最多的应该是如果我们需要让 Database MySQL 来帮助我们从数据库层面过滤掉对应字段的重复数据我们会选择唯一索引,如果没有前者的需求,一般都会使用普通索引. ...
- 定时任务调度工作(学习记录 二)timer定时函数的用法
schedule的四种用法: 1.schedule(task,time) 参数: task----所安排的任务 time----执行任务的时间 作用: 在时间等于或超过time的时候执行且仅执行一次t ...
- codeforces675D
Tree Construction CodeForces - 675D During the programming classes Vasya was assigned a difficult pr ...