题解 P1082 【同余方程】
这里给出非递归的 exgcd 做法
【基础】
( 只需要非递归的同学麻烦跳过 )
由于欧几里德算法 ( 又名辗转相除法 ) 可以帮助我们求出最大公约数,并且提出对于
\(\forall a,b\in Z_+,gcd(a,b)|c\) 则 $ ax+by=c $ 一定有整数解
因此,在 \(gcd\) 的基础上,我们提出可以求解 \(x,y\) 的算法:拓展欧几里德算法(exgcd)
如果我们要求解 \(ax+by=gcd(a,b)\)
我们可以知道 \(gcd(a,b)=gcd(b,a\%b)\)
而还有 \(gcd(b,a\%b)=bx'+(a\%b)y'\)
所以,很快可以得到:\(ax+by=bx'+(a\%b)y'\)
\(\because\) 很显然 \(a\%b=a-\lfloor{a\over b}\rfloor\times b\)
\(\therefore\) 我们把上式打开可得到:
\(\qquad ax+by=bx'+ay'-\lfloor{a\over b}\rfloor\times by'\)
\(\qquad ax+by=ay'+b(x'-\lfloor{a\over b}\rfloor y')\)
因此,我们要求得 \(x,y\) 就要先求出 \(x',y'\)
而对于 \(ax+0y=gcd(a,0)=a\) ,显然有解为 \(x=1,y=0\)
这就是递归边界
因此,我们可以写出求 \(x,y\) 的程序:
void exgcd(int a,int b,int &x,int &y){
if(!b) { x=1; y=0; return ; }
exgcd(b,a%b,x,y);
x-=a/b*y;
swap(x,y);
}
那这根这一题有什么关系呢?
你想想,我们要求 \(ax\equiv1(\mod b)\)
是不是相当于求 \(ax=1+by\)
即 \(ax+b(-y)=1\) ?
看到了吗,exgcd 出现了!
所以就是再加一句
inline int ny(int a,int b){
int x,y;
exgcd(a,b,x,y);
x=(x%b+b)%b;
return x;
}
【分析】
这边讲一下如何非递归的实现 exgcd :
有一部分原理麻烦看一下 这篇文章
已知我们只是需要每次出现的 \(a/b\) 这个值,而且要后面出现的先算,所以我们就把它们压到一个栈里面
算的时候再弹出来就可以了:
int a=read(),b=read(),m=b,tmp[10000]={0},cur=1,x=1,y=0;
tmp[0]=a/b;
while(b^=a^=b^=a%=b) tmp[cur++]=a/b;
while(cur--) y^=x^=y^=x-=tmp[cur]*y;
x=(x%m+m)%m;
【代码】
核心代码我已经放在上面了,现在放出来的是本蒟蒻 码风极丑的 代码:
#include<cstdio>
using namespace std;
inline int read(){
register int ans=0;register char c=getchar(); while(c<48||c>57) c=getchar();
while(c>=48&&c<=57) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();
return ans;
}
int main(){
register int a=read(),b=read(),m=b,tmp[10000]={0},cur=1,x=1,y=0; tmp[0]=a/b;
while(b^=a^=b^=a%=b) tmp[cur++]=a/b;
while(cur--) y^=x^=y^=x-=tmp[cur]*y; x%=m;
printf("%d",(x<0)?(x+m):x);
return 0;
}
最后安利一下 本蒟蒻的博客
题解 P1082 【同余方程】的更多相关文章
- 洛谷——P1082 同余方程
P1082 同余方程 题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输 ...
- 洛谷P1082 同余方程 [2012NOIP提高组D2T1] [2017年6月计划 数论06]
P1082 同余方程 题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输 ...
- 洛谷P1082 同余方程 题解
题目链接:https://www.luogu.com.cn/problem/P1082 题目大意: 求关于 \(x\) 的同余方程 ax≡1(mod b) 的最小正整数解. 告诉你 \(a,b\) 求 ...
- Luogu P1082 同余方程(NOIP 2012) 题解报告
题目传送门 [题目大意] 求关于x的同余方程 ax≡1(mod b)的最小整数解. [思路分析] 由同余方程的有关知识可得,ax≡1(mod b)可以化为ax+by=1,此方程有解当且仅当gcd(a, ...
- 【luogu P1082 同余方程】 题解
最近一直在学习数论,讲得很快,害怕落实的不好,所以做一道luogu的同余方程练练手. 关于x的同余方程 ax ≡ 1 mod m 那么x其实就是求a关于m的乘法逆元 ax + my = 1 对于这个不 ...
- 洛谷 P1082 同余方程 题解
每日一题 day31 打卡 Analysis 题目问的是满足 ax mod b = 1 的最小正整数 x.(a,b是正整数) 但是不能暴力枚举 x,会超时. 把问题转化一下.观察 ax mod b = ...
- [NOIP2012] 提高组 洛谷P1082 同余方程
题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输出只有一行,包含一个正 ...
- 洛谷P1082 同余方程
题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输出只有一行,包含一个正 ...
- [Luogu P1082]同余方程
题目链接 这道题求关于x的同余方程ax≡1(mod b)的最小正整数解.换而言之方程可以转换为ax+by=1,此时有y为负数.此时当且仅当gcd(a,b)|1时,方程有整数解. 于是乎这道题就变成了a ...
随机推荐
- 045、Java中使用if语句进行判断
01.代码如下: package TIANPAN; /** * 此处为文档注释 * * @author 田攀 微信382477247 */ public class TestDemo { public ...
- JAVA字符串比较问题
在java中值类型通过==来进行比较值是否相等 而字符串作为一种引用类型,通过==是用来比较其内存位置的,使用equals才是用来比较其值是否相等 使用equals时养成将字符串放在前面的好习惯 字符 ...
- 20200119日志 EPLAN高压房 VFD 单线图 心得
提纲: EPLAN 画单线图的方法,可以先绘制原理图,然后在一个柜子里面的器件 用方框圈起来.方框名称一样.注意 一个工程里面的器件编号是唯一的. 断路器选型. PT 手车 接地刀作用 具体内 ...
- java 三羊献瑞
三羊献瑞 观察下面的加法算式: 其中,相同的汉字代表相同的数字,不同的汉字代表不同的数字. 请你填写"三羊献瑞"所代表的4位数字(答案唯一),不要填写任何多余内容. public ...
- python 4个人中有一人做了好事,一直有三个人说了真话,根据下面的对话判断是谁做的好事
题目:4个人中有一人做了好事,一直有三个人说了真话,根据下面的对话判断是谁做的好事. A:不是我 B:是C C:是D D:C胡说 for inum in ['A','B','C','D']: prin ...
- 《Python爬虫技术:深入理解原理、技术与开发》已经出版,送Python基础视频课程
好消息,<Python爬虫技术:深入理解原理.技术与开发>已经出版!!! JetBrains官方推荐图书!JetBrains官大中华区市场部经理赵磊作序!送Python基础视频课程!J ...
- 10.swoole学习笔记--进程队列通信
<?php //进程仓库 $workers=[]; //最大进程数 $worker_num=; //批量创建进程 ;$i<$worker_num;$i++){ //创建子进程 $proce ...
- 003、mysql输出多个结果
SELECT VERSION(); SELECT NOW(); 结果1: 结果2: 不忘初心,如果您认为这篇文章有价值,认同作者的付出,可以微信二维码打赏任意金额给作者(微信号:382477247)哦 ...
- 一个swift版简单的用户名和密码输入textField
http://www.code4app.com/thread-31992-1-1.html 常见的动画提交按钮 http://www.code4app.com/thread-32239-1-1.htm ...
- [APIO2012]派遣 可并堆
Background 在一个忍者的帮派里,一些忍者们被选中派遣给顾客,然后依据自己的工作获取报偿. Description 在这个帮派里,有一名忍者被称之为Master.除了Master以外,每名忍者 ...