扩展欧几里得算法
思路分析:
根据裴蜀定理知:
gcd(a, b): a 和 b 的最大公约数
对于任意一对正整数a, b, 那么一定存在整数x, y, 使得ax + by = gcd(a, b)
ax + by = d, d 一定是 gcd(a, b)的倍数,且其最小正整数就是gcd(a,b)
证明:
因为a 是 gcd(a,b)的倍数,b 也是gcd(a,b)的倍数,所以ax + by 也一定是gcd(a,b)倍数,正整数倍数最小为1倍即gcd(a,b)
扩展欧几里得算法(exgcd):
用于求解ax + by = gcd(a,b)中x, y 的解。
根据b = 0 与否进行分类讨论:
1)b == 0
ax + by = gcd(a, b)
ax + 0 *y = gcd(a, 0)
ax = a
故而求得解:x = 1, y = 0
2 )b != 0
根据辗转相除法知:gcd(a, b) = gcd(b, a %b)
题目:
给定 n 对正整数 ai,bi,对于每对数,求出一组 xi,yi,使其满足 ai×xi+bi×yi=gcd(ai,bi)。
输入格式
第一行包含整数 n。
接下来 n 行,每行包含两个整数 ai,bi。
输出格式
输出共 n 行,对于每组 ai,bi,求出一组满足条件的 xi,yi,每组结果占一行。
本题答案不唯一,输出任意满足条件的 xi,yi 均可。
数据范围
1≤n≤105,
1≤ai,bi≤2×109
输入样例:
2
4 6
8 18
输出样例:
-1 1
-2 1
实现代码:
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
int n,a,b,x,y;
void exgcd(int a,int b,int &x,int &y){
if(b==0){
x=1;y=0;
return;
}
int x1,y1;
exgcd(b,a%b,x1,y1);
x=y1;y=x1-a/b*y1;
}
int main()
{
cin>>n;
while(n--){
cin>>a>>b;
exgcd(a,b,x,y);
cout<<x<<" "<<y<<endl;
}
return 0;
}
扩展:
1、对于一般方程的求解:ax + by = c
这个是不是和我们的扩展欧几里得算法很像,我们只需要将欧几里得算出来的答案x和y的解,都乘以c/gcd(a,b)即可。
证明:
求解ax’ + by‘ = c
因为:ax + by = d
(ax + by) * c/d = d * c/d
(ax + by) * c/d = c
a(x * c/d) + b(y * c/d) = c
得x’和y‘的解:x‘ = x * c/d, y‘ = y * c/d
2、求解一次同余方程ax ≡ b (mod p)
ax ≡ b (mod p)
ax = p * (-y) + b
ax + py = b
当gcd(a, p)| b,然后用扩展欧几里得求解即可
特别当b = 1,且a和p互质,所求x为a的逆元。