給定兩個正整數m和n,我們計算它們的最大公因子d和兩個整數a和b,使得a*m+b*n=d

算法流程
　　E1.置a'=b=1;a=b'=0;c=m,d=n;

　　E2.計算d和r,使得c=q*d+r;

　　E3.若r==0;則退出,當前已有a*m+b*n=d;

　　E4;c=d;d=r;t=a';a'=a;a=t-q*a;t=b';b'=b;b=t-q*b;返回E2.

證明

　　對于已有的m和n,假設m>n;如果刨除變量a,b,a',b';算法與歐幾里得算法完全一樣,為計算最大公約數的算法.

　　最終要求的為a*m+b*n=d=GCD(m,n);如果改式子成立由歐幾里得算法可推出a'*n+b'*(m%n)=GCD(n,m%n);

　　因為GCD(m,n)=GCD(n,m%n);

　　所以a*m+b*n=a'*n+b'*(m%n)

　　　　　　　　=a'*n+b'*(m-(m/n)*n)

　　　　　　　　=a'*n+b'*m-b'*(m/n)*n

　　　　　　　　=b'*m+(a'-b'*(m/n))*n

　　所以a=b';b=a'-b'*(m/n);

　　可以推出根據a‘、b'可以計算a、b。

代碼實現