用C++写一个计算圆周率的代码

以下是使用蒙特卡罗方法计算圆周率的C++代码：

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

using namespace std;

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

    int n = 1000000; // 抛点次数
    int cnt = 0; // 落在圆内的点数

    for (int i = 0; i < n; i++)
    {
        double x = (double)rand() / RAND_MAX * 2.0 - 1.0; // 生成随机点的横坐标
        double y = (double)rand() / RAND_MAX * 2.0 - 1.0; // 生成随机点的纵坐标

        if (x * x + y * y <= 1.0) // 判断该随机点是否落在单位圆内
            cnt++;
    }

    double pi = 4.0 * cnt / n; // 计算圆周率

    cout << "Pi = " << pi << endl;

    return 0;
}

该程序通过随机抛点的方式，统计落在单位圆内的点数，从而估计圆的面积。根据圆的面积公式 $A=\pi r^2$，可以计算出圆周率 $\pi$ 的近似值。在上述代码中，抛点次数为 1000000 次，可以根据需要适当调整。