看完数论数的这一段,觉得这个算法太精辟了,其实扩展欧几里德算法就是书上的公式翻译成代码,实在是太棒了,我还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);
}
}