话说题目是,用随机算法求PI,就是圆周那个PI。随机算法神马的从来都是用过,没写过,这里随便贴个线性同余法
![](https://i-blog.csdnimg.cn/blog_migrate/3f7db8baea3c96dcb466faf6a47ad760.jpeg)
其中b >= 0,c >= 0,d <= m。d称为该随机序列的种子。如何选取该方法中的常数b、c和m直接关系到所产生的随机序列的随机性能。这是随机性理论研究的内容,已超出本书讨论的范围。从直观上看,m应取得充分大,因此可取m为机器大数,另外应取gcd(m, b) = 1,因此可取b为一素数。
理论的东西,不折腾了,老实说我也不会,这里直接用c库里的rand,种子也不初期化了。
![image image](https://i-blog.csdnimg.cn/blog_migrate/be47d6670169278c3756fcb83553b50c.png)
原理是很简单了,以圆心为原点,往上面扔随机点,如果随机点的个数足够多,无疑圆心内的点数除以总点数 = (PI * pow(r,2)) / (4 * pow(r,2)) = PI / 4 。所以当n足够大时,k与n之比就逼近这一概率。从而,PI 约等于 (4*k)/n。看实现:
double get_rand_num() {
//以圆心作为原点,随机点的x轴,y轴的坐标
return rand() / RAND_MAX;
}
double get_pi(long num = 200000000) {
int k = 0;
for (int i=0; i<num; i++) {
double x = get_rand_num();
double y = get_rand_num();
if ((x*x + y*y) < 1.0) {
k++;
}
}
return (double)4*k / (double)num;
}
不过么,结果的精确度实在是不咋的,我的机器上跑出来3.9几,可能是计算中的误差问题,也肯能是点数不够大,。。。就这样吧。