蒙特卡洛法多线程求圆周率

最新推荐文章于 2025-03-27 16:30:45 发布

我说我是程序猿

最新推荐文章于 2025-03-27 16:30:45 发布

阅读量3.5k

点赞数

分类专栏：一些代码文章标签： c++ 多线程圆周率

本文链接：https://blog.youkuaiyun.com/yy64578537/article/details/51543193

版权

这是一个C++程序，通过多线程应用蒙特卡洛方法来估算圆周率。程序首先定义了最大迭代次数MAX，然后创建指定数量的线程。每个线程负责计算一部分迭代，确保在并发环境中对共享变量result进行互斥访问。最后，程序输出计算得到的π值及其计算所用的时间。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<iostream>
#include<string>
#include<pthread.h>
#include<cmath>
#include<cstdlib>
#include<iomanip>
#include<cstdio>
#include<sys/time.h>
#define MAX 100000000
using namespace std;
double result = 0;
int thread_num;
int Every_P;
pthread_mutex_t mut;
void *t_hread(void *arg){
    int m = *(int*)arg;
    if(m != thread_num - 1){
        for(int i = m * Every_P;i < (m + 1) * Every_P;i++){
	    	double x = sqrt(MAX) * sqrt(MAX - ((i + 0.5)/MAX)*(i + 0.5));
	    	int xx = x;
	    	double count;
	    	if((x - xx) > 0.5){
				count = xx + 1; 		
	    	}
	    	else{
				count = xx;
	    	}
	    	pthread_mutex_lock(&mut);
	    	result = result + (double)count/MAX;
	    	pthread_mutex_unlock(&mut);							
		}	
    }
    if(m == thread_num - 1){
		for(int i = m * Every_P;i < MAX;i++){
	   		double x = sqrt(M

最低0.47元/天解锁文章