c++用蒙特卡洛方法计算圆周率

使用蒙特卡洛方法计算圆周率的基本思路是，假设有一个正方形，其内切一个半径为1的圆形，然后在正方形内随机生成若干个点，并计算有多少个点落在圆内，最后根据公式 π = 4 * 圆内点数 / 总点数 计算圆周率。

以下是使用 C++ 实现该算法的代码示例：

#include <iostream>
#include <cmath>
#include <cstdlib>
#include <ctime>

using namespace std;

int main()
{
    // 设置随机数种子
    srand(time(NULL));

    int total_points = 1000000; // 总点数
    int inside_points = 0; // 圆内点数

    for (int i = 0; i < total_points; i++) {
        // 在正方形内随机生成点
        double x = (double)rand() / RAND_MAX;
        double y = (double)rand() / RAND_MAX;

        // 判断该点是否在圆内
        if (sqrt(x * x + y * y) <= 1) {
            inside_points++;
        }
    }

    // 计算圆周率
    double pi = 4.0 * inside_points / total_points;

    cout << "圆周率的近似值为：" << pi << endl;

    return 0;
}

注意，在实际应用中，我们需要通过增加随机点的数量来提高计算精度。