c语言 水仙花数

水仙花数是指一个N位正整数（N≥3），它的每个位上的数字的N次幂之和等于它本身。例如：153=13+53+33。 本题要求编写程序,计算所有N位水仙花数。

输入格式:

输入在一行中给出一个正整数N（3≤N≤7）。

输出格式:

按递增顺序输出所有N位水仙花数，每个数字占一行。

输入样例:

3

输出样例:

153
370
371
407

源代码

#include <stdio.h>
#include <math.h>
int main(){
    int n,i;
    
    scanf("%d",&n);               //输入要判断水仙花数的位数。
    int a=pow(10,n-1),b=pow(10,n);        //定义a,b为n位数的范围用以历遍

    for(i=a;i<b;i++){
        int arr[50]={0};  //不能确定有几位数，所以定义数组用来保存数字个，十，百等位的数值。
        int sum=0,j;                    
        int ti=i;                          //设置变量ti保存i值
        
        for(j=0;j<n;j++){
            arr[j]=ti%10;             //求出各个位上的值
            ti/=10;            
        }
        for( j=0;j<n;j++){
            sum=sum+pow(arr[j],n);
            
        }   
        if(sum==i){                //判断水仙花数
                printf("%d\n",i);
            }
    }

    return 0;
}

结果：

可以看到当N的值大些时，程序运行超时了。

而原因在于这个pow函数！

通过Stack Overflow和CSDN其他博主那我了解到：

实际上pow()函数是一个浮点函数，计算pow(x,y)时 int 在调用 pow 之前会被转化为double，这样算法就复杂多了，这不能像整那样循环或迭代来计算，因为指数有可能是小数。在运算中要调用符点运算的乘法，除法，log。

因此，使用速度比两个整型相乘慢，并且pow是函数，传递参数，调用函数的代价远大于直接相乘。如果指数是整数，则不要使用 pow。所以非必要一般不用pow函数，直接写x*x他不香吗。（还可以少打个调用math库）

因此，我们将pow函数进行修改（这里我直接自定义函数，有更好的方法可以提出帮助我改进）

#include <stdio.h>
int pow(int o,int w);         //自定义pow函数
int main(){                        //一样的流程
    int n,i;
    
    scanf("%d",&n);
    int a=pow(10,n-1),b=pow(10,n);

    for(i=a;i<b;i++){
        int arr[50]={0};
        int sum=0,j;
        int ti=i;
        
        for(j=0;j<n;j++){
            arr[j]=ti%10;
            ti/=10; 
        }
        for( j=0;j<n;j++){
            sum=sum+pow(arr[j],n);
            
        }
        if(sum==i){
                printf("%d\n",i);
            }
        
    }
    return 0;
}
int pow(int o,int w){            // 利用循环达成o的w次方
    int sum=1,i;
    for(i=0;i<w;i++){
        sum=sum*o;
    }
    return sum;
}

结果：