GCD和EGCD

#include <iostream>

using namespace std;

int gcd(int a, int b)
 {
    while (b != 0)
 {
        int temp = b;
        b = a % b;
        a = temp;
    }
    return a;
}

int main() {
    int num1, num2;
    cout << "请输入两个整数: ";
    cin >> num1 >> num2;
    int result = gcd(num1, num2);
    cout << "最大公约数为: " << result << endl;
    return 0;
}

gcd算法如上

二

#include <iostream>
using namespace std;
int egcd(int a, int b, int& r, int& s, int& d) 
{
    if (b == 0) 
    {
        r = 1;
        s = 0;
        d = a;
        return d;
    }
    else 
    {
        int r1, s1;
        d = egcd(b, a % b, r1, s1, d);
        r = s1;
        s = r1 - (a / b) * s1;
        return d;
    }
}

int main() 
{
    int a, b, r, s, d;
    cout << "请输入两个整数a和b：" << endl;
    cin >> a >> b;
    egcd(a, b, r, s, d);
    cout << "r = " << r << ", s = " << s << ", d = " << d << endl;
    return 0;
}

egcd算法如上

算法的递归步骤如下：

1.如果 b等于 0，那么结果就是 a，此时 r 被设为 1，s 被设为 0，即 ar + 0s = a。

2.否则，递归调用 egcd 函数，将 b 和 a % b 作为新的参数。这样会一直进行gcd操作，直到求出a和b的最大公因数，最终得到的返回值 d 就是 a 和 b 的最大公因数，而r1和s1就被赋为了1和0.

3.在递归返回后，根据递推关系式 ar + bs = d，重新为 r 和 s 赋值：r = s1，s = r1 - (a / b) * s1，r1起始值为1，s1起始值为0，递归n次就进行n次赋值运算，最终得出的r和s就是系数.

原理如下：

假设在某个递归调用中，我们已经计算得到了参数为 (a', b') 的 r1 和 s1，其中 a' = b， b' = a % b。我们可以将原始的递归关系式 ar + bs = d 应用于 (a', b')，得到：

a' * r1 + b' * s1 = d

将 a' = b 和 b' = a % b 替换到上式中，得到：

b * r1 + (a % b) * s1 = d

我们希望找到一对整数 r 和 s，使得上式依然成立，并满足 ar + bs = d。为了达到这个目标，我们可以将上式中的 a % b 替换为 a - (a // b) * b，其中 // 表示整数除法。

将替换后的表达式代入上式，整理之后可以得到：

b * r1 + (a - (a // b) * b) * s1 = d

进一步整理，得到：

a * s1 + b * (r1 - (a // b) * s1) = d

从上式可以看出，s = s1，而 r = r1 - (a // b) * s1。这样，我们就找到了新的 r 和 s 的值，满足递归关系式，并继续进行下一步的递归调用，直到满足终止条件。

因此，层层递推计算 r1 和 s1 的原理就是基于递归调用的展开和替代，并根据递推关系式确定每一层递归中 r1 和 s1 的值，最终得到满足条件的 r 和 s。

三

证明存在性：

构造集合

S = {a − bk : k ∈ Z 且 a − bk ≥ 0}.

显然，集合 S 非空。由良序原则，存在一个最小元 r ∈ S，且 r = a − qb。因此，

a = qb + r, r ≥ 0。

若r>=b,则可设r=b+m,其中m ∈ Z 。则a=qb+b+m=(q+1)b+m。不符合商和余数的原则。

证明唯一性：

假设存在不止一个整数对满足条件，则可设：

a=q1b+r1

a=q2b+r2

不妨设q1>q2,显然可得r1<r2.

q1b+r1=q2b+r2

(q1-q2)*b=r2-r1

b=(r2-r1)/(q1-q2)

因为0<=r1<b,0<=r2<b,r1<r2;所以0<r2-r1<b.

因为q1,q2为整数，q1>q2;所以q1-q2>=1.

因为0<r2-r1<b,q1-q2>=1,所以(r2-r1)/(q1-q2)<b,与b=(r2-r1)/(q1-q2)矛盾。

因此，假设不成立。唯一性得证。

四