#include <iostream>
#include <cstdlib>
#include <ctime>
using namespace std;
int main() {
srand(time(0));
int N = 100000, 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;
}
I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y3N0ZGxpYj4KI2luY2x1ZGUgPGN0aW1lPgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKaW50IG1haW4oKSB7CgkKCXNyYW5kKHRpbWUoMCkpOwoJaW50IE4gPSAxMDAwMDAsIE0gPSAwOwoJCglmb3IgKGludCBpID0gMTsgaSA8PSBOOyBpKyspIHsKCQkKCQlkb3VibGUgeCA9IChkb3VibGUpIHJhbmQoKSAvIFJBTkRfTUFYOwoJCWRvdWJsZSB5ID0gKGRvdWJsZSkgcmFuZCgpIC8gUkFORF9NQVg7CgkJCgkJaWYgKHggKiB4ICsgeSAqIHkgPCAxKSBNKys7Cgl9CgkKCWNvdXQgPDwgInBpIH4gIiA8PCA0LjAgKiBNIC8gTiA8PCBlbmRsOwoJCglyZXR1cm4gMDsKfQ==