[C/C++ 演算法]- 蒙地卡羅法求 PI

[C/C++ 演算法]- 蒙地卡羅法求 PI

剛才找資料時發現一個C/C++的教學網站，趕快發揮(C/P)的長才將它備份來，有需要的同好，歡迎來(C/P)一下^^。

拷貝來源：
http://openhome.cc/Gossip/AlgorithmGossip/
http://openhome.cc/Gossip/AlgorithmGossip/MathPI.htm

#include <stdio.h> #include <stdlib.h> #include <time.h> #define N 50001 int main(void) { srand(time(NULL)); int sum = 0; int i; for(i = 1; i < N; i++) { double x = (double) rand() / RAND_MAX; double y = (double) rand() / RAND_MAX; if((x * x + y * y) < 1) { sum++; } } printf("PI = %f\n", (double) 4 * sum / (N - 1)); return 0; }