#include <stdio.h>
#include <math.h>
#define PARTS 740000000
int main(void) {
double old_y,new_y,old_x,new_x,d,S;
int i;
d = 1.0 / PARTS;
S = 0;
old_y = 1;
old_x = 0;
for (i = 1;i <= PARTS;i++) {
new_x = i * d;
new_y
= sqrt(1 - new_x
* new_x
);
S += ((old_y + new_y) * d / 2);
old_y = new_y;
old_x = new_x;
}
printf ("分割数 = %d,π = %0.9f\n",PARTS
,S
* 4);
return 0;
}
I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxtYXRoLmg+CiNkZWZpbmUgUEFSVFMgNzQwMDAwMDAwCiAKaW50IG1haW4odm9pZCkgewpkb3VibGUgb2xkX3ksbmV3X3ksb2xkX3gsbmV3X3gsZCxTOwppbnQgaTsKIApkID0gMS4wIC8gUEFSVFM7ClMgPSAwOwpvbGRfeSA9IDE7Cm9sZF94ID0gMDsKIApmb3IgKGkgPSAxO2kgPD0gUEFSVFM7aSsrKSB7Cm5ld194ID0gaSAqIGQ7Cm5ld195ID0gc3FydCgxIC0gbmV3X3ggKiBuZXdfeCk7CiAKUyArPSAoKG9sZF95ICsgbmV3X3kpICogZCAvIDIpOwogCm9sZF95ID0gbmV3X3k7Cm9sZF94ID0gbmV3X3g7Cn0KIApwcmludGYgKCLliIblibLmlbAgPSAlZCzPgCA9ICUwLjlmXG4iLFBBUlRTLFMgKiA0KTsKIApyZXR1cm4gMDsKfQo=