From 5b9155aabef4a72ebe123ef62dcfe13b71d00689 Mon Sep 17 00:00:00 2001
From: Zengtudor <Zengtudor@outlook.com>
Date: Wed, 5 Feb 2025 13:13:39 +0800
Subject: [PATCH] update

---
 .../index/P4549.cpp.952993C4FCC24847.idx      | Bin 0 -> 414 bytes
 src/2/P4549.cpp                               |  27 ++++
 src/2/day1.md                                 | 152 +++++++++++++++++-
 3 files changed, 178 insertions(+), 1 deletion(-)
 create mode 100644 .cache/clangd/index/P4549.cpp.952993C4FCC24847.idx
 create mode 100644 src/2/P4549.cpp

diff --git a/.cache/clangd/index/P4549.cpp.952993C4FCC24847.idx b/.cache/clangd/index/P4549.cpp.952993C4FCC24847.idx
new file mode 100644
index 0000000000000000000000000000000000000000..3a11b6487f52b6ea93877de28761ce18b592c357
GIT binary patch
literal 414
zcmWIYbaR`=$iU#7;#rZKT9U}Zz`!63#Kk2=nJ0mC1W>MGPVa<`S%)2X+P)gi?OF4q
zpiN|tNx%z%7uSBS3GQ9-hTF>e{CsZTg;T`1W-alrvA0)XN-Emh@p#T;p5V(2`zAc!
z6EunY<4>Q9%QkRVYn@xdeR~?q2i;kWMh{qKFvfb_o5MQ!mqT>nrZX2$l&Eam;d4cn
z;q)xFyjdHX@AE!hHtpYzyKfrfe%ID#yq4SN=@x6JYWnA+i2ajD-&1qfHrD;yCV5Z$
zqsJHftKv1r_qY}6IQ&ajyuN**Y{||wdWCJ4ub9H!k7XTWD6Y&+k_Y<8K>Uit<99`o
zJPZO1416psEF540C<6gSscFToKqgEbBNvK%MBp~3%C$Dj7+Lt(c-Z)0s{G0~Yj*x-
zSjWi2$HvLV0~2N_O3g_GnOj_xTr3V^@Pi02VB|a5?ZG?c;5=1Eeum`SlpH>wFc%LG
P6Ehnp3p+0>2O~EC7!Qye

literal 0
HcmV?d00001

diff --git a/src/2/P4549.cpp b/src/2/P4549.cpp
new file mode 100644
index 0000000..c04f78a
--- /dev/null
+++ b/src/2/P4549.cpp
@@ -0,0 +1,27 @@
+#include <cstdlib>
+#include <iostream>
+
+int gcd(int a,int b){
+    if(b==0)return a;
+    return  gcd(b,a%b);
+}
+
+template<class T>
+void abss(T &t){
+    t=abs(t);
+}
+
+int main(){
+    int n;
+    std::cin>>n;
+
+    int tmpa,tmpb;
+    std::cin>>tmpa;
+    abss(tmpa);
+    for(int i{1};i<n;i++){
+        std::cin>>tmpb;
+        abss(tmpb);
+        tmpa = gcd(tmpa,tmpb);
+    }
+    std::cout<<tmpa<<"\n";
+}
\ No newline at end of file
diff --git a/src/2/day1.md b/src/2/day1.md
index 64b63fd..4464fb7 100644
--- a/src/2/day1.md
+++ b/src/2/day1.md
@@ -34,4 +34,154 @@ GCD CRT 逆元
 3. 回溯到 $\gcd(30, 12)$，计算 $x = 1$, $y = -2$，满足 $30 \cdot 1 + 12 \cdot (-2) = 6$。
 
 ### 总结
-代码通过递归分解问题，逐步缩小规模。每次递归后，利用子问题的解 $x'$ 和 $y'$ 推导当前层的解 $x$ 和 $y$，最终得到满足 $ax + by = \gcd(a, b)$ 的整数解。
\ No newline at end of file
+代码通过递归分解问题，逐步缩小规模。每次递归后，利用子问题的解 $x'$ 和 $y'$ 推导当前层的解 $x$ 和 $y$，最终得到满足 $ax + by = \gcd(a, b)$ 的整数解。
+
+好的，下面是对这段 C++ 代码的逐行详细解释：
+
+```cpp
+int Exgcd(int a, int b, int &x, int &y) {
+```
+- 定义了一个名为 `Exgcd` 的函数。
+- 函数参数包括两个整数 `a` 和 `b`，以及两个引用参数 `x` 和 `y`。
+- `x` 和 `y` 用来输出结果，表示满足方程 $ ax + by = \text{GCD}(a, b) $ 的系数。
+
+```cpp
+  if (!b) {
+```
+- 检查 `b` 是否为 0。`!b` 表示 `b` 为 0 的时候，条件为真。
+
+```cpp
+    x = 1;
+    y = 0;
+    return a;
+```
+- 如果 `b` 是 0，则最大公约数是 `a`，同时根据贝祖定理，我们可以得到系数：
+  - 设置 `x = 1`（表示 $a \cdot 1$），
+  - 设置 `y = 0`（表示 $b \cdot 0$）。
+- 返回 `a`，也就是最大公约数。
+
+```cpp
+  }
+```
+- 结束了对 `b` 是否为 0 的检查。
+
+```cpp
+  int d = Exgcd(b, a % b, x, y);
+```
+- 递归调用 `Exgcd` 函数，参数为 `b` 和 `a % b`（即 `a` 除以 `b` 的余数）。
+- 这个递归调用会继续进行，直到 `b` 为 0，返回的最大公约数会被存储在 `d` 中，同时计算得到的 $x$ 和 $y$ 值也会在递归返回的过程中更新。
+
+```cpp
+  int t = x;
+```
+- 创建一个临时变量 `t`，存储 `x` 的当前值。这个值在接下来的代码中会被更新。
+
+```cpp
+  x = y;
+```
+- 将 `y` 的值赋给 `x`。在执行 `Exgcd(b, a % b)` 时，新的 $x$ 是旧的 $y$。
+
+```cpp
+  y = t - (a / b) * y;
+```
+- 更新 `y` 的值。根据公式：
+  - 我们把旧的 `y` 乘以 $ -\left(\frac{a}{b}\right) $（即取整运算向下取整），再减去 `t`（旧的 `x` 值）。
+- 这个步骤是用于反向计算贝祖系数。
+
+```cpp
+  return d;
+```
+- 返回最大公约数 `d`。
+
+整个函数的目的是找到两个数 $a$ 和 $b$ 的最大公约数，同时计算出系数 $x$ 和 $y$，使得 $ax + by = \text{GCD}(a, b)$ 成立。这个算法在许多数字和代数问题中都非常有用，尤其是在计算模反元素和求解线性同余方程时。
+
+好的！以下是逐行注释的 **扩展欧几里得算法（Exgcd）** 的 C++ 代码解释，帮助理解递归过程和解的推导逻辑：
+
+```cpp
+int Exgcd(int a, int b, int &x, int &y) {
+    // 递归终止条件：当 b = 0 时，gcd(a, 0) = a
+    if (!b) {
+        // 此时方程 ax + 0*y = a，显然 x = 1, y 可以是任意整数（此处设为 0）
+        x = 1;   // 基础解 x = 1
+        y = 0;   // 基础解 y = 0（为了后续递归回溯时计算方便）
+        return a; // 返回最大公约数
+    }
+
+    // 递归调用：计算 gcd(b, a % b)，并得到对应的解 x' 和 y'
+    // 注意：这里的 x 和 y 会被递归过程不断修改，最终得到的是下一层的解
+    int d = Exgcd(b, a % b, x, y);
+
+    // 保存递归返回后的 x（此时 x 已经是下一层的解 x'）
+    int t = x; 
+
+    // 当前层的 x 等于下一层的 y'
+    // 根据公式推导，x = y'
+    x = y; 
+
+    // 当前层的 y 等于下一层的 x' - (a / b) * y'
+    // 公式推导：y = x' - floor(a / b) * y'
+    y = t - (a / b) * y; 
+
+    // 返回最大公约数（递归过程中一直传递 d，最终返回给最外层）
+    return d;
+}
+```
+
+---
+
+### **关键步骤图解**
+假设输入 `a = 30`, `b = 12`，递归过程如下：
+
+1. **第一层递归** `Exgcd(30, 12, x, y)`:
+   - 计算 `a % b = 30 % 12 = 6`
+   - 进入第二层递归 `Exgcd(12, 6, x, y)`
+
+2. **第二层递归** `Exgcd(12, 6, x, y)`:
+   - 计算 `a % b = 12 % 6 = 0`
+   - 进入第三层递归 `Exgcd(6, 0, x, y)`
+
+3. **第三层递归** `Exgcd(6, 0, x, y)`:
+   - 终止条件触发：`x = 1`, `y = 0`, 返回 `d = 6`
+
+4. **回溯到第二层** `Exgcd(12, 6, x, y)`:
+   - `t = x = 1`（保存下一层的 x' = 1）
+   - `x = y = 0`（当前层 x = y' = 0）
+   - `y = t - (12/6)*y = 1 - 2*0 = 1`（当前层 y = 1）
+   - 返回 `d = 6`
+
+5. **回溯到第一层** `Exgcd(30, 12, x, y)`:
+   - `t = x = 0`（保存下一层的 x' = 0）
+   - `x = y = 1`（当前层 x = y' = 1）
+   - `y = t - (30/12)*y = 0 - 2*1 = -2`（当前层 y = -2）
+   - 返回 `d = 6`
+
+最终得到解：`x = 1`, `y = -2`，满足 `30*1 + 12*(-2) = 6`（即 gcd(30, 12) = 6）。
+
+---
+
+### **核心公式推导**
+递归过程基于以下恒等式：
+$$
+\gcd(a, b) = \gcd(b, a \% b)
+$$
+假设递归后得到的解为 $ x' $ 和 $ y' $，满足：
+$$
+b \cdot x' + (a \% b) \cdot y' = \gcd(b, a \% b)
+$$
+利用 $ a \% b = a - \lfloor a/b \rfloor \cdot b $，可重写为：
+$$
+a \cdot y' + b \cdot \left( x' - \lfloor a/b \rfloor \cdot y' \right) = \gcd(a, b)
+$$
+因此，当前层的解为：
+$$
+x = y', \quad y = x' - \lfloor a/b \rfloor \cdot y'
+$$
+
+---
+
+### **总结**
+- **递归逻辑**：通过不断缩小问题规模（`a → b`, `b → a%b`），直到 `b = 0` 时终止。
+- **解的回溯**：利用递归返回的解 `x'` 和 `y'`，推导当前层的解 `x` 和 `y`。
+- **变量交换**：通过临时变量 `t` 保存中间值，确保计算顺序正确。
+
+这个算法高效地求解了 **贝祖等式（Bézout's identity）** 的整数解，同时计算最大公约数。
\ No newline at end of file