% 辗转相处法求最大公因数
% 定理:gcd(a,b) = gcd(b,a mod b)
% 证明
% a可以表示成a = kb + r,则r = a mod b
% 假设d是a,b的一个公约数,则有
% d|a, d|b,而r = a - kb,因此d|r
% 因此d是(b,a mod b)的公约数
%
% 假设d 是(b,a mod b)的公约数,则
% d | b , d |r ,但是a = kb +r
% 因此d也是(a,b)的公约数
%
% 因此(a,b)和(b,a mod b)的公约数是一样的,其最大公约数也必然相等,得证
function nResult = my_gcd(m, n)
% 调整次序
if m < n, swap(m, n),end
% 第一次
nResult = mod(m, n);
last_n = n;
last_n = nResult;
nResult = mod(n, nResult);
n = last_n;
end
end
end
% 定理:gcd(a,b) = gcd(b,a mod b)
% 证明
% a可以表示成a = kb + r,则r = a mod b
% 假设d是a,b的一个公约数,则有
% d|a, d|b,而r = a - kb,因此d|r
% 因此d是(b,a mod b)的公约数
%
% 假设d 是(b,a mod b)的公约数,则
% d | b , d |r ,但是a = kb +r
% 因此d也是(a,b)的公约数
%
% 因此(a,b)和(b,a mod b)的公约数是一样的,其最大公约数也必然相等,得证
function nResult = my_gcd(m, n)
% 调整次序
if m < n, swap(m, n),end
% 第一次
nResult = mod(m, n);
last_n = n;
% 不使用递归
last_n = nResult;
nResult = mod(n, nResult);
n = last_n;
end
if nResult == 1
% 没有最大公因数,1为他们的最大公因数
else
% 找到
end
end