用C语言帮我实现以下蒲丰投针实验

蒲丰投针实验是一个著名的概率论实验，用以估计π的值。在这个实验中，人们将很多根长度相等的针随机投到一张纸上，纸上画有间距相等的平行线。通过计算针与平行线的交叉次数，可以估算出π的值。其原理是几何概率的应用。

以下是一个简化版的C语言程序实现：

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>

// 定义针的长度和线之间的距离
#define NeedleLength 1.0
#define LineDistance NeedleLength

// 生成[0, 1]之间的随机数
double random_double() {
    return (double)rand() / RAND_MAX;
}

// 实现蒲丰投针实验
double buffonNeedleExperiment(int total, double lineDistance) {
    int hits = 0;
    double piEstimate = 0.0;

    for (int i = 0; i < total; i++) {
        // 生成两个随机数，分别代表针的中心位置和针与水平线的夹角
        double center = random_double() * lineDistance;
        double angle = random_double() * (M_PI / 2.0); // 生成[0, π/2]之间的随机数

        // 判断针是否与线相交
        if (center < NeedleLength / 2.0 || sin(angle) < (lineDistance * 2.0) / NeedleLength) {
            hits++;
        }
    }

    // 根据公式计算π的估计值
    piEstimate = 2 * NeedleLength * total / (hits * lineDistance);
    return piEstimate;
}

int main() {
    int totalNeedles; // 需要投掷的针的数量
    printf("请输入投掷的针的数量: ");
    scanf("%d", &totalNeedles);

    srand(time(NULL)); // 初始化随机数种子
    double estimatedPi = buffonNeedleExperiment(totalNeedles, LineDistance);
    printf("通过蒲丰投针实验估算得到的π值为: %f\n", estimatedPi);

    return 0;
}

这个程序首先定义了针的长度和线间距，然后通过循环模拟投针过程，每次投掷针时生成两个随机数，一个代表针的中心位置，另一个代表针与水平线的夹角。根据这个夹角判断针是否与线相交，并记录相交次数。最后，使用公式`π ≈ 2 * 针的长度 * 投掷次数 / (相交次数 * 线间距)`来估算π的值。