蒙特卡洛法串行求解Pi值C++程序

以下是串行求解 Pi 值的 C++ 代码，使用蒙特卡洛方法：

#include <iostream>
#include <random>

int main() {
    const int num_points = 1000000;    // 模拟点数
    int num_in_circle = 0;             // 圆内点数

    std::random_device rd;
    std::mt19937 gen(rd());
    std::uniform_real_distribution<> dis(-1.0, 1.0);

    for (int i = 0; i < num_points; ++i) {
        double x = dis(gen);
        double y = dis(gen);
        if (x * x + y * y <= 1) {
            ++num_in_circle;
        }
    }

    double pi = 4.0 * num_in_circle / num_points;
    std::cout << "Pi = " << pi << std::endl;

    return 0;
}

该程序使用 `std::random_device` 和 `std::mt19937` 生成器生成均匀分布在 $[-1,1]$ 区间内的随机数，模拟在单位正方形内随机投点，判断每个点是否在以原点为圆心、半径为 $1$ 的圆内，最后根据圆内点数和总点数的比例计算 Pi 值。