二分法、牛顿迭代法求sqrt()

最新推荐文章于 2022-07-12 21:20:18 发布

starbug

最新推荐文章于 2022-07-12 21:20:18 发布

阅读量1k

点赞数

分类专栏：数据结构与算法

本文链接：https://blog.csdn.net/u010590166/article/details/17120199

版权

数据结构与算法专栏收录该内容

95 篇文章 0 订阅

订阅专栏

给定一个正数a，不用库函数求其平方根。

设其平方根为x，则有x²=a，即x²-a=0。设函数f(x)= x²-a,则可得图示红色的函数曲线。在曲线上任取一点(x₀,f(x₀))，其中x₀≠0那么曲线上该点的切线方程为

(1-1)

求该切线与x轴的交点得

（1-2）

因为1-2式中x₀作为分母，所以在之前限定了一下初始值不要选0。那么得到的这个与x轴的交点其实是最终要求得的x的一次逼近，我们再以这个x基准继续迭代就可以求得更逼近的x，至于逼近到什么时候才算完，这个取决于你自己设定的精度。整个过程的迭代只需要几步就可以求得最终的结果。

double nt_sqrt(double n)
{
	double x = 1.0;
	while(abs(x*x-n) > 1e-5)
	    x = (x+n/x)/2;
	return x;
}

http://blog.csdn.net/wumuzi520/article/details/7026808

惊人的代码（求1/sqrt())http://www.cnblogs.com/pkuoliver/archive/2010/10/06/1844725.html

二分法：

int sqrt(int x)
{
    int left = 1, right = x / 2;
    int mid;
    int last_mid; //记录最近一次mid
    if (x < 2) return x;
    while(left <= right)
    {
        mid = left + (right - left) / 2;
        if(x / mid > mid)   // 不要用x > mid * mid，会溢出
        {
            left = mid + 1;
            last_mid = mid;
        }
        else if(x / mid < mid)
        {
            right = mid - 1;
        }
        else
        {
            return mid;
        }
    }
    return last_mid;
}

double sqrt(double n) //用二分法
{
    if(abs(n-(1e-5))<0) //小于0的按照你需要的处理
        return n;
    double mid,last;
    double low,up;
    low=0,up=n;
    mid=(low+up)/2;
    do
    {
        if(mid*mid>n)
            up=mid;
        else
            low=mid;
        last=mid;
        mid=(up+low)/2;
    }
    while(abs(mid-last) > eps); //精度控制
    return mid;
}