使用蒙地卡拉法求PI

        最近公司没有安排什么事做，于是研究了一些算法，其中求PI的方法值得思考。

        蒙地卡罗的基本原理为以乱数配合面积公式来进行解题，这种以机率来解题的方式带有赌博的意味，虽然在精确度上有所疑虑，但其解题的思考方向却是个值得学习的方式。

        蒙地卡罗的解法适用于与面积有关的题目，例如求PI值或椭圆面积，这边介绍如何求PI值；假设有一个圆半径为1，所以四分之一圆面积就为PI/4，而包括此四分之一圆的正方形面积就为1,于是便可以得到PI=c*4/n,如果随意的在正方形中投射飞标（点），则这些飞标（点）有些会落于四分之一圆内，假设所投射的飞标（点）有n点，在圆内的飞标（点）有c点，则依比例来算，就会得到上图中最后的公式，落入圆内的点满足条件x2+y2=1.解法如下：

#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;
}
使用这种方法很巧妙，使用了概率的知识，相信生活中也有很多问题可以用概率的方法来得到结果，只有在以后实践中才知道了！