#include <stdio.h>
#include <math.h>
int main(void) {
double x =1e15;
double x1 =1e16;
double result1
= sqrt(x1
+ 1) - sqrt(x1
);
printf("sqrt(x + 1) - sqrt(x)\n"); printf("result(10^15):%.25f\n",result
); printf("result(10^16):%.25f\n",result1
);
double x2 =1e15;
double x3 =1e16;
double result2a
= 1 / (sqrt(x2
+ 1) + sqrt(x2
)); double result2b
= 1 / (sqrt(x3
+ 1) + sqrt(x3
));
printf("1 / (sqrt(x + 1) + sqrt(x))\n"); printf("result2(10^15):%.25f\n",result2a
); printf("result2(10^16):%.25f\n",result2b
); return 0;
}
I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxtYXRoLmg+CgppbnQgbWFpbih2b2lkKSB7CgkKCWRvdWJsZSB4ID0xZTE1OwoJZG91YmxlIHgxID0xZTE2OwogCglkb3VibGUgcmVzdWx0ID0gc3FydCh4ICsgMSkgLSBzcXJ0KHgpOwoJZG91YmxlIHJlc3VsdDEgPSBzcXJ0KHgxICsgMSkgLSBzcXJ0KHgxKTsKCQoJcHJpbnRmKCJzcXJ0KHggKyAxKSAtIHNxcnQoeClcbiIpOwoJcHJpbnRmKCJyZXN1bHQoMTBeMTUpOiUuMjVmXG4iLHJlc3VsdCk7CglwcmludGYoInJlc3VsdCgxMF4xNik6JS4yNWZcbiIscmVzdWx0MSk7CglwcmludGYoIlxuIik7CgkKCWRvdWJsZSB4MiA9MWUxNTsKCWRvdWJsZSB4MyA9MWUxNjsKCQoJZG91YmxlIHJlc3VsdDJhPSAxIC8gKHNxcnQoeDIgKyAxKSArIHNxcnQoeDIpKTsKCWRvdWJsZSByZXN1bHQyYj0gMSAvIChzcXJ0KHgzICsgMSkgKyBzcXJ0KHgzKSk7CgkKCXByaW50ZigiMSAvIChzcXJ0KHggKyAxKSArIHNxcnQoeCkpXG4iKTsKCXByaW50ZigicmVzdWx0MigxMF4xNSk6JS4yNWZcbiIscmVzdWx0MmEpOwoJcHJpbnRmKCJyZXN1bHQyKDEwXjE2KTolLjI1ZlxuIixyZXN1bHQyYik7CglyZXR1cm4gMDsKfQo=