求助：实现通用实数结果根函数，解决std::pow负数开根失效问题

针对你遇到的std::pow在处理负数非整数次根时返回NaN的问题，这里提供两种实现方案，分别对应严格数学定义和满足你示例需求的扩展逻辑：

方案一：严格遵循实数域数学定义 #

该方案仅在实数域存在有效根时返回结果，否则返回NaN：

#include <cmath>
#include <limits>

double myroot(double n, double k)
{
    // 0次方根无定义，返回NaN
    if (k == 0.0) {
        return std::nan("");
    }

    // 处理0的次方根
    if (n == 0.0) {
        return k > 0.0 ? 0.0 : std::numeric_limits<double>::infinity();
    }

    // 正数的任何非零次方根都存在唯一正实根
    if (n > 0.0) {
        return std::pow(n, 1.0 / k);
    }

    // 处理负数的情况：仅当k是奇数整数时存在实根
    const double k_int = std::round(k);
    // 判断k是否为整数（考虑浮点数精度误差）
    if (std::fabs(k - k_int) > 1e-9) {
        return std::nan("");
    }

    // 判断是否为偶数整数
    if (static_cast<long long>(k_int) % 2 == 0) {
        return std::nan("");
    }

    // 奇数整数次根，返回负的绝对值开方结果
    return -std::pow(-n, 1.0 / k);
}

逻辑说明：

• 负数仅在开奇数整数次根时存在实数解，其他情况（如非整数次根、偶数整数次根）均返回NaN，完全符合实数域的数学定义。
• 使用1e-9作为浮点数整数判断的阈值，避免因精度误差导致的误判。

方案二：满足你示例需求的扩展实现 #

该方案扩展了判断逻辑，允许负数开非整数次根时返回近似结果（即-|n|^(1/k)），仅在开偶数整数次根时返回NaN：

#include <cmath>
#include <limits>

double myroot(double n, double k)
{
    // 0次方根无定义，返回NaN
    if (k == 0.0) {
        return std::nan("");
    }

    // 处理0的次方根
    if (n == 0.0) {
        return k > 0.0 ? 0.0 : std::numeric_limits<double>::infinity();
    }

    // 正数的任何非零次方根都存在唯一正实根
    if (n > 0.0) {
        return std::pow(n, 1.0 / k);
    }

    // 处理负数的情况：仅排除偶数整数次根
    const double k_int = std::round(k);
    if (std::fabs(k - k_int) <= 1e-9) {
        // k是整数，判断是否为偶数
        if (static_cast<long long>(k_int) % 2 == 0) {
            return std::nan("");
        }
    }

    // 非偶数整数的情况，返回负的绝对值开方结果
    return -std::pow(-n, 1.0 / k);
}

逻辑说明：

• 对于myroot(-125,3.001)这类场景，会返回-125^(1/3.001)的近似值，符合你给出的预期结果。
• 严格排除负数开偶数整数次根的情况（如myroot(-4,2)返回NaN），避免明显违反数学常识的结果。

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