#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",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+CgppbnQgbWFpbih2b2lkKSB7CgkKCWRvdWJsZSB4ID0xZTE1OwoJZG91YmxlIHgxID0xZTE2OwogCglkb3VibGUgcmVzdWx0ID0gc3FydCh4ICsgMSkgLSBzcXJ0KHgpOwoJZG91YmxlIHJlc3VsdDEgPSBzcXJ0KHgxICsgMSkgLSBzcXJ0KHgxKTsKCWRvdWJsZSByZXN1bHQyYT0gMSAvIChzcXJ0KHggKyAxKSArIHNxcnQoeCkpOwoJZG91YmxlIHJlc3VsdDJiPSAxIC8gKHNxcnQoeDEgKyAxKSArIHNxcnQoeDEpKTsKCQoJcHJpbnRmKCJzcXJ0KHggKyAxKSAtIHNxcnQoeClcbiIpOwoJcHJpbnRmKCJyZXN1bHQoMTBeMTUpOiUuMjVmXG4iLHJlc3VsdCk7CglwcmludGYoInJlc3VsdCgxMF4xNik6JS4yNWZcbiIscmVzdWx0MSk7CglwcmludGYoIlxuIik7CgkKCXByaW50ZigiMSAvIChzcXJ0KHggKyAxKSArIHNxcnQoeCkpXG4iKTsKCXByaW50ZigicmVzdWx0MigxMF4xNSk6JS4yNWZcbiIscmVzdWx0MmEpOwoJcHJpbnRmKCJyZXN1bHQyKDEwXjE2KTolLjI1ZlxuIixyZXN1bHQyYik7CglyZXR1cm4gMDsKfQo=