Algorithm Gossip(14) 蒙地卡罗法求PI

    xiaoxiao2021-04-14  64

    前言

    This Series aritcles are all based on the book 《经典算法大全》; 对于该书的所有案例进行一个探究和拓展,并且用python和C++进行实现; 目的是熟悉常用算法过程中的技巧和逻辑拓展。

    提出问题

    14.Algorithm Gossip:求 蒙地卡罗法求 PI 送分算法题, 非常简单, 就是一个正方形,里面内切一个比例, 随机描点,问在圆里面的点的比例 得到一个值 k (面积比)。 那么我们可以知道 (k = PI /4 ). 所以 PI = 4*k;

    分析和解释

    代码

    #include <stdio.h> #include <stdlib.h> #include <time.h> #define N 50000 int main(void) { int i, sum = 0; double x, y; srand(time(NULL)); for(i = 1; i < N;i++) { x = (double)rand()/ RAND_MAX; y = (double)rand()/ RAND_MAX; if((x * x + y * y) < 1) sum++; } printf("PI = %f\n", (double)4 * sum / N); return 0; }

    拓展和关联

    Buffon 投针试验

    布丰投针试验和这个原理类似, 有兴趣可以看看, 也是测算 PI

    后记

    参考书籍

    《经典算法大全》维基百科
    转载请注明原文地址: https://ju.6miu.com/read-670286.html

    最新回复(0)