求解10位数以内的水仙花数：普通的穷举算法 另：优化的深度优先搜索算法。...

/*
水仙花数是指一个 n 位数 ( n≥3 )，它的每个位上的数字的 n 次幂之和等于它本身。
（例如：1^3 + 5^3 + 3^3 = 153）
普通的穷举算法，
CPU0: Intel(R) Pentium(R)
CPU P6200  @ 2.13GHz
求解8位数以内的水仙花数需要45s左右的时间，

*/

/**
@author : 东海陈光剑 2013.4.26 Friday ,01：02

*/

#include <stdio.h>
#include <stdlib.h>
#include <time.h>

//#include <math.h>

/* 求1-N 以内的水仙花数 */
#define N 1e8
//typedef enum {true = 0, false = !true} bool;

/**
打印日志时间
*/
void print_time()
{
    time_t my_time; // time_t
/*
    时间类型（time.h 定义）
struct tm -- 时间结构，time.h 定义如下：
int tm_sec;
int tm_min;
int tm_hour;
int tm_mday;
int tm_mon;
int tm_year;
int tm_wday;
int tm_yday;
int tm_isdst;

*/
    struct tm* timeinfo;
    time( &my_time );//获取时间，以秒计，从1970年1月一日起算，存于my_time
    timeinfo = localtime( &my_time );//转为当地时间，tm 时间结构
    //printf ( "\007The current date is: %s", asctime (timeinfo) );
/*asctime()-- 转为标准ASCII时间格式：星期 月 日 时：分：秒 年*/

    printf ( "[%4d-%02d-%02d %02d:%02d:%02d]",
            1900+timeinfo->tm_year,
            1+timeinfo->tm_mon,
            timeinfo->tm_mday,
            timeinfo->tm_hour,
            timeinfo->tm_min,
            timeinfo->tm_sec);
/*
打印tm，tm_year 从1900年计算，所以要加1900，
月tm_mon，从0计算，所以要加1 */

}


/**
取数字位数
*/
int get_power_index(int n)
{

    int power_index = 0;
    while (n > 0)
    {
        power_index ++;
        n /= 10;
    }

    return power_index;
}

/**计算数字各位上数字的k次方（k=数字的位数）*/

int my_pow(int* a, int* b)
{
    int i=0;
    int p=1;
    //int pow = 1;
    for (i = 0; i < *b; i++)
    {
      p *= *a;
    }

    return p;
}

/**判断x是否是水仙花数*/
int is_flower(int x)
{
    int sum = 0;
    int a = 0;
    //int* pindex =NULL;
    int pindex = get_power_index( x );
   // int bits = power_index;
    int n=x;

    while (n > 0)
    {
        a = n % 10;
        sum += my_pow( & a , & pindex);
        n = (( n - n % 10 )/10);// n去尾
    }


    if (x == sum )
    {
       // printf(" %d 位数: "，*pindex);
        return 0;
    }
    else
    {
        return 1;
    }
}



int main()
{
    // test
    print_time();
    printf("Testing....\nThe result is 9,625,010001000\n");
    //printf("%c\n",ctime(gmtime()));
    int a=5;
    int b=4;
    printf("%d\n",get_power_index(124667459));// 9
    printf("%d\n",my_pow(&a,&b));// 625

    printf("%d\n",is_flower(153));//0
    printf("%d\n",is_flower(875));//1
    printf("%d\n",is_flower(370));//0
    printf("%d\n",is_flower(371));//0
    printf("%d\n",is_flower(407));//0
    printf("%d\n",is_flower(934));//1
    printf("%d\n",is_flower(93084));//0
    printf("%d\n",is_flower(92727));//0
    printf("%d\n",is_flower(54748));//0


    print_time();
    printf("求1- %lf 以内的水仙花数.....\n",N);


    int i;

    clock_t start, finish;
    double duration;

    //printf("%c\n",ctime(gmtime()));

    start = clock(); // 计时开始

    // 穷举算法
    for(i=1;i<N;i++)
    {
        if(is_flower(i) == 0)
        {
            printf("%d位数的水仙花：%d\n",get_power_index(i),i);
        }
    }

    print_time();

    finish = clock(); //计时结束

    //printf("%c\n",ctime(gmtime()));

       /* 测量一个事件持续的时间*/
    duration = (double)(finish - start) / CLOCKS_PER_SEC;
    printf( "Time using is %f seconds\n", duration );
      //system("pause");

    return 0;
}

/*
#include <stdio.h>
#include <time.h>

void main ()
{
time_t rawtime;
struct tm * timeinfo;

time ( &rawtime );
timeinfo = localtime ( &rawtime );
printf ( "\007The current date/time is: %s", asctime (timeinfo) );

exit(0);
}

=================
#include <time.h> -- 必须的时间函数头文件
time_t -- 时间类型（time.h 定义）
struct tm -- 时间结构，time.h 定义如下：
int tm_sec;
int tm_min;
int tm_hour;
int tm_mday;
int tm_mon;
int tm_year;
int tm_wday;
int tm_yday;
int tm_isdst;

time ( &rawtime ); -- 获取时间，以秒计，从1970年1月一日起算，存于rawtime
localtime ( &rawtime ); -- 转为当地时间，tm 时间结构
asctime （）-- 转为标准ASCII时间格式：
星期 月 日 时：分：秒 年
=========================================
你要的格式可这样输出：
printf ( "%4d-%02d-%02d %02d:%02d:%02d\n",1900+timeinfo->tm_year, 1+timeinfo->tm_mon,
timeinfo->tm_mday,timeinfo->tm_hour,timeinfo->tm_min,timeinfo->tm_sec);

就是直接打印tm，tm_year 从1900年计算，所以要加1900，
月tm_mon，从0计算，所以要加1
其它你一目了然啦。
*/



/*
#include<math.h>

函数原型是：
1.double pow(double _X,double _Y);
2.double pow(double _X,int _Y);
3.long double pow(long double _X,long double _Y);
4.long double pow(long double _X,int _Y);
5.float pow(float _X,float _Y);
6.float pow(float _X,int _Y);

原型：extern float pow(float x, float y);

用法：#include

功能：计算x的y次幂。

说明：x应大于零，返回幂指数的结果。

举例：

// pow.c

#include
#include
#include
void main()
{
printf("4^5=%f",pow(4.,5.));
getchar();
}
#include
#include

void main( void )
{
double x = 2.0, y = 3.0, z;

z = pow( x, y );
printf( "%.1f to the power of %.1f is %.1f ", x, y, z );
}


Output

2.0 to the power of 3.0 is 8.0


*/



/*
整型 [signed]int -2147483648~+2147483648
无符号整型unsigned[int] 0~4294967295
短整型 short [int] -32768~32768
无符号短整型unsigned short[int] 0~65535
长整型 Long int -2147483648~+2147483648
无符号长整型unsigned [int] 0~4294967295
字符型[signed] char -128~+127
无符号字符型 unsigned char 0~255
单精度 float 3.4 x 10^（-38）~ 3.4 x 10^（+38）
双精度double 1.7 x 10^（-308）~ 1.7 x 10^（+308）
长双精度 long double 1.7 x 10^（-308）~ 1.7 x 10^（+308）
*/

/*
time() 取得本地时间(日期时间函数)
settimeofday() 设置当前时间戳
mktime() 将时间结构数据转换成经过的秒数
localtime() 获取当地目前时间和日期
gmtime() 获取当前时间和日期
gettimeofday() 获取当前时间
ctime() 将时间日期以字符串格式表示
asctime() 将时间日期以字符串格式表示

printf("%s\n",ctime(gmtime()));





*/



/*
常见水仙花数
水仙花数又称阿姆斯特朗数。
水仙花数

水仙花数
三位的水仙花数共有4个：153，370，371，407；四位的水仙花数共有3个：1634，8208，9474；
五位的水仙花数共有3个：54748，92727，93084；
六位的水仙花数只有1个：548834；
七位的水仙花数共有4个：1741725，4210818，9800817，9926315；
八位的水仙花数共有3个：24678050，24678051，88593477
……
……
使用高精度计算，可以得到超过INT类型上限的水仙花数：
5: 93084
5: 92727
5: 54748
6: 548834
7: 9800817
7: 4210818
7: 1741725
7: 9926315
8: 24678050
8: 24678051
8: 88593477
9: 146511208
9: 912985153
9: 472335975
9: 534494836
10: 4679307774
11: 32164049650
11: 40028394225
11: 42678290603
11: 49388550606
11: 32164049651
11: 94204591914
11: 44708635679
11: 82693916578
14: 28116440335967
16: 4338281769391370
16: 4338281769391371
17: 35875699062250035
17: 21897142587612075
19: 3289582984443187032
19: 4929273885928088826
19: 4498128791164624869
20: 63105425988599693916
21: 449177399146038697307
21: 128468643043731391252
23: 27907865009977052567814
23: 35452590104031691935943
23: 27879694893054074471405
23: 21887696841122916288858
24: 174088005938065293023722
24: 188451485447897896036875
(为环保起见，24位以上的水仙花数略)
理论上，最大的水仙花数不超过34位。
*/