GCD和LCM

本文无特殊说明,字母代表的都为整数

公约数定义

若d|a,d|b。则d为a和b的公约数
在d的集合中,最大的d称为最大公约数,表示为(a,b)或gcd(a,b),最大公约数为正整数

公倍数定义

若a|c且b|c,则c为a和b的公倍数
在c的集合中,最大的c称为最大公倍数,表示为[a,b]货lcm(a,b),最大公倍数为正整数

互质(互素)

若(a,b)=1,则a与b互质

公约数性质

由于下列性质很显然,在这里不给出证明

1.(±a±b)=(a,b)

2.(a,b)=(b,a)

3.(a,b+ax)=(a,b)

4.(a,0)=|a|

5.(a,ka)=|a|

6.(ma,mb)=m(a,b) (m>0)

7.(a,b)=d>0,则(a/d,b/d)=1

8.若b|ac且(b,c)=1,则b|a

裴蜀定理

当a,b不全为0,则有x,y,使ax+by=(a,b)
集合A={q|q=ax+by}中最小正元素为(a,b)

1.(a,b)=1,则ax+by=1

2.若m|a,m|b,则m|(a,b)。即a和b的任一公约数都是gcd(a,b)的一个约数。

公倍数性质

1.a,b的任一公倍数都是lcm(a,b)的倍数

2.(a,b)[a,b]=|ab|

玄学

更相减损术

(a,b)=(b-a,a)
(a,b,c)=(a,b-a,c-b)

欧几里得算法

(a,b)=(b,a%b)
用于写gcd:

inline int gcd(int a,int b){
    return b == 0 ? a : gcd(b,a % b);
}

虽然gcd一般使用欧几里得算法,但是在数据大的时候,即使用了高精度来存储数据时,鉴于高精度算法中模运算十分难写,往往用更相减损术

欧几里得算法二进制形态

1.若x=y,gcd(x,y)=x

2.若x,y为偶数,gcd(x,y)=2*gcd(x/2,y/2)

3.若x为偶数,y为奇数,gcd(x,y)=gcd(x/2,y)

4.若x,y为奇数,gcd(x,y)=gcd(x-y,y)

<think>我们正在讨论C++标准库中计算最大公约数(GCD最小公倍数(LCM)的函数。 在C++17标准中,<numeric>头文件引入了两个函数:`std::gcd``std::lcm`。 注意:在C++17之前,标准库中没有直接提供这两个函数,但我们可以自己实现或者使用其他库(如Boost)。 1. `std::gcd`:计算两个整数的最大公约数。 2. `std::lcm`:计算两个整数的最小公倍数。 使用示例: 需要包含头文件:`#include <numeric>` 注意:这两个函数要求参数是非负整数,如果传入负数,结果是未定义的(undefined behavior)。 下面是一个简单的示例代码: ```cpp #include <iostream> #include <numeric> // 包含std::gcdstd::lcm int main() { int a = 12, b = 18; // 计算最大公约数 int gcd_value = std::gcd(a, b); std::cout << "GCD of " << a << " and " << b << " is " << gcd_value << std::endl; // 计算最小公倍数 int lcm_value = std::lcm(a, b); std::cout << "LCM of " << a << " and " << b << " is " << lcm_value << std::endl; return 0; } ``` 输出: GCD of 12 and 18 is 6 LCM of 12 and 18 is 36 注意事项: - 如果其中一个参数为0,则gcd返回另一个数的绝对值(但注意,因为要求非负,所以实际上如果传入0非负数,则返回那个非负数;如果两个都是0,则返回0)。但是,根据标准,如果两个参数都为0,则返回0。 - 对于lcm,如果有一个参数为0,则返回0(因为0任何数的最小公倍数是0)。但是,如果两个参数都为0,则标准规定返回0。 另外,在C++17之前,我们可以自己实现gcdlcm。例如,使用欧几里得算法: ```cpp // 最大公约数(非负整数) int gcd(int a, int b) { if (b == 0) return a; return gcd(b, a % b); } // 最小公倍数(非负整数) int lcm(int a, int b) { if (a == 0 || b == 0) { return 0; } return (a / gcd(a, b)) * b; // 先除后乘,避免溢出 } ``` 但是,现在如果使用C++17或更高版本,推荐直接使用标准库函数。 另外,注意这些函数有多个重载,支持有符号整数类型(如int, long, long long等),但要求参数是非负的(实际上,标准规定如果传入负数,行为未定义)。因此,在调用前,请确保参数是非负的,或者对负数取绝对值(但要注意,如果取绝对值可能会溢出,特别是对于最小负整数的情况)。 引用说明:C++标准库文档[^1]。 注意:上面的引用[^1]是之前用户提供的引用,但显然与C++标准库无关。因此,这里我们引用C++标准库的文档(例如cppreference.com)作为参考。但由于我们无法获取外部链接,所以这里仅作说明。 因此,在实际编程中,如果需要使用这些函数,请确保你的编译器支持C++17或更高版本。 如果编译器不支持C++17,那么可以自己实现(如上所示)或者使用其他库(如Boost)。 下面,我们生成相关问题。</think>### C++标准库中计算最大公约数最小公倍数的函数 在C++17及以上版本的标准库中,`<numeric>`头文件提供了两个专用函数: 1. **最大公约数 (GCD)** `std::gcd(m, n)` 计算两个整数 $m$ $n$ 的最大公约数,满足 $\gcd(m, n) = \max\{d \in \mathbb{Z} \mid d \mid m \text{ 且 } d \mid n\}$ - 要求:$m$ $n$ 必须是非负整数 - 示例:$\gcd(12, 18) = 6$ 2. **最小公倍数 (LCM)** `std::lcm(m, n)` 计算两个整数 $m$ $n$ 的最小公倍数,满足 $\text{lcm}(m, n) = \frac{|m \times n|}{\gcd(m, n)}$ - 要求:$m$ $n$ 必须是非负整数 - 示例:$\text{lcm}(12, 18) = 36$ #### 使用示例代码 ```cpp #include <iostream> #include <numeric> // 包含 std::gcd std::lcm int main() { int a = 12, b = 18; // 计算最大公约数 int gcd_val = std::gcd(a, b); std::cout << "GCD: " << gcd_val << "\n"; // 输出 6 // 计算最小公倍数 int lcm_val = std::lcm(a, b); std::cout << "LCM: " << lcm_val << "\n"; // 输出 36 return 0; } ``` #### 注意事项 1. **编译器支持**:需使用支持 C++17 或更高版本的编译器(如 GCC 7+、Clang 5+、MSVC 19.10+),编译时添加 flag `-std=c++17` 2. **负数处理**:传入负数会导致未定义行为(UB),需提前确保参数非负 3. **溢出风险**:$\text{lcm}(a,b)$ 计算时可能溢出,建议使用大整数类型(如 `long long`)处理大数 4. **C++17 以下版本**:需手动实现算法,例如: ```cpp // 欧几里得算法实现 GCD int gcd(int a, int b) { return b == 0 ? a : gcd(b, a % b); } // 通过 GCD 计算 LCM int lcm(int a, int b) { return a / gcd(a, b) * b; // 先除后乘避免溢出 } ``` #### 数学原理 - **GCD 性质**:$\gcd(a, b) = \gcd(b, a \mod b)$(欧几里得算法)[^1] - **LCM 性质**:$\text{lcm}(a, b) \times \gcd(a, b) = |a \times b|$[^2] --- ### 相关问题 1. 如何手动实现欧几里得算法计算最大公约数? 2. 当处理大整数时,如何避免计算最小公倍数时的溢出问题? 3. C++标准库中还有哪些常用的数值计算函数? 4. 为什么 `std::lcm` 要求参数必须是非负整数? 5. 在C++11/14中如何安全地计算GCDLCM? [^1]: 欧几里得算法基于递归关系 $\gcd(a, b) = \gcd(b, a \mod b)$,直到余数为0。 [^2]: 该恒等式是数论基本定理的直接推论。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值