C++ 玩转数学 - L19 高级概念

📚 编程与数学的结合

为什么用编程解决数学问题？

编程可以自动化处理复杂的数学计算，让原本需要手工计算的 tedious 工作变得简单高效。

✅ 快速计算：计算机可以在毫秒级完成复杂运算
✅ 精确无误：避免人工计算的错误
✅ 批量处理：一次性处理大量数据
✅ 可视化：将抽象的数学概念直观呈现
✅ 算法思维：培养逻辑推理和问题分解能力

💡 学习价值：
• 加深对数学概念的理解
• 提高编程实践能力
• 培养算法思维
• 为竞赛和实际应用打基础

常用的数学库函数

#include <iostream>
#include <cmath>  // 数学函数库
using namespace std;

int main() {
    // 常用数学函数
    cout << "平方根：" << sqrt(16) << endl;        // 4
    cout << "幂运算：" << pow(2, 10) << endl;         // 1024
    cout << "绝对值：" << abs(-5) << endl;             // 5
    cout << "向上取整：" << ceil(3.2) << endl;          // 4
    cout << "向下取整：" << floor(3.8) << endl;         // 3
    cout << "四舍五入：" << round(3.5) << endl;         // 4
    cout << "π的值：" << M_PI << endl;                  // 3.14159...
    
    return 0;
}

💡 重要提示：
• 使用数学函数需要包含<cmath>头文件
• sqrt()：平方根
• pow(x, y)：x的y次方
• abs()：绝对值（整数）
• M_PI：圆周率常量

⚙️ 基础数学算法

判断质数（素数）

质数：只能被1和自身整除的大于1的自然数，如2、3、5、7、11...

#include <iostream>
using namespace std;

// 函数：判断是否为质数
bool isPrime(int n) {
    if (n <= 1) return false;
    if (n <= 3) return true;
    if (n % 2 == 0 || n % 3 == 0) return false;
    
    // 只需检查到√n
    for (int i = 5; i * i <= n; i += 6) {
        if (n % i == 0 || n % (i + 2) == 0) {
            return false;
        }
    }
    return true;
}

int main() {
    int num = 17;
    
    if (isPrime(num)) {
        cout << num << " 是质数" << endl;
    } else {
        cout << num << " 不是质数" << endl;
    }
    
    // 输出100以内的所有质数
    cout << "100以内的质数：";
    for (int i = 2; i <= 100; i++) {
        if (isPrime(i)) {
            cout << i << " ";
        }
    }
    cout << endl;
    
    return 0;
}

💡 优化技巧：
• 只需检查到√n，因为如果n有因子，必有一个≤√n
• 跳过偶数和3的倍数，每次增加6
• 时间复杂度从O(n)降到O(√n)

求最大公约数（GCD）

最大公约数：两个或多个整数共有约数中最大的一个。

#include <iostream>
using namespace std;

// 方法1：辗转相除法（欧几里得算法）
int gcd(int a, int b) {
    while (b != 0) {
        int temp = b;
        b = a % b;
        a = temp;
    }
    return a;
}

// 方法2：递归实现
int gcd_recursive(int a, int b) {
    if (b == 0) return a;
    return gcd_recursive(b, a % b);
}

int main() {
    int a = 48, b = 18;
    
    cout << a << " 和 " << b << " 的最大公约数：";
    cout << gcd(a, b) << endl;  // 6
    
    // 应用：化简分数
    int numerator = 48, denominator = 18;
    int divisor = gcd(numerator, denominator);
    cout << numerator << "/" << denominator << " = ";
    cout << numerator/divisor << "/" << denominator/divisor << endl;  // 8/3
    
    return 0;
}

求最小公倍数（LCM）

最小公倍数：两个或多个整数公有的倍数中最小的一个。

#include <iostream>
using namespace std;

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

// 利用公式：LCM(a,b) = a * b / GCD(a,b)
int lcm(int a, int b) {
    return a / gcd(a, b) * b;  // 先除后乘，避免溢出
}

int main() {
    int a = 12, b = 18;
    
    cout << a << " 和 " << b << " 的最小公倍数：";
    cout << lcm(a, b) << endl;  // 36
    
    return 0;
}

🎮 互动实验：数学计算器

输入数字，体验各种数学算法

数字1：数字2：

输入数字后点击"计算"...

🛠️ 经典数学序列

斐波那契数列

斐波那契数列：0, 1, 1, 2, 3, 5, 8, 13, 21...，每个数是前两个数之和。

#include <iostream>
using namespace std;

// 方法1：迭代法（推荐）
long long fibonacci_iterative(int n) {
    if (n <= 0) return 0;
    if (n == 1) return 1;
    
    long long a = 0, b = 1, c;
    for (int i = 2; i <= n; i++) {
        c = a + b;
        a = b;
        b = c;
    }
    return b;
}

// 方法2：递归法（简洁但效率低）
long long fibonacci_recursive(int n) {
    if (n <= 0) return 0;
    if (n == 1) return 1;
    return fibonacci_recursive(n - 1) + fibonacci_recursive(n - 2);
}

int main() {
    // 输出前20个斐波那契数
    cout << "前20个斐波那契数：";
    for (int i = 0; i < 20; i++) {
        cout << fibonacci_iterative(i) << " ";
    }
    cout << endl;
    
    // 计算第50个斐波那契数
    cout << "第50个斐波那契数：" << fibonacci_iterative(50) << endl;
    
    return 0;
}

💡 两种方法对比：
• 迭代法：时间复杂度O(n)，效率高，推荐使用
• 递归法：代码简洁，但时间复杂度O(2^n)，效率极低
• 实际应用中优先使用迭代法或动态规划

阶乘计算

阶乘：n! = 1 × 2 × 3 × ... × n

#include <iostream>
using namespace std;

// 方法1：迭代法
long long factorial_iterative(int n) {
    long long result = 1;
    for (int i = 2; i <= n; i++) {
        result *= i;
    }
    return result;
}

// 方法2：递归法
long long factorial_recursive(int n) {
    if (n <= 1) return 1;
    return n * factorial_recursive(n - 1);
}

int main() {
    // 计算1到10的阶乘
    for (int i = 1; i <= 10; i++) {
        cout << i << "! = " << factorial_iterative(i) << endl;
    }
    
    return 0;
}

判断完全数

完全数：所有真因子（除了自身以外的因子）之和等于自身的数。如6 = 1 + 2 + 3。

#include <iostream>
using namespace std;

bool isPerfect(int n) {
    if (n <= 1) return false;
    
    int sum = 1;  // 1是所有数的因子
    for (int i = 2; i * i <= n; i++) {
        if (n % i == 0) {
            sum += i;
            if (i != n / i) {  // 避免重复加平方根
                sum += n / i;
            }
        }
    }
    
    return sum == n;
}

int main() {
    cout << "10000以内的完全数：";
    for (int i = 2; i <= 10000; i++) {
        if (isPerfect(i)) {
            cout << i << " ";
        }
    }
    cout << endl;  // 6 28 496 8128
    
    return 0;
}

🎯 数学序列生成器

选择要生成的数学序列

项数：

选择序列类型后查看结果...

🚀 高级数学应用

实际应用：简单的计算器

#include <iostream>
#include <cmath>
using namespace std;

int main() {
    double a, b;
    char op;
    
    cout << "=== 简易计算器 ===" << endl;
    cout << "请输入：数字1 运算符 数字2" << endl;
    cout << "支持的运算符：+ - * / ^ s(平方根)" << endl;
    cout << "示例：5 + 3 或 s 16" << endl;
    cout << "输入：";
    
    cin >> op;
    
    if (op == 's' || op == 'S') {
        // 平方根运算
        cin >> a;
        if (a < 0) {
            cout << "错误：不能对负数开平方" << endl;
        } else {
            cout << "sqrt(" << a << ") = " << sqrt(a) << endl;
        }
    } else {
        // 二元运算
        cin >> a >> b;
        
        switch (op) {
            case '+':
                cout << a << " + " << b << " = " << (a + b) << endl;
                break;
            case '-':
                cout << a << " - " << b << " = " << (a - b) << endl;
                break;
            case '*':
                cout << a << " * " << b << " = " << (a * b) << endl;
                break;
            case '/':
                if (b == 0) {
                    cout << "错误：除数不能为0" << endl;
                } else {
                    cout << a << " / " << b << " = " << (a / b) << endl;
                }
                break;
            case '^':
                cout << a << " ^ " << b << " = " << pow(a, b) << endl;
                break;
            default:
                cout << "不支持的运算符" << endl;
        }
    }
    
    return 0;
}

数学在编程中的应用场景

数学知识在实际编程中有广泛应用：

🎮 游戏开发：碰撞检测、物理引擎、图形变换
📊 数据分析：统计计算、概率分布、回归分析
🔐 密码学：质数、模运算、加密算法
🤖 人工智能：线性代数、微积分、优化算法
🌐 图形学：几何计算、矩阵变换、光照模型
💰 金融计算：复利计算、风险评估、期权定价

性能优化技巧

// 技巧1：预计算并缓存结果
const int MAX_N = 1000;
long long fib_cache[MAX_N];

long long fibonacci_memo(int n) {
    if (n <= 1) return n;
    if (fib_cache[n] != 0) return fib_cache[n];  // 已计算过
    fib_cache[n] = fibonacci_memo(n - 1) + fibonacci_memo(n - 2);
    return fib_cache[n];
}

// 技巧2：使用位运算优化
// 判断奇偶：n % 2 等价于 n & 1
// 乘以2：n * 2 等价于 n << 1
// 除以2：n / 2 等价于 n >> 1

// 技巧3：减少重复计算
double result = sqrt(x);  // 只计算一次
for (int i = 0; i < n; i++) {
    // 使用result，而不是每次都调用sqrt(x)
}

📝 实战练习

📝 理解测试

以下代码的输出是什么？

#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() {
    cout << gcd(36, 24) << endl;
    return 0;
}

选择答案查看解析...

💻 综合挑战

题目：创建一个程序，找出1000以内的所有孪生素数对。孪生素数是指相差2的两个质数，如(3,5)、(5,7)、(11,13)。

点击"生成代码"查看完整实现...