#include <stdio.h>
#include <math.h>
#include <time.h>
int main() {
const int n = 100000;
const double b = 1.0;
const double a = 0.70710678118;
const double h = (b - a) / n;
double sqrx, x, s;
clock_t start, end;
double total_time = 0.0;
for (int j = 0; j < 1000; j++) {
s = 0.0;
for (int i = 1; i < n; i ++) {
x = a + h * i;
sqrx = x * x;
s
+= sqrt(1 - sqrx
) / sqrx
; }
total_time += (double)(end - start) / CLOCKS_PER_SEC;
}
printf("Average time: %lf seconds\n", total_time
/ 1000.0);
}
I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxtYXRoLmg+CiNpbmNsdWRlIDx0aW1lLmg+CgppbnQgbWFpbigpIHsKICAgIGNvbnN0IGludCBuID0gMTAwMDAwOwogICAgY29uc3QgZG91YmxlIGIgPSAxLjA7CiAgICBjb25zdCBkb3VibGUgYSA9IDAuNzA3MTA2NzgxMTg7CiAgICBjb25zdCBkb3VibGUgaCA9IChiIC0gYSkgLyBuOwogICAgZG91YmxlIHNxcngsIHgsIHM7CiAgICAKICAgIGNsb2NrX3Qgc3RhcnQsIGVuZDsKICAgIGRvdWJsZSB0b3RhbF90aW1lID0gMC4wOwoKICAgIGZvciAoaW50IGogPSAwOyBqIDwgMTAwMDsgaisrKSB7CiAgICAJcyA9IDAuMDsKCSAgICBzdGFydCA9IGNsb2NrKCk7CgkgICAgCgkgICAgZm9yIChpbnQgaSA9IDE7IGkgPCBuOyBpICsrKSB7CgkgICAgICAgIHggPSBhICsgaCAqIGk7CgkgICAgICAgIHNxcnggPSB4ICogeDsKCSAgICAgICAgcyArPSBzcXJ0KDEgLSBzcXJ4KSAvIHNxcng7CgkgICAgfQogICAgCiAgICAJZW5kID0gY2xvY2soKTsKICAgICAgICB0b3RhbF90aW1lICs9IChkb3VibGUpKGVuZCAtIHN0YXJ0KSAvIENMT0NLU19QRVJfU0VDOwogICAgfQoKICAgIHByaW50ZigiQXZlcmFnZSB0aW1lOiAlbGYgc2Vjb25kc1xuIiwgdG90YWxfdGltZSAvIDEwMDAuMCk7CgogICAgcHJpbnRmKCJzID0gJWYiLCBzICogaCk7Cn0=