数论---扩展欧几里得算法，线性同余，求逆元

裴蜀定理

内容：对于任意的整数a，b，一定存在两个未知数 x，y ，使得 ax+by =d (d=(a,b))

对于任意的x，y，所得到的 ax+by = d的倍数

对于任意的x，y，所得到的最小正整数是d

推论：如果两个x，y互质，ax+by 最小正整数 =1 ，任意的x，y可以表达所有的数

扩展欧几里得算法

内容：

已知整数a、b，扩展欧几里得算法可以在求得a、b的最大公约数的同时，能找到整数x、y（其中一个很可能是负数），使它们满足裴蜀定理

证明：

#include <iostream>
#include <algorithm>
using namespace std;

int exgcd(int a,int b,int &x,int &y)
{
    if(!b)
    {
        x=1,y=0;
        return a;
    }
    int x1,y1,gcd;
    gcd=exgcd(b,a%b,x1,y1);
    x=y1;
    y=x1-a/b*y1;
    return gcd;
}
int main()
{
    int n;
    cin>>n;
    while(n--)
    {
        int a,b,x,y;
        cin>>a>>b;
        exgcd(a,b,x,y);
        cout<<x<<' '<<y<<endl;
    }
    return 0;
}

线性同余

内容： 给定a，b，m ， 求出最大公约数d和扩展欧几里得算法的x’

当a*x 余数= b (mod m) 时，当有解时，x = (b/d)*x'%m

证明：

 

#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
int exgcd(int a,int b,int &x,int &y)
{
    if(!b)
    {
        x=1,y=0;
        return a;
    }
    int d=exgcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
int main()
{
    int n;
    cin>>n;
    while(n--)
    {
        int a,b,m;
        cin>>a>>b>>m;
        int x,y,d;
        d=exgcd(a,m,x,y);
        if(b%d) cout<<"impossible\n";
        else cout<<1ll*b/d*x%m<<endl;
    }
    return 0;
}

求逆元

求逆元就让上面的b=1，就是求a的逆元