洛谷 P1082 同余方程(同余&&exgcd)
嗯...
题目链接:https://www.luogu.org/problem/P1082
这道题很明显涉及到了同余和exgcd的问题,下面推导一下:
首先证明有解情况:
ax + by = m有解的必要条件是 m mod gcd(a, b) = 0
a为gcd(a, b)的倍数,b为gcd(a, b)的倍数,x、y为整数,
所以ax + by是gcd(a, b)的倍数,所以m是gcd(a, b)的倍数
然后证明a、b互质(下面会用到):
本题中1 mod gcd(a, b) = 0,所以gcd(a, b) = 1,所以a、b互质
同余:
a≡b(mod n) --> 含义:a和b关于模n同余,即 a mod n = b mod n。
所以不难推出,a≡b的充要条件是a-b是n的整数倍(a > b)。
因为a-b是n的整数倍,所以a-b = ny(y为倍数)
所以,根据同余,我们可以把本题中的同余式转化为(注:这里的a.b与上文不同):
ax≡1(modb) --> ax % b = 1 % b --> ax - 1 = by --> ax - by = 1
下一步,便进行exgcd(关于exgcd的证明见:https://www.luogu.org/problemnew/solution/P1082),分别求出ax - by = 1中x和y的值。
最后进行答案处理:
因为答案要求是x的最小正整数,所以我们进行一个答案处理:x = (x % b + b) % b
证明其正确性:
设新得到的x为xn,x = kb + q(q < b)则:
x % b = q ,把x % b = q带入 xn = (x % b + b) % b,得
xn = (x % b + b) % b = (q + b) % b = (q % b + b % b) % b = q % b = q
把xn = q带入x = kb + q,得,x = kb + xn, 所以xn = x - kb,然后再根据下面的推导得知它是正确的...
证明:
x批量地减去或加上 b,能保证 ax + by = ax + by = 1:
ax + by = 1
ax + by + k*ba - k*ba = 1
a (x + kb) + (y - ka) b = 1
1.显然这并不会把方程中任何整数变成非整数。
2.“加上或减去 b”不会使 x 错过任何解。可以这么理解:
已经求出一组整数 x,y 使得 ax+by =1 ,也就是(1 - ax) / b = y。y 是整数,可见目前 1−ax 是 b 的倍数。现在想改变 x并使得方程仍然成立。
已知 a,b 互质,假若x的变化量Δx不是b的倍数,则1−ax 的变化量−a*Δx也不是 bb 的倍数,这会使得1-ax不再是b的倍数,则y不是整数了。
仅当x的变化量是b的倍数时,1−ax能保持自己是b的倍数,此时就出现新的解了。
AC代码:
#include<cstdio>
#include<iostream> using namespace std;
//ax % b == 1 % b --> ax - 1 = y * b --> ax - yb == 1 long long d, x, y; inline void exgcd(long long a, long long b, long long &d, long long &x, long long &y){
if(!b) {d = a; x = ; y = ;}
else{ exgcd(b, a % b, d, y, x); y -= x * (a / b);}
}
int main(){
long long a,b;
scanf("%lld%lld", &a, &b);
exgcd(a, b, d, x, y);
x = (x % b + b) % b;
printf("%lld", x);
return ;
}
AC代码
洛谷 P1082 同余方程(同余&&exgcd)的更多相关文章
- 洛谷——P1082 同余方程
P1082 同余方程 题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输 ...
- 洛谷P1082 同余方程 [2012NOIP提高组D2T1] [2017年6月计划 数论06]
P1082 同余方程 题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输 ...
- 洛谷 P1082 同余方程 —— exgcd
题目:https://www.luogu.org/problemnew/show/P1082 用 exgcd 即可. 代码如下: #include<iostream> #include&l ...
- 洛谷 P1082 同余方程(exgcd)
题目传送门 解题思路: 因为推导过程过于复杂,懒得写,所以题解传送门 AC代码: #include<iostream> #include<cstdio> using names ...
- [NOIP2012] 提高组 洛谷P1082 同余方程
题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输出只有一行,包含一个正 ...
- 洛谷P1082 同余方程
题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输出只有一行,包含一个正 ...
- 洛谷 P1082 同余方程
题目描述 求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式 输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输出只有一行,包含一个正 ...
- 洛谷P1082 同余方程 题解
题目链接:https://www.luogu.com.cn/problem/P1082 题目大意: 求关于 \(x\) 的同余方程 ax≡1(mod b) 的最小正整数解. 告诉你 \(a,b\) 求 ...
- 洛谷 P1082 同余方程 题解
每日一题 day31 打卡 Analysis 题目问的是满足 ax mod b = 1 的最小正整数 x.(a,b是正整数) 但是不能暴力枚举 x,会超时. 把问题转化一下.观察 ax mod b = ...
随机推荐
- Pandas 数据分析——超好用的 Groupby 详解
在日常的数据分析中,经常需要将数据根据某个(多个)字段划分为不同的群体(group)进行分析,如电商领域将全国的总销售额根据省份进行划分,分析各省销售额的变化情况,社交领域将用户根据画像(性别.年龄) ...
- rabbitmq 命令行与控制台
命令行和管控台 rabbitmqctl stop_app 关闭应用 rabbitmqctl start_app 打开应用 rabbitmqctl status 节点状态 rabbitmqctl add ...
- google插件跨域含用户请求WebApi解决的方案
问题描述: google插件跨域请求WebApi相关解决方案 1.ajax解决含登录用户信息 $.ajax({ url: url, type: "POST", timeout: 6 ...
- ACM-ICPC实验室20.2.21测试-图论(二)
A.患者的编号 给出一个有向图,要求你输出字典序最小的拓扑排序. 常规拓扑排序是做不了的,正解是反向建图,同时用大根堆的优先队列维护,保证每次优先访问编号大的结点,再反向输出~ #include< ...
- element中的分页
在template <template> <div class="Terminal" v-loading="loading"> < ...
- vue中子组件调用父组件里面的数据和方法 父组件调用子组件的数据和方法
1.子组件直接调用父组件的数据和方法 在父组件father,vue <template> <div> <!-- 父组件里面的数据 --> <p>父组件里 ...
- logstash output到kafka记录与总结( No entry found for connection 2)
简述 本文记录logstash的output配置为kafka的过程.这里是简单的例子,输入为stdin,本文主要目的是为了记录在这次配置过程中遇到的问题和解决的过程及总结. 关于kafka集群的搭建可 ...
- JS高级---原型的简单的语法
原型的简单的语法 构造函数,通过原型添加方法,以下语法,手动修改构造器的指向 实例化对象,并初始化,调用方法 <!DOCTYPE html> <html lang="en& ...
- java8date
Java 8:新的时间和日期API 在Java 8之前,所有关于时间和日期的API都存在各种使用方面的缺陷,因此建议使用新的时间和日期API,分别从旧的时间和日期的API的缺点以及解决方法.Java ...
- springboot @Configuration @bean注解作用
@Configuration注解可以达到在Spring中使用xml配置文件的作用 @Bean就等同于xml配置文件中的<bean> 在spring项目中我们集成第三方的框架如shiro会在 ...