0.简介
根据上一篇的内容,今天又看到了比较的好的一个优化方法,在这里记录一下。
1.判断条件优化
之前的二分法和牛顿迭代法的判断条件都有乘法,今天发现一个比较的好的判断方法。
二分法
long double l = 0, r = (n>=1 ? n:1), mid=0;
while (fabs(l-r)>1e-8)
{
mid = (l + r) / 2;
if (mid * mid > n) r = mid;
else l = mid;
}
牛顿迭代法
long double x = n,y = -1,k,b,lx = n/2;
while (fabs(lx - x)>1e-6)//计算两次x的相近程度
{
lx = x;//保存x
y = x * x - n;
k = 2 * x;
b = y - k * x;
x = -(b / k);//新的x
}
这样就没有了乘法在判断条件中。
牛顿迭代法的代码再简单整理下。
long double x = n,y = -1,k,b,lx = n/2;
while (fabs(lx - (x = ((lx = x) * x + n) / (2 * x))) > 1e-6);