题解 | #求解立方根#

求解立方根

http://www.nowcoder.com/practice/caf35ae421194a1090c22fe223357dca

# 牛顿迭代法求解立方根的思路:
# 令f(x) = x^3 - a,求解f(x) = x^3 - a = 0。
# 利用泰勒公式展开,即f(x)在x0处的函数值为:
# f(x) = f(x0) +f'(x0)(x-x0) = (x0^3-a) + (3x0^2)(x-x0) = 0,
# 解之得:x = x0 - (x0^3 - a) / (3x0^2)。
#     即 x = x - ((x*x*x - n) / (3*x*x));

# 拓展:求平方根用一个套路:
# 令f(x) = x^2 - a,求解f(x) = x^2 - a = 0。
# 利用泰勒公式展开,即f(x)在x0处的函数值为:
# f(x) = f(x0) +f'(x0)(x-x0) = (x0^2-a) + 2x0(x-x0) = 0,
# 解之得:x = x0 - (x0^2 - a) / 2x0
#     即 x = x - (x*x-a)/2x 可进一步化简为:=(x+a/x) / 2。

# 总结:
# 平方根与立方根的求解迭代公式:
# 新x = 旧x - f(x)/f'(x)
# 新x = 旧x - (x平方或者立方与输入数a的差)/f(x)求导数
'''
# 法一:牛顿迭代法
while True:
    try:
        a = float(input().strip())  # 获取输入的实数a
        e = 0.0001  # 设定一个精度值
        t = a  # 初始化立方根t的值为输入的值a
        while abs(t*t*t - a) > e:  # 差值没有达到精度,便一直更新立方根
    # x(i+1) = x(i) - f(xi)/f'(xi)
    # 更新后的x = 原x - (原x的立方-a)/f(原x)导数
            t = t - (t*t*t - a) * 1.0 / (3 * t*t)
        print("%.1f" % t)  # 当精度达到要求时,此时的立方根t便为输入实数的立方根解
    except:
        break
'''

# 法二:二分法

while True:
    try:
        a = float(input().strip())
        epsilon = 0.0001
        low = min(-1.0, a)
        high = max(1.0, a)
        ans = (low + high)/2
        while abs(ans**3 - a) >= epsilon:
            if ans**3 < a:
                low = ans
            else:
                high = ans
            ans = (low + high)/2.0
        print('%.1f' % ans)
    except:
        break

	


全部评论
#41 42 行 ## 同时处理了abs<1 以及负数的情况 ## 注意minmax前后两个数字要同类型
点赞 回复 分享
发布于 2022-07-29 22:27
这题我没看懂,直接对输入x求它的1/3次幂不就行了吗
2 回复 分享
发布于 2022-08-07 10:02
佬讲的太清楚了,一下子就学会了两种方法。
1 回复 分享
发布于 2024-10-03 14:49 新疆
看不懂,并不妨碍我膜拜大神
1 回复 分享
发布于 2022-04-25 22:49
何不偷鸡?
点赞 回复 分享
发布于 2024-09-12 20:48 陕西
为何不直接 x ** (1/3)? 又没用库
点赞 回复 分享
发布于 2024-08-25 19:11 上海
二分法其实也使用了三次方这个函数,所以二分法个人认为也不行。
点赞 回复 分享
发布于 2023-04-15 21:11 上海
泰勒展开 熊弟们高数没学过嘛
点赞 回复 分享
发布于 2023-03-04 16:00 湖北
没有考虑输入的值为负数或者输入的值小于1的情况,不够严谨...
点赞 回复 分享
发布于 2023-02-06 12:29 广东
厉害,根本想不到
点赞 回复 分享
发布于 2022-06-05 10:56
点赞 回复 分享
发布于 2022-04-03 16:15

相关推荐

03-06 18:20
门头沟学院 Java
点赞 评论 收藏
分享
评论
87
28
分享

创作者周榜

更多
牛客网
牛客网在线编程
牛客网题解
牛客企业服务