gcd(a,b)=gcd(b,a%b);
终止条件a=gcd b=0;
(gcd为a,b的最大公约数)
扩展欧几里得算法: a 和 b 的最大公约数是 gcd ,一定能够找到这样的 x 和 y ,使得: a*x + b*y = gcd 成立
我们只需要找到特殊解x0,y0;
则通解为 x = x0 + (b/gcd)*t y = y0 – (a/gcd)*t
那如何求出下一组解呢
仿照欧几里得算法a=b,b=a%b代入。
a%b = a - (a/b)*b(这里的 “/” 指的是整除)可以得到:
gcd = b*x1 + (a-(a/b)*b)*y1
= b*x1 + a*y1 – (a/b)*b*y1
= a*y1 + b*(x1 – a/b*y1)
对比易得 x = y1
y = x1 – a/b*y1
终止为 a*1 + b*0 = gcd即a的系数为1,b的系数为0或者其他值(*0=0)
我们再看过河这个问题 设a=s,b=s+1,ax+by=n a,b的最大公约数为1因为是相邻的自然数
ax+by=n是由ax+by=1(1为最大公约数)de系数*n得到的。(a=s,b=s+1)
设根分别为x0,y0
令x1=x0+bt;
y1=y0-at;
证明存在整数解x1≥0,y1≥0即可
x0+bt≥0->t≥-x0/b
y0-at≥0->t≤y0/a
综合-x0/b≤t≤y0/a
等于证明y0/a-(-x0/b)≥1
化简得n/ab≥1
原命题得证,存在整数解
下面贴上老师给我写的过程图
未完待续。