为何C++中pow()函数求区间阿姆斯特朗数结果错误，连乘却正确？

这是典型的浮点运算精度问题导致的结果差异，咱们一步步拆解清楚：

1. pow()的本质是浮点运算 #

pow()是C标准库中用于计算幂的函数，它的返回值是double（双精度浮点型），而非整数。虽然从数学角度看pow(3,3)绝对等于27，但计算机的浮点存储机制会让部分整数幂的结果变成近似值——比如pow(3,3)实际可能被存储为26.999999999999996，而非精确的27。

2. 浮点转整数的截断坑 #

你的代码里sum是int类型，当把pow(temp,3)的结果加到sum时，会自动触发浮点转整数的截断行为：直接丢弃小数部分，只保留整数位。

• 当temp=3时，pow(3,3)的近似值26.999999999999996会被截断成26，而非正确的27。
• 计算153的各位立方和时，原本的1³+5³+3³=1+125+27=153，就变成了1+125+26=152，和原数153不相等，自然就被漏掉了。

而temp*temp*temp是纯整数运算，所有计算结果都是精确的整数，完全不会有精度损失，所以能输出包含153在内的所有正确阿姆斯特朗数。

3. 快速验证这个问题 #

你可以在代码里加一行测试代码，看看pow(3,3)的实际值：

cout << "pow(3,3)的实际值：" << pow(3,3) << endl;

运行后你会发现，输出的大概率不是精确的27，而是一个无限接近27的小数。

解决办法（如果一定要用pow()） #

如果坚持要用pow()函数，记得用round()函数对结果做四舍五入后再转整数，避免截断误差：

sum = sum + round(pow(temp,3));

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