基础数论——EXGCD
1.前言
\(皆さん、こんにちは。\)今天我们来讲 \(EXGCD\) 。(扩展欧几里得)
既然是扩展嘛,那肯定有不扩展的,也就是 \(GCD\) 。
我们都知道 \(GCD\) 怎么写:
ll GCD(ll X,ll Y){ return Y==0?X:GCD(Y,X%Y);}
然后今天就要讲它的升级版: \(EXGCD\) 。
2.讲解
上面一段就是 \(GCD\) 的写法,它能够快速解决最大公约数问题。
那么 \(EXGCD\) 就是要解决下面的这个问题:
\(a\times q\equiv c \pmod{p}\)
已知整数 \(q\) , \(p\) , \(c\) ,求整数 \(a\) 。
那么我们现在把问题转化一下,显然我们不可能直接对这个问题进行求解。
\(a\times q + b\times p = c\)
模 \(p\) 相当于就是加上一些 \(p\) 或者减去一些,我们用一个未知量 \(b\) 来表示 \(p\) 的数量。但相较于第一个式子,这个式子更加简洁。
所以现在我们就是要求出 \(a\) 和 \(b\) 的值。
由于是模 \(p\) 意义下的,所以我们令 \(a\) , \(b\) 都是非负整数。
我们设 \(d=gcd(q,p)\) ,那么可得以下式子:
\(d \times \dfrac{q}{d} \times a + d \times \dfrac{p}{d} \times b = d \times \dfrac{c}{d}\)
由于 \(gcd(q,p) \mid q\) 且 \(gcd(q,p) \mid p\) ,那么 \(\dfrac{q}{d}\) 和 \(\dfrac{p}{d}\) 就是整数,那么 \(\dfrac{c}{d}\) 也是整数,所以 \(gcd(q,p) \mid c\) 。
那么如果我们求解了以下式子:
\(x\times q + y\times p = gcd(q,p)\)
然后令 \(x\) 乘上 \(\dfrac{c}{gcd(q,p)}\) 就是 \(a\) 的值了。
我们设 \(x_1\) , \(y_1\) 为 \(gcd(q,p)\) 的解。
\(x_1 \times q + y_1 \times p = gcd(q,p)\) ①
然后设 \(x_2\) , \(y_2\) 为 \(gcd(p,q\ mod\ p)\) 的解。
\(x_2 \times p + y_2 \times (q\ mod\ p) = gcd(p,q\ mod\ p)\)
显然 \(gcd(p,q\ mod\ p)\) 和 \(gcd(q,p)\) 是相等的,都是求 \(p\) 和 \(q\) 的最大公约数,所以得:
\(x_2 \times p + y_2 \times (q\ mod\ p) = gcd(q,p)\) ②
接下来把 ③ 式的 “\(q\ mod\ p\)” 展开得:
\(x_2 \times p + y_2 \times q - y_2 \times \left\lfloor\dfrac{q}{p}\right\rfloor \times p = gcd(q,p)\) ③
由 ①、③ 二式得, \(x_1=y_2\) 、 \(y_1=x_2-y_2 \times \left\lfloor\dfrac{q}{p}\right\rfloor\) 。
所以我们发现, \(p\) 和 \(q\) 的解是由 \(q\) 和 \((p\ mod\ q)\) 推过来的,所以我们可以在求 \(q\) 和 \(p\) 的最大公约数递归返回的过程中求解 \(x\) 和 \(y\) 。
当 \(p\) 为 \(0\) 的时候 \(q\) 为最大公约数,此时的 \(x=1\) , \(y=0\) 。
这样 \(EXGCD\) 就讲完了。
3.代码
#include<bits/stdc++.h>
using namespace std;
long long Q,P,X,Y;
long long Exgcd(long long Q,long long P,long long &X1,long long &Y1){
if(!P){ X1=1;Y1=0; return Q;}
long long X2,Y2,D=Exgcd(P,Q%P,X2,Y2);
X1=Y2;Y1=X2-(Q/P)*Y2;
return D;
}
int main(){
scanf("%lld%lld",&Q,&P);
Exgcd(Q,P,X,Y);
while(X<0) X=X+P;
printf("%lld\n",X%P);
return 0;
}
\(Thanks for Watch.\)
\(皆さん、さよなら。\)
基础数论——EXGCD的更多相关文章
- LightOJ1214 Large Division 基础数论+同余定理
Given two integers, a and b, you should check whether a is divisible by b or not. We know that an in ...
- HDU-1576 A/B 基础数论+解题报告
HDU-1576 A/B 基础数论+解题报告 题意 求(A/B)%9973,但由于A很大,我们只给出n(n=A%9973) (我们给定的A必能被B整除,且gcd(B,9973) = 1). 输入 数据 ...
- RSA算法原理——(2)RSA简介及基础数论知识
上期为大家介绍了目前常见加密算法,相信阅读过的同学们对目前的加密算法也算是有了一个大概的了解.如果你对这些解密算法概念及特点还不是很清晰的话,昌昌非常推荐大家可以看看HTTPS的加密通信原理,因为HT ...
- POJ 2891 Strange Way to Express Integers 中国剩余定理 数论 exgcd
http://poj.org/problem?id=2891 题意就是孙子算经里那个定理的基础描述不过换了数字和约束条件的个数…… https://blog.csdn.net/HownoneHe/ar ...
- ACM&OI 基础数论算法专题
ACM&OI 基础数学算法专题 一.数论基础 质数及其判法 (已完结) 质数的两种筛法 (已完结) 算数基本定理与质因数分解 (已完结) 约数与整除 (已完结) 整除分块 (已完结) 最大公约 ...
- HDU 1299 基础数论 分解
给一个数n问有多少种x,y的组合使$\frac{1}{x}+\frac{1}{y}=\frac{1}{n},x<=y$满足,设y = k + n,代入得到$x = \frac{n^2}{k} + ...
- NOIp 基础数论知识点总结
推荐阅读 NOIp 数学知识点总结: https://www.cnblogs.com/greyqz/p/maths.html Basic 常用素数表:https://www.cnblogs.com/g ...
- HDU 1060 Leftmost Digit 基础数论
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1060 这道题运用的是数学方法. 假设S=n^n.两边同时取对数,得到lgS=nlgn.即有S=10 ...
- hdoj1905 Pseudoprime numbers (基础数论)
Problem Description Fermat's theorem states that for any prime number p and for any integer a > 1 ...
随机推荐
- 记录一下关于在工具类中更新UI使用RunOnUiThread犯的极其愚蠢的错误
由于Android中不能在子线程中更新ui,所以平时在子线程中需要更新ui时可以使用Android提供的RunOnUiThread接口,但是最近在写联网工具类的时候,有时候会出现联网异常,这个时候为了 ...
- LeetCode--Sort Array By Parity && N-Repeated Element in Size 2N Array (Easy)
905. Sort Array By Parity (Easy)# Given an array A of non-negative integers, return an array consist ...
- iOS的内存管理和引用计数规则、Block的用法以及三种形式(stack、malloc、global)
学习内容 iOS的内存管理和引用计数规则 内存管理的思考方式 自己生成的对象自己持有 非自己生成的对象自己也能持有 自己持有的对象不需要时释放 非自己持有的对象不能释放 ARC有效时,id类型和对象类 ...
- 用Python快速实现一个垃圾分类APP|附带微信小程序
最近北京开始实行垃圾分类,导致大家对垃圾的研究热度突然涨高,垃圾们也纷纷表示从来没有获得过这么高的关注度.其实,上海市去年已经开始实行,网上已经有不少成熟的教程了,像什么<垃圾分类从入门到精通& ...
- 【matlab 基础篇 01】快速开始第一个程序(详细图文+文末资源)
快速入门matlab,系统地整理一遍,如何你和我一样是一个新手,那么此文很适合你: 文章目录 1 软件安装 2 打开软件 3 编写程序 3.1 基础步骤 3.2 添加PATH 3.3 命令行模式 4 ...
- [hdoj5192] 树状数组
枚举所有的区间.对于确定的区间,假设最终的高度为h, 代价是max(∑(Hi−h),∑(h−Hj))(Hi>h,Hj≤h) 等价于max(∑Hi−cnt(i)∗h,cnt(j)∗h−∑Hj) ( ...
- mybatis-generator生成的mapper中的
int updateByExampleSelective(@Param("record") Shop record, @Param("example") Sho ...
- xshell密钥登录服务器
其实很简单 1 xshell 生成pub key . 在工具 -> 用户密钥管理. 生成 .另存为id_rsa_1024.pub 2.服务器上ssh-keygen 3.将生成的文件id_rsa_ ...
- 轻松扩展机器学习能力:如何在Rancher上安装Kubeflow
随着机器学习领域不断发展,对于处理机器学习的团队来说,在1台机器上训练1个模型已经有些难以为继,并且现在业界的共识是机器学习已经不仅仅是简单的模型训练. 在模型训练之前.过程中和之后,需要进行许多活动 ...
- jupyter notebook 修改前端样式
目录 jupyter notebook主题 修改css和js 最终效果 jupyter notebook主题 作者的GitHub地址:https://github.com/dunovank/jupyt ...