C++编程用蒙特卡洛方法计算圆周率。

蒙特卡洛方法是一种基于随机抽样的统计方法，可以用来估计圆周率。具体实现步骤如下：

1. 首先，我们在一个正方形内部画一个圆，使得圆心位于正方形的中心，半径为正方形边长的一半。

2. 接着，我们随机生成 N 个点，这些点的 x 和 y 坐标都在正方形内部均匀分布。

3. 对于每一个点，我们计算其到圆心的距离 d，如果 d 小于等于圆的半径，则该点位于圆内。

4. 最后，我们统计落在圆内的点数 M，以及总共生成的点数 N，根据圆周率的定义，圆的面积为 π * r^2，其中 r 为圆的半径，正方形的面积为 (2r)^2，因此圆周率的估计值为 4 * M / N。

下面是 C++ 代码实现：

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

using namespace std;

int main()
{
    int N = 1000000; // 生成的点数
    int M = 0; // 落在圆内的点数
    double r = 0.5; // 圆的半径
    srand(time(NULL)); // 初始化随机数生成器

    for (int i = 0; i < N; i++) {
        double x = (double)rand() / RAND_MAX; // x 坐标在 [0,1) 内随机生成
        double y = (double)rand() / RAND_MAX; // y 坐标在 [0,1) 内随机生成
        double d = sqrt((x - 0.5) * (x - 0.5) + (y - 0.5) * (y - 0.5)); // 计算到圆心的距离

        if (d <= r) {
            M++;
        }
    }

    double pi = 4.0 * M / N; // 计算圆周率的估计值
    cout << "圆周率的估计值为：" << pi << endl;

    return 0;
}

注意，这种方法只是一种估计圆周率的方法，其精度取决于生成的随机数数量。当 N 越大时，估计值越接近真实值，但同时计算时间也会增加。