拓展欧几里得算法

因为博主是个sb，所以下文内容有可能不对，推荐大家去看这篇博文，对不起了
欧几里德算法部分我们好像只能用来求解最大公约数，但是扩展欧几里德算法就不同了，我们既可以求出最大公约数，还可以顺带求解出使得： a*x + b*y = gcd 的通解 x 和 y
扩展欧几里德定律：　
对于不完全为0的非负整数a，b，gcd(a, b)表示a, b的最大公约数，必定存在整数对x，y，满足a*x+b*y==gcd(a, b)。
证明：（转）
a*x1+b*y1=gcd(a, b);

b*x2+(a%b)*y2=gcd(b, a%b);

因为由欧几里德定理知：gcd(a, b)==gcd(b, a%b)

所以a*x1+b*y1=b*x2+(a%b)*y2; 　　因为r=a%b, r =a-k*b所以==>

a*x1+b*y1=b*x2+(a-k*b)*y2; 　　因为k=a/b;所以　==>

a*x1+b*y1=b*x2+(a-(a/b)*b)*y2; 　　展开得到　　==>　　　　

a*x1+b*y1=b*x2+a*y2-b*(a/b)*y2;　　转换得到 ==>

a*x1+b*y1=a*y2+b*(x2-(a/b)*y2);

观察上式可知 x1=y2, y1=x2-a/b*y2;

由此可知x1，y1是由x2，y2得出来的，由此类推x2，y2是由x3，y3得出来的，

那什么时候是终止呢？也就是递归gcd(a, b)中b=0时；也就是说此时a的值就是要求得最大公约数

即gcd(a, 0)此时由扩展欧几里得定律a*x+b*y==gcd(a, b)

知 a*x+b*y=a;

解出x=1, y=0;

此时就是递归终止的地方：
代码如下

#include <cstdio>
using namespace std;
int a,b,x,y;
int e_gcd(int a,int b,int &x,int &y)
{
    if(!b)
    {
        x=1;
        y=0;
        return a;
    }

    int ans=e_gcd(b,a%b,x,y);
    int tmp=x;

    x=y;
    y=tmp-a/b*y;

    return ans;
} 
int main()
{


    scanf("%d%d",&a,&b);

    e_gcd(a,b,x,y);

    printf("%d %d",x,y);

    return 0;
} 

那么问题来了，这破东西有啥用呢
问的好，它可以用来求一个同余方程的解，也就是逆元
ax ≡ 1 (mod b)，现在找x能够使这个式子成立
这个式子等价于ax+by=1
这不就是拓展欧几里得的表达式吗，从这里我们可以看出如果gcd(a,b)!=1那么这个方程就不会有解
所以说呢
形如a*x + b*y = c这样的式子，若c%gcd(a,b)==0，那么这个方程就会有解
但是一般题目里会让你求一个最小的x，当你用拓欧求出一个解时，一般会让你去找一个最小解，我们只需要对这个数取模b就行了（如果求正数，你只需要先加一个b，再取模行了，应该都知道吧）
代码如下（我做的板子题爆int了，不要介意开了一串long long）

#include <cstdio>
using namespace std;
long long a,b,x,y;
long long e_gcd(long long a,long long b,long long &x,long long &y)
{
    if(!b)
    {
        x=1;
        y=0;
        return a;
    }

    long long ans=e_gcd(b,a%b,x,y);
    long long tmp=x;

    x=y;
    y=tmp-(a/b)*y;

    return ans;
} 
int main()
{
    scanf("%lld%lld",&a,&b);

    e_gcd(a,b,x,y);

    printf("%lld",(x+b)%b);

    return 0;
}