蒙特卡洛算法,简而言之,就是初中学的概率;对于一个圆的面积,我们设置一个正方形,然后不断向里面抛入点,那个点落在圆内的概率 乘以正方形的面积,就是圆的面积!简单吧,哈哈,上个图给你们理解一下:
#include<iostream>
#include<stdio.h>
#include<stdlib.h>
#include<time.h>
using namespace std;
#define N 50000
//抛出点设置为5万个,当然,点越多,那么结果越精确
//这个地方,由于是电脑自动产生的随机数,每次结果不一样,可能差的有一点点大,那么你就把抛出点设置多一点
int main() {
//sum是点在圆内的个数;
int sum = 0;
//x,y是坐标
double x, y;
srand(time(NULL));
for (int i = 0; i < N; i++) {
//RAND_MAX是 0-1.0的浮点数,自带的
x = (double)rand() / RAND_MAX;
y = (double)rand() / RAND_MAX;
if ((x*x + y * y) < 1) {
sum++;
}
}
cout << sum << endl;
//这里一定要小心,输出值请转化我浮点数,不然就是3!
cout << "圆的面积:" << (double)4 * sum / N << endl;
system("pause");
return 0;
}