c++统计的方法求PI

使用蒙特卡洛方法可以估算圆周率π。该方法的基本思想是，在一个正方形内绘制一个圆，然后随机投掷大量的点到这个正方形中，统计落在圆内的点的数量与总点数的比例。这个比例将近似于圆的面积与正方形面积的比例，而这个比例又近似于π/4。

以下是一个简单的C++程序，演示了如何使用蒙特卡洛方法估算π：

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

int main() {
    std::srand(std::time(0)); // 使用当前时间作为随机数种子

    int totalPoints = 1000000; // 总点数
    int pointsInsideCircle = 0; // 圆内的点数

    for (int i = 0; i < totalPoints; ++i) {
        double x = static_cast<double>(std::rand()) / RAND_MAX;  // 在 [0, 1) 范围内生成随机 x 坐标
        double y = static_cast<double>(std::rand()) / RAND_MAX;  // 在 [0, 1) 范围内生成随机 y 坐标

        // 判断点是否在单位圆内
        if (x * x + y * y <= 1.0) {
            pointsInsideCircle++;
        }
    }

    // 使用圆内点数与总点数的比例估算 π
    double estimatedPi = 4.0 * static_cast<double>(pointsInsideCircle) / totalPoints;

    std::cout << estimatedPi << std::endl;

    return 0;
}

请注意，这只是一个估算，精确度取决于投掷的点数。增加totalPoints的值可以提高估算的准确性。