HDU 2669 扩展欧几里德求二元不定方程，超水..

看完数论数的这一段，觉得这个算法太精辟了，其实扩展欧几里德算法就是书上的公式翻译成代码，实在是太棒了，我还SB似的去纠结过程，其实自己已经会笔算了

可是书上的第二种解二元一次不定方程方法，目前还在纠结中...

运用的定理中的推论得出一定存在整数 x , y 使得 ax + by =gcd( a,b); 成立. 同样有定理得到 如果 ax + by = n 有解，则 gcd( a,b ) | n 一定成立

令 a1 = a | gcd( a,b ); b1 = b | gcd( a,b ) ; n1 = n | gcd( a,b ) ; ( a1 , b1 ) =1 ; 所以 ax +by = n; 和 a1x + b1y = n1 ; 的解是一样的， 所以要求 a1 x + b1y = n1 的通解

则必须求一特解，则就解方程 a1x + b1y=1; 解的解（x' , y' ) ,扩展解（ x'+t*b1 , y'-t*a1 ）t = 1,2,3,4,.... 利用数论书上第一章的定理来解方程..

则翻译成代码就成了 所谓的扩展欧几里得算法...

 

注意题意最后是要你求 x 的最小正整数解...

#include<iostream>
using namespace std;

int kzgcd(int a,int b,int &x,int &y){
    if(b==0){ x=1,y=0; return a; }
    int res = kzgcd(b,a%b,y,x);
    y -= a/b*x;
    return res;
}

int main(){
    int n,m;
    while(cin>>n>>m){
        int x,y;
        int flag=kzgcd(n,m,x,y);
        if(flag!=1) {
            printf("sorry\n");
            continue;
        }
        while(x<0){  对解的要求
            x += m;  注意解的转化
            y -= n;   
        }
        printf("%d %d\n",x,y);
    }
}