4.用C语言和OpenMP编写Π值的计算：1-\frac{1}{3}+\frac{1}{5}-\frac{1}{7}+\cdots=\frac{\pi}{4}

下面是使用C语言和OpenMP并行计算Π值的示例代码：

#include <stdio.h>
#include <omp.h>

int main() {
    int n = 1000000000; // 计算n项
    double pi = 0.0;
    double step = 1.0 / (double) n;

    #pragma omp parallel for reduction(+:pi)
    for (int i = 0; i < n; i++) {
        double x = (i + 0.5) * step;
        pi += 4.0 / (1.0 + x * x);
    }

    pi *= step;

    printf("π ≈ %f\n", pi);

    return 0;
}

这段代码使用了OpenMP的并行for指令，将n个任务并行执行。每个任务计算一个矩形的面积，并将结果累加到pi变量中。最后根据公式π/4 = 1 - 1/3 + 1/5 - 1/7 + ... 计算π的值。使用reduction(+:pi)指令将pi变量的值在所有线程中进行求和。