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