C语言pow()函数的计算精度问题

编程计算 a+aa+aaa+…+aa…a（n个a）的值，n和a的值由键盘输入。例如，当n=4,a=2,表示计算2+22+222+2222的值。

程序运行结果示例：

Input a,n:

2,4↙

sum=2468

一开始以为这是个非常简单的问题。C实现代码如下：

 

#include<stdio.h>
#include<math.h>
int main(void)
{
    int a = 0;
    int n = 0;
    int sum = 0;
    int i,j;
    int sum1 = 0;

    printf("Input a,n:\n");
    scanf("%d,%d",&a, &n);

    for(i = 0; i < n; i++)
    {
       for(j = 0; j <= i; j++)
       {           
           sum += a*pow(10,j);           
       }     

    }

   // sum1 = 2*pow(10,2);
    //printf("%d\n",2*pow(10,2));
    printf( "sum=%ld\n", sum);
    return 0;
}

 

程序的输出为2467，比正确值少1。debug发现pow（10，j），当进入循环j=2时，计算结果为99。单独计算pow（10,2）结果为100，非常奇怪。查看函数原型发现pow函数的入口和返回值都是double类型，

即double pow( double, double )，基本明白问题所在。我的奇葩编译器在数据类型转换时候出现了问题。改用其他编译器发现没有问题。不再使用pow库函数，改用递归解决此问题。

#inlcude <stdio.h>
int powCaculate(int a, int n);
int main(void)
{
    int a = 0;
    int n = 0;
    int i = 0;
    int sum=0;

    printf("Input a,n:\n");
    scanf("%d,%d",&a, &n);

    for(i = 0; i < n; i++)
    {
        printf("%d\n",sum += powCaculate(a, i));
    }

    return 0;
}

int powCaculate(int a, int n)
{
    if(n == 0)
        return a;
    else if(n == 1)
        return a*10 + a;
    else
        return powCaculate(a, n-1)*10 + 2;
}