#include <iostream>
#include <stdlib.h>
#include <ctime>
using namespace std;
//创建长浮点型函数使用蒙特卡洛算法计算圆周率
long double jisuanpi(int *n){
//定义坐标x,y均为长浮点型,以提高精度
long double x,y;
//定义圆周率为长浮点型数
long double pi;
//定义一个临时变量用来接收落在圆形区域内的点的个数
int sum=0;
//创建时间随机种子
srand(time(NULL));
//模拟算法在区域内撒点
for(int i=0;i<*n;i++){
x=(long double)rand()/RAND_MAX;
y=(long double)rand()/RAND_MAX;
//如果满足圆形区域的条件,则加1
if(x*x+y*y<=1){
sum++;
}
}
//pi的计算方式
pi=(4.0*sum)/(*n);
return pi;
}
int main(){
int n=0;
cout<<"蒙罗卡特算法计算圆周率:"<<endl;
cout<<"请输入生成的随机数个数:"<<endl;
cin>>n;
cin.get();
cout<<"模拟得到的圆周率为:"<<jisuanpi(&n)<<endl;
cin.get();
}
蒙洛卡特算法近似计算圆周率pi
最新推荐文章于 2024-08-22 17:27:46 发布