数论（总）

一.求最大公约数

辗转相除法

int gcd(int a,int b){
    if(a%b==0){
        return b;
    }else{
        return gcd(b,a%b);
    }
}

二.裴蜀定理与扩展欧几里得算法

裴蜀定理：

总存在整数x，y使 ax+by=gcd(a，b) 成立。

证明方法就是辗转相除法逆过来。

裴蜀定理的推论：

1.a，b最大公约数的另一种理解方法是：找到一组整数x，y使 ax+by 得到最小正数解。

2.假如方程 ax+by=c 有整数解，则c一定是gcd(a，b)的倍数。

3.可扩展到多项式。

扩展欧几里得算法：

寻找 ax+by=gcd(a，b) 的一个解。

原理（公式推导）：

设g为最大公约数，gcd调用时，第k次输入的参数为a[k],b[k];

在最后一次调用中，g=x*a[n]+y*b[n]。

(实际上最后一次x==1，y==0)

则可推导上一次调用时的x，y：g=x*b[n-1]+y*((a[n-1]-(a[n-1]/b[n-1])*b[n-1])

(a[n]对应b[n-1]，b[n]对应(a[n-1]/b[n-1])*b[n-1]，即a%b)

整理得：

g = y*a[n-1] + (x-y*(a[n-1]/b[n-1]))*b[n-1];

这样一直推到初始的两个a，b，就得出一个解x，y。

代码：

int a,b,x,y,g;

void exgcd(int a,int b){
    if(b==0){
        g=a;
        x=1;
        y=0;
    }else{
        exgcd(b,a%b);
        int px=x,py=y;
        x=py;
        y=px-py*(a/b);
    }
}