【数论】扩展欧几里得

大家一定都知道欧几里得算法gcd(a,b)，它能够在O(lg(ab))的时间内求出整数a和b的最大公约数，是一个非常优秀的算法。这么一个优秀的算法，能不能在不增加时间复杂度的情况下，让它能处理更多的信息，以便解决更多的问题呢？今天，我们来介绍一下 扩展欧几里得算法。

算法解决对象

扩展欧几里得算法可以用来求形如ax+by=GCD(a,b)的方程的 整数解。（当然，a，b是整数）

算法证明

首先，一定存在解符合题意。（我也不知道为什么？）

其次，因为GCD(a,b)=GCD(b,a%b)，所以x*a+y*b=GCD(a,b)=GCD(b,a%b)=p*b+q*a%b=p*b+q*(a-a/b*b)=q*a+(p-a/b*q)*b. 说成人话就是： 假设方程 b*x+a%b*y=GCD(b,a%b) 的解为x0,y0;那么原方程的解就是 x=y0,y=x0-(a/b)*y0.

基本情况（Base Case）：当b=0时，x=1,y=0.

算法实现

#include<bits/stdc++.h>
using namespace std;
int extended_gcd(int a, int b, int &x, int &y){
 int res,x0,y0;
 if (!b){
  x=1;
  y=0;
  return a;
 }
 res=extended_gcd(b,a%b,x0,y0);
 x=y0; y=x0-a/b*y0;
 return res;
} 
int main(){
 int a,b,x,y,z;
 scanf("%d%d", &a, &b);
 z=extended_gcd(a,b,x,y);
 printf("%d %d %d\n",z,x,y);
 return 0;
} 
/*
模板名称：扩展欧几里得
模板作用：求方程ax+by=GCD(a,b)的一组整数解
*/