#include <stdio.h>
#include <stdlib.h>
double estimate_pi(long long n){
double randomx, randomy, equation, pi;
long long i, incircle = 0;
for(i = 0; i < n; i++){
randomx = (double)(rand() % (1+1-0) + 0);
randomy = (double)(rand() % (1+1-0) + 0);
printf("%d %d\n", randomx, randomy)
equation = randomx * randomx + randomy * randomy;
if(equation <= 1){
incircle++;
}
}
pi = (long double)4 * (long double)incircle / (long double)n;
return pi;
}
int main()
{
printf("%d\n", estimate_pi(100));
return 0;
}
I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxzdGRsaWIuaD4KCgpkb3VibGUgZXN0aW1hdGVfcGkobG9uZyBsb25nIG4pewogICAgZG91YmxlIHJhbmRvbXgsIHJhbmRvbXksIGVxdWF0aW9uLCBwaTsKICAgIGxvbmcgbG9uZyBpLCBpbmNpcmNsZSA9IDA7CgogICAgZm9yKGkgPSAwOyBpIDwgbjsgaSsrKXsKICAgICAgICByYW5kb214ID0gKGRvdWJsZSkocmFuZCgpICUgKDErMS0wKSArIDApOwogICAgICAgIHJhbmRvbXkgPSAoZG91YmxlKShyYW5kKCkgJSAoMSsxLTApICsgMCk7CgoJCXByaW50ZigiJWQgJWRcbiIsIHJhbmRvbXgsIHJhbmRvbXkpCgogICAgICAgIGVxdWF0aW9uID0gcmFuZG9teCAqIHJhbmRvbXggKyByYW5kb215ICogcmFuZG9teTsKCiAgICAgICAgaWYoZXF1YXRpb24gPD0gMSl7CiAgICAgICAgICAgIGluY2lyY2xlKys7CiAgICAgICAgfQogICAgfQoKICAgIHBpID0gKGxvbmcgZG91YmxlKTQgKiAobG9uZyBkb3VibGUpaW5jaXJjbGUgLyAobG9uZyBkb3VibGUpbjsKCiAgICByZXR1cm4gcGk7Cn0KCmludCBtYWluKCkKewoJcHJpbnRmKCIlZFxuIiwgZXN0aW1hdGVfcGkoMTAwKSk7CglyZXR1cm4gMDsKCn0=