随机算法之求π

用随机方法求π,具体方法如下:
易想到用 S(正方形内切圆 )与S(正方形）比值求得
我们首先对内切圆圆心为坐标原点建立坐标系, 考虑到随机函数rand()
产生0~RAND_MAX的随机数,所以产生的点均落在第一象限,
第一象限1/4圆和圆外 的点的个数比值即为 S(正方形内切圆 )与S(正方形）比值

?

#include <stdio.h>

#include <stdlib.h> //srand,rand,RAND_MAX

#include <time.h> //time

#include <math.h>

 

int main(){

    //比较10次运算结果

    for(int i = 0 ; i < 10 ; i ++ ) {

        double x,y;

        int n,count;

        srand((unsigned int)time(NULL));//初始化伪随机数生成器

        printf("产生随机数个数(尽量较大):");

        scanf("%d",&n);

        count = 0;

        for(int j = 0 ; j < n ; j ++ ){

            x = rand()/(double)RAND_MAX;//0.0~1.0之间

            y = rand()/(double)RAND_MAX;

            if(sqrt(1-x*x) >= y) //圆内

                count ++;

        }

        printf("pai:%5.8f\n",4.0*count/n);

    }

    return 0;

}

 参照:http://www.acmerblog.com/random-%cf%80-2512.html