#include <stdio.h>
#include <math.h>
int main(void) {
double x =1e15;
double x1 =1e16;
double result1
= sqrt(x1
+ 1) - sqrt(x1
); double result2a
= 1 / (sqrt(x
+ 1) + sqrt(x
)); double result2b
= 1 / (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\n",result1
);
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+CgppbnQgbWFpbih2b2lkKSB7CgkKCWRvdWJsZSB4ID0xZTE1OwoJZG91YmxlIHgxID0xZTE2OwogCglkb3VibGUgcmVzdWx0ID0gc3FydCh4ICsgMSkgLSBzcXJ0KHgpOwoJZG91YmxlIHJlc3VsdDEgPSBzcXJ0KHgxICsgMSkgLSBzcXJ0KHgxKTsKCWRvdWJsZSByZXN1bHQyYT0gMSAvIChzcXJ0KHggKyAxKSArIHNxcnQoeCkpOwoJZG91YmxlIHJlc3VsdDJiPSAxIC8gKHNxcnQoeDEgKyAxKSArIHNxcnQoeDEpKTsKCQoJcHJpbnRmKCJzcXJ0KHggKyAxKSAtIHNxcnQoeClcbiIpOwoJcHJpbnRmKCJyZXN1bHQoMTBeMTUpOiUuMjVmXG4iLHJlc3VsdCk7CglwcmludGYoInJlc3VsdCgxMF4xNik6JS4yNWZcblxuIixyZXN1bHQxKTsKCQoJcHJpbnRmKCIxIC8gKHNxcnQoeCArIDEpICsgc3FydCh4KSlcbiIpOwoJcHJpbnRmKCJyZXN1bHQyKDEwXjE1KTolLjI1ZlxuIixyZXN1bHQyYSk7CglwcmludGYoInJlc3VsdDIoMTBeMTYpOiUuMjVmXG4iLHJlc3VsdDJiKTsKCXJldHVybiAwOwp9Cg==