C++ 概率算法 利用蒙特卡羅算法計算圓周率

概率算法大致可分為4種形式：

數值概率算法；

蒙特卡羅算法；

拉斯維加斯算法；

舍伍德算法；

計算蒙特卡羅概率的算法實現：

#include "stdio.h"
#include "time.h"
#include "stdlib.h"

double MontePI(int n)
{
    double PI;
    double x, y;
    int i, sum;

    sum = 0;
    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++;
    }
    PI = 4.0*sum / n;
    return PI;
}

int main()
{
    int n;
    double PI;

    scanf("%d",&n);
    PI = MontePI(n);
    printf("PI=%f\n",PI);
    system("pause");
    return 0;
}

注意 思想 的學習。