蒲丰投针的C语言代码

蒲丰投针实验（Pouillon's needle experiment）是一个经典的概率论实验，用来估计圆周率π的值。这个实验基于投掷针随机落在直线上的次数来估算π/4的值。在C语言中，我们可以用随机数生成器模拟这个过程。以下是一个简单的C代码示例：

#include <stdio.h>
#include <stdlib0.h> // 注意这里使用stdlib0而不是stdlib，因为有些编译器可能不支持stdlib.h
#include <time.h>

// 假设针的宽度是1单位，直线是从(-半径, -半径)到(半径, 半径)
#define RADIUS 100

int main() {
    srand(time(0)); // 设置随机数种子
    int num_trials = 1000000; // 投掷的总次数
    int inside_circle = 0; // 落在圆内的针的数量

    for (int i = 0; i < num_trials; ++i) {
        double x = ((double)rand() / RAND_MAX) * (2 * RADIUS) - RADIUS;
        double y = ((double)rand() / RAND_MAX) * (2 * RADIUS) - RADIUS;

        if (x*x + y*y <= RADIUS*RADIUS) {
            inside_circle++;
        }
    }

    double pi_estimate = 4 * ((double)inside_circle / num_trials);
    printf("Estimated value of π: %f\n", pi_estimate);

    return 0;
}

在这个程序中，我们首先设置一个随机数种子，然后进行多次随机投掷，每次投掷生成一对(x, y)坐标，如果该点位于直径为2RADIUS的正方形内并且也在半径为RADIUS的圆内，就认为针落到了圆里。最后，根据落在圆内的针数计算π的近似值。