#include <iostream>
#include <cstdlib>
#include <ctime>
using namespace std;
int main() {
srand(time(0));
int N = 1000;
int count = 0;
for(int i = 0; i < N; i++) {
double x = (double)rand()/RAND_MAX;
double y = (double)rand()/RAND_MAX;
if(x*x + y*y <= 1) count++;
}
double pi = 4.0 * count / N;
cout << "Estimated Pi = " << pi << endl;
return 0;
}
I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y3N0ZGxpYj4KI2luY2x1ZGUgPGN0aW1lPgp1c2luZyBuYW1lc3BhY2Ugc3RkOwppbnQgbWFpbigpIHsKICAgIHNyYW5kKHRpbWUoMCkpOwogICAgaW50IE4gPSAxMDAwOwogICAgaW50IGNvdW50ID0gMDsKCiAgICBmb3IoaW50IGkgPSAwOyBpIDwgTjsgaSsrKSB7CiAgICAgICAgZG91YmxlIHggPSAoZG91YmxlKXJhbmQoKS9SQU5EX01BWDsKICAgICAgICBkb3VibGUgeSA9IChkb3VibGUpcmFuZCgpL1JBTkRfTUFYOwogICAgICAgIGlmKHgqeCArIHkqeSA8PSAxKSBjb3VudCsrOwogICAgfQogICAgZG91YmxlIHBpID0gNC4wICogY291bnQgLyBOOwogICAgY291dCA8PCAiRXN0aW1hdGVkIFBpID0gIiA8PCBwaSA8PCBlbmRsOwogICAgcmV0dXJuIDA7Cn0K