Python如何实现分数立方根计算 返回分母有理化的分数结果

需求说明 #

开发任意三次方程求解程序时，需要输出根的精确非十进制形式，核心需要实现支持分数输入、返回分母有理化的精确立方根结果的函数，替代原有返回浮点十进制值的立方根逻辑。

原有浮点立方根实现 #

原有代码仅返回十进制近似值，无法满足精确输出要求：

third = 1./3. #Cube root
def cubic_root(x):
        if x >= 0:
            return x**third
        else:
            return -(-x)**third

原有尝试代码的典型问题 #

之前编写的版本存在几类影响运行和结果正确性的错误：

• 变量名拼写错误：比如负数分支里把cubeRoot_num写成cubedRoot_num、cubeRoot写成root，会直接触发运行报错
• 负数立方根逻辑错误：实数范围内负数立方根是实数，不需要加虚数单位i，原代码错误给负数结果加了i
• 浮点数精度问题：用int(x**(1/3))取整时，因为浮点计算误差，比如8**(1/3)可能得到1.999999999999，取整后变成1，漏判完全立方数
• 分母有理化缺失：分开处理分子分母后没有做分母去根号的有理化操作，不符合精确表达式的书写规范
• 重复冗余代码：整数、分子、分母的化简逻辑完全重复，没有封装成通用函数，维护成本高

可直接集成的实现方案 #

先封装两个通用工具函数，再实现最终的分数立方根函数，返回标准化的精确表达式字符串：

import math
from fractions import Fraction

RADICAL = "³√"

def get_prime_factors(n):
    """返回正整数n的质因数分解字典，key为质因数，value为指数"""
    factors = {}
    while n % 2 == 0:
        factors[2] = factors.get(2, 0) + 1
        n = n // 2
    i = 3
    while i*i <= n:
        while n % i == 0:
            factors[i] = factors.get(i, 0) + 1
            n = n // i
        i += 2
    if n > 2:
        factors[n] = 1
    return factors

def simplify_radical(n):
    """
    化简正整数的三次根号，返回(外系数, 内剩余数)
    比如simplify_radical(24)返回(2, 3)，对应2³√3
    """
    if n == 1:
        return (1, 1)
    factors = get_prime_factors(n)
    outer = 1
    inner = 1
    for prime, exp in factors.items():
        outer *= prime ** (exp // 3)
        inner *= prime ** (exp % 3)
    return (outer, inner)

def format_cubic_root(num):
    """
    对分数/整数取三次根，返回分母有理化后的精确表达式字符串
    支持输入整数、浮点数（会自动转最简分数）、Fraction分数对象
    """
    # 统一转成最简分数处理
    if not isinstance(num, Fraction):
        num = Fraction(num).limit_denominator()
    sign = 1
    numerator = num.numerator
    denominator = num.denominator
    # 处理负号，全部放到结果最外层
    if numerator < 0:
        sign = -1
        numerator = -numerator
    
    # 分别化简分子分母的三次根
    num_outer, num_inner = simplify_radical(numerator)
    den_outer, den_inner = simplify_radical(denominator)

    # 分母有理化：给分子分母同乘den_inner²，把分母根号完全去掉
    if den_inner != 1:
        new_num_inner = num_inner * den_inner * den_inner
        new_num_outer, new_num_inner = simplify_radical(new_num_inner)
        final_num_outer = num_outer * new_num_outer
        final_num_inner = new_num_inner
        final_den = den_outer * den_inner
    else:
        final_num_outer = num_outer
        final_num_inner = num_inner
        final_den = den_outer
    
    # 拼接符合数学书写规范的结果字符串
    res = ""
    if sign == -1:
        res += "-"
    # 处理分子部分
    if final_num_inner == 1:
        num_part = str(final_num_outer)
    else:
        if final_num_outer == 1:
            num_part = f"{RADICAL}{final_num_inner}"
        else:
            num_part = f"{final_num_outer}{RADICAL}{final_num_inner}"
    # 处理分母部分，分母为1时直接返回分子结果
    if final_den == 1:
        res += num_part
    else:
        res += f"{num_part}/{final_den}"
    return res

使用示例 #

# 测试用例
print(format_cubic_root(8))      # 输出 2
print(format_cubic_root(-1/8))   # 输出 -1/2
print(format_cubic_root(2/27))   # 输出 ³√2/3
print(format_cubic_root(16/9))   # 输出 2³√6/3 （分母有理化后的标准结果）
print(format_cubic_root(-24))    # 输出 -2³√3

实现说明 #

• 用质因数分解代替浮点开三次方取整，完全规避浮点数精度误差，不会漏判完全立方数
• 自动统一处理负号，实数范围内负数立方根直接输出负的实数结果，不会错误引入虚数单位
• 自动完成分母有理化，输出的表达式符合常规数学书写规范
• 所有化简逻辑封装成通用函数，没有重复代码，后续可以直接传入三次方程求解过程中得到的分数中间值调用

如果三次方程求解过程中出现负数开平方的情况（即判别式小于0的不可约情形），再单独引入虚数单位i处理即可，立方根本身的实数计算不需要加i。

内容的提问来源于stack exchange，提问作者Patman1O1