思维之花-方程

题目描述

  有形如：ax³+bx²+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数(a，b，c，d 均为实数)，并约定该方程存在三个不同实根(根的范围在-100至100之间)，且根与根之差的绝对值> =1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后2位。
  提示：记方程f(x)=0，若存在2个数x1和x2，且x1< x2，f(x1)*(x2)< 0，则在(x1，x2)之间一定有一个根。

输入数据

  输入该方程中各项的系数 (a ， b ， c ， d 均为实数)，

输出数据

  由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后 2 位。

样例输入

1 -5 -4 20

样例输出

-2.00 2.00 5.00

程序分析¹

• 这是一道经典的分治问题。分治怎么分，以什么作为端点来分治？首先采用的便是二分策略，左端点x到中间点mid，中间点mid再到右端点y，以此类推。

• rootNum ==3时，就停止运算，消减了很多的运算步骤。

• 还有就是当f（x）==0时，我们也是选择直接输出，多判断了一步，却可能省下很多步，十分的划算

class solveEquation:
    rootNum = 0
    rootList = []

    def f(self, x):
        return a * pow(x, 3) + b * (x ** 2) + c * x + d

    def BiSearch(self, ld, rd):  # 调用这个函数说明(ld，rd)区间内一定有解
        # 有两个递归出口
        if abs(ld - rd) < 0.001:  # 若解没有落在端点上，而是再(ld,rd)之间，在满足精度的条件下，就可以用ld或rd当作答案
            self.rootNum += 1
            # print("%.2f" % ld)
            self.rootList.append(format(ld, ".2f"))
            return
        mid = (ld + rd) / 2
        if self.f(mid) == 0:  # 解落在了端点上
            self.rootNum += 1
            # print("%.2f" % mid)
            self.rootList.append(format(mid, ".2f"))
            return
        if self.f(ld) * self.f(mid) < 0:
            self.BiSearch(ld, mid)
        else:
            self.BiSearch(mid, rd)

a, b, c, d = map(float, input().split())
SE = solveEquation()
for i in range(-100, 100):
    ld, rd = i, i + 1
    if SE.f(ld) == 0:  # 判断左端点是解，通过for循环可以依次遍历所有端点
        # print("%.2f" % ld)
        SE.rootList.append(format(ld, ".2f"))
        SE.rootNum += 1
        continue
    if SE.f(ld) * SE.f(rd) < 0:  # <0表示区间(ld,rd)内一定存在一个解，但是没有考虑到端点可能是解，需要额外判断
        SE.BiSearch(ld, rd)
    if SE.rootNum == 3:
        break
for i in range(2):
    print(SE.rootList[i], end=' ')
print(SE.rootList[2])

---

1. https://blog.csdn.net/qq_43049569/article/details/82084767 ↩︎