HJ107 求解立方根 泰勒+二分 两种解法 C/C++

题目链接

主要思路

方法一：泰勒
假设输入为n，即要求n的立方根，我们设n的立方根为x，即$x^3=n$。设函数$f(x)=x^3-n，令f(x)=0$，将$f(x)$泰勒展开：
$f(x)=f(x_0)+f^{\prime }(x_0)(x-x_0)+...+o(x-x_0)$
代入得：$f(x)=x_0^3-n+3x_0^2(x-x0)=0$
化简得：$x=x_0-(x_0^3-n)/(3x_0^2)$
再化简：$x=2/3\ast x_0+n/3x_0^2$
初始设定x为一个真正答案附近的值，不要太大，也不要太小，然后按照上述求解x的公式，循环逼近，最后在误差允许的范围内求出结果。

---

方法二:二分
初始设定一个区间，答案一定在这个区间内，每次取该区间中点的值进行运算比较，根据与实际值的大小关系，调整左右边界值，不停缩短区间，最终找到答案。

---

#include <iostream>
#include <algorithm>
#include <string>
#include <cmath>
#define inf 0x3f3f3f3f
using namespace std;
const double eps = 1e-6;//精度

//方法一 泰勒
/*
假设输入为n，即要求n的立方根，我们设n的立方根为x，
即x^3=n。
设函数f(x) = x^3-n，令f(x) = 0，将f(x)泰勒展开：
f(x) = f(x0)+f'(x0)(x-x0)+...+o(x-x0)
代入得：f(x) = x0^3-n+3x0^2(x-x0) = 0
化简得：x = x0 - (x0^3-n)/(3x0^2)
再化简：x = 2/3*x0+n/(3*x0*x0);
*/
double cube1(double n) {
    double x = 2.5;//初值x随便取，不要太大，不要太小
    while(fabs(x*x*x-n)>eps) {
        //x = x-(pow(x,3)-n)/(3*x*x);
        x = 2.0/3*x+n/(3*x*x);
    }
    return x;
}

//方法二：二分
double cube2(double n) {
    double low = min(-1.0,n);
    double high = max(1.0,n);
    double x = (low+high)/2;
    while(fabs(x*x*x-n)>eps) {
        if(x*x*x<n)low = x;
        else high = x;
        x = (low+high)/2.0;
    }
    return x;
}
int main(){
    double n;
    while(cin>>n) {
        printf("%.1lf\n",cube1(n));//调用法一
        //printf("%.1lf\n",cube2(n));//调用法二
    }
    return 0;
}