pow 函数误差分析及解决方案

我们都知道 （int）99.99999 的结果是99，取其整数部分

对于幂函数pow，是用数值逼近的方法给出的，并不保证严格精确，执行效率也不高，在不同平台精确度可能有一定的差异，所以也就会导致结果不是我们预期想要的。

分析:

#include<stdio.h>
#include<math.h>
int main(){
    int c=4;
    int w=1234;
    int s;
    for(int i=0;i<c;i++){
        s=w%(int)(pow(10,c-i)+0.5)/(int)(pow(10,c-i-1)+0.5);
        printf("%d ",s);
    }
    return 0;
}

这段代码目的是想输出数字 1234 的每一位，一次输出一位。
本来想着的结果会是 1 2 3 4
然而实际的结果却是 1 2 4 4
调整其他数字，发现每次都是倒数第二位数字出错
用 codeblocks 单步调试，追踪变量，也没发现为什么

---

解决方案：

s=w%(int)(pow(10,c-i))/(int)(pow(10,c-i-1));

改为

s=w%(int)(pow(10,c-i)+0.5)/(int)(pow(10,c-i-1)+0.5);

其实其实质是进行四舍五入，然后再取整数部分，因为该函数有一定的准确性，误差不可能太大。

最后应该注意：如果按照四舍五入的方法来处理，最后必须对结果进行强制类型转换，不然当结果乘上一个比较大的数时，在通常情况下结果必定会出错！