引例:
求满足ax+by+c=0的x的最小正整数解;
扩展欧几里德算法利用了后半句结论;
证明可以在看完引例求解之后随便找一篇讲证明的文章就能看懂,因为思路是一致的,是引例求解思路的逆用,用数学归纳法证明,这里不赘述;
引例求解思路:
代码实现:
细节都在注释里了,上面理解了看完注释应该可以理解代码的写法
//扩展欧几里德,给定a,b,则一定存在整数x,y,使得ax+by=gcd(a,b);可以用于二元一次方程求整数解
//整个过程是,先得到了题目条件的a,b,然后不断向下exgcd,直到bn等于0,此时才满足了if(!b),给x,y赋值,
//此时跳出了最底层的exgcd,那么倒数第二层的exgcd这一行执行完毕,向下执行y -= x * (a / b),然后每一层往上都是这么执行的;
void Exgcd(int a, int b, int& d, int& x, int& y)//int& d这种写法是引用传参,类似于指针的原理可以对函数外的同名实参进行修改
{
if (!b)//因为gcd(a,b)=gcd(b,a%b),所以出口就是gcd(b,0)
{
d = a;//d就是a,b的最大公约数
x = 1;//xn=1
y = 0;//yn=0
//也就是xn=1,yn=0是回溯的起点
}