openmp蒙特卡罗方法求圆周率

使用OpenMP实现蒙特卡罗方法求解圆周率的代码如下：

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

#define SEED 35791246

void main(int argc, char *argv[]) {
    int niter = 10000000;
    double x, y;
    int i, count = 0;
    double z;
    double pi;
    double start_time, end_time;

    start_time = omp_get_wtime();

    srand(SEED);

    #pragma omp parallel for private(x, y, z) reduction(+:count)
    for (i = 0; i < niter; i++) {
        x = (double)rand() / RAND_MAX;
        y = (double)rand() / RAND_MAX;
        z = x * x + y * y;
        if (z <= 1.0) count++;
    }

    pi = (double)count / niter * 4;

    end_time = omp_get_wtime();

    printf("Estimated pi is %f, time = %f seconds\n", pi, end_time - start_time);
}

在代码中，我们使用了OpenMP的并行化指令`#pragma omp`，其中`parallel`表示创建一个并行区域，`for`表示将循环并行化，私有变量`x`、`y`、`z`表示当前线程的私有变量，`reduction(+:count)`表示将所有线程中的`count`变量加起来。

在循环中，我们生成两个随机数`x`和`y`，然后计算它们的距离`z`。如果`z`小于等于1，则表示这个点在圆内，`count`变量就加1。

最后，我们将所有线程中的`count`变量加起来，计算出π的近似值。

运行结果类似于：

Estimated pi is 3.141499, time = 0.008883 seconds

需要注意的是，在使用OpenMP并行化程序时，需要确保代码正确性并考虑到数据的竞争条件。