Buffon投针（近似计算π）

问题定义：

如果地板的宽度是针的两倍，则向地板投针，与地板缝隙相交的概率为1/π

参考资料：

http://www.hudong.com/wiki/Buffon%E6%8A%95%E9%92%88%E9%97%AE%E9%A2%98

程序如下：

#include<stdio.h> #include<ctime> #include<cstdlib> #include <math.h> // 设地板的宽度为2，针的长度为1 void rand_seed(); float Randomf( float a, float b); int Randomi( int i, int j); int main(void) { float x1, y1, x2, y2; const int MAX= 20000000; int nCase = MAX; int k = 0; // 针与缝隙交叉的次数 double high, width; int degree; while( nCase-- ) { x1 = Randomf(2.0,100.0); y1 = Randomf(2.0,100.0); degree = Randomi( 0, RAND_MAX); high = sin( degree ); width = cos( degree ); x2 = x2 + width; y2 = y1 + high; //确定x2以后，由于两个点之间的距离是恒定的，所以y2也求出来了，但是这里y2可以去两个值，为了公平起见，交替取以过（x1，y1）的水平线以上的点，和水平线以下的点。 int min,max; if( y1 > y2 ) { min = (int)y2; max = (int)y1; } else { min = (int)y1; max = (int)y2; } if( max % 2 == 0 && min % 2 != 0) k++; } printf("%g\n", MAX * 1.0 / k ); return 0; } /** 设置随机数生成的种子 */ void rand_seed() { int seed=static_cast<int>(time(0)); srand(seed); } float Randomf( float a, float b) { int na = (int)(a * 10000); int nb = (int)(b * 10000); int temp = rand() % ( nb - na + 1 ); temp += na; return (temp / 10000.0); } int Randomi( int i, int j) { int temp; temp =( rand() % ( j - i + 1 )) + i; return temp; }