优化版勾股数暴力算法在b=969时异常的技术求助

核心问题：浮点数精度坑了你的判断 #

你遇到的异常完全是浮点数精度丢失导致的：

• C语言里的double类型只有53位有效二进制精度，对应十进制约15-17位。当a、b的数值增大到一定程度，pow(a,2)计算整数平方时，会因为double无法精确表示大整数的平方产生误差。
• 你用sqrt(pow(a,2)+pow(b,2))算出c后，靠ceil(c) == c判断是否为整数——这个逻辑完全依赖浮点数的精确性，一旦计算过程中出现舍入误差，就会把非整数的c误判成整数，导致所有数对都被当成勾股数。

为什么刚好是969？ #

969只是触发误差的一个临界点：当b=969时，pow(b,2)的计算误差开始明显叠加，和pow(a,2)的误差凑在一起，让平方和的平方根刚好被舍入成整数，触发了大规模误判。至于“勾股数数量是上限的0.969倍”纯属巧合，本质是误差积累到一定程度后，大量非勾股数组合被错误识别。

另外提一句：你的旧算法虽然没出明显问题，但用pow做整数运算本身也有精度隐患，只是三重循环的遍历逻辑刚好抵消了部分误差而已，而且效率极低。

解决办法：彻底抛弃浮点数，用整数运算 #

勾股数的判断完全可以用整数实现，从根源上避免精度问题：

1. 核心逻辑修改 #

把浮点数运算全部换成整数计算，用更大的整数类型避免溢出：

for(unsigned int b = 1; b < limit; b++) {
    for(unsigned int a = 1; a <= b; a++) {
        // 用unsigned long long存平方和，避免溢出
        unsigned long long a_sq = (unsigned long long)a * a;
        unsigned long long b_sq = (unsigned long long)b * b;
        unsigned long long sum_sq = a_sq + b_sq;
        
        // 先算整数c，再验证平方是否等于sum_sq
        unsigned int c = sqrtl(sum_sq); // 用long double版sqrt提高精度
        if((unsigned long long)c * c == sum_sq && c < limit) {
            triples[i][0] = a;
            triples[i][1] = b;
            triples[i][2] = c;
                
            i++;
            printf("found: %d, %d, %d\n", a, b, c);
        }
    }
}

2. 其他优化点 #

• 别用pow算平方：pow是给浮点数设计的，算整数平方既慢又容易出问题，直接用a*a就好。
• 避免栈溢出：原代码里unsigned int triples[limit][3]如果limit很大，会直接把栈撑爆，改成动态内存分配：

  unsigned int (*triples)[3] = malloc(limit * sizeof(*triples));
  // 用完记得释放内存
  free(triples);
  

• 旧算法直接淘汰：三重循环效率太低，而且同样有pow的精度问题，直接用上面的整数版逻辑就行。

效果验证 #

修改后不管是小范围还是大范围测试，都不会再出现误判，而且运行速度比旧算法快N倍（少了一层循环）。

内容的提问来源于stack exchange，提问作者a human being