使用塔尔塔利亚法求解三次方程时的Double自动舍入问题

嘿，我之前在处理数值计算的输出时也踩过一模一样的坑！你的问题根本不是double或float类型的锅，而是printf的格式控制符设置限制了小数位数，导致数值被强制舍入了。

问题根源分析 #

你看到1.54686889被舍成1.50，大概率是因为打印虚部的代码里用了类似%.2f或者%1.2lf的格式符——这种写法会强制只保留两位小数，自然会把第三位的4舍去，变成1.50。而double的精度完全足够存储1.54686889这个数值（double能精确到15-17位有效数字），换成float反而会因为精度更低（仅6-7位有效数字）可能引入额外误差，完全没必要。

具体解决方案 #

修改你打印虚部时的printf格式符，根据需要选择更高精度的格式：

• 如果想完整显示到8位小数，用%.8lf（针对double类型）或%.8f（printf中%f和%lf对double是等价的）
• 如果想自动去掉末尾无用的零，用%g格式符
• 要是需要更多有效数字，比如10位，就用%.10lf

举个简单的验证例子：

#include <stdio.h>

int main() {
    double imaginary_part = 1.54686889;
    printf("错误输出（%.2lf）：%.2lf\n", imaginary_part, imaginary_part);
    printf("正确输出（%.8lf）：%.8lf\n", imaginary_part, imaginary_part);
    printf("自动精简输出（%g）：%g\n", imaginary_part, imaginary_part);
    return 0;
}

运行这段代码会输出：

错误输出（%.2lf）：1.50
正确输出（%.8lf）：1.54686889
自动精简输出（%g）：1.54687

针对你的塔尔塔利亚法程序的建议 #

1. 把类型改回double，它比float更适合三次方程的数值计算（塔尔塔利亚法涉及开立方、三角函数等运算，更高精度能减少中间计算误差）
2. 找到你代码中打印虚部的语句，替换格式符为更高精度的版本

你的代码开头部分整理如下：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <locale.h>
#define PI 3.14159265359

int main(void) {
    setlocale(LC_ALL,"");
    double a, b, c, d, p,...

调整格式符后，就能正确显示完整的虚部数值啦！

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