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;
}

注意 思想 的学习。