扩展欧几里得算法【原理+编码】

一、问题描述

已知两个整数m和n，求mx + ny = gcd(m, n)。

二、知识引入

贝祖定理：m>=0, n>=0, 且m,nZ    =>   E x, y Z, 使 mx + ny =gcd(m,n)

定理证明：1）设m=n=0，则gcd(m, n)=0，显然定理成立。
                  2）设m != 0 ，集合A={ mx+ny | x, y   Z }，则 A Z.
                        =>   令集合P=AN，则P={z | z >= 0 }. 【证明P】
                        =>   设d是P中最小的正整数，则 E x0, y0 Z , 使 mx0 + ny0 = gcd(m, n) = d.
                        =>   设d|m = q...r，即m=d*q+r,  d>r>=0.
                        =>   r = m - d * q = m - (mx0 + ny0)*q = m(1-x0) - n(qy0) 
                        =>   r   P, 又d>r与d是P中最小的正整数矛盾
                        =>   r = 0.
                        =>   m=d*q, 同理可证n=d*p，即d是m和n的公约数