ZOJ 1133 Smith Numbers && POJ 1142 Smith Numbers

史密斯数

   美国有一位数字家名叫阿尔伯特·威兰斯基，他姐夫史密斯非常喜欢研究数学，所以两人经常在一起研讨各种数学问题。有时，两人碰不到一起，就习惯性地用电话交流。

　　一天，两人刚结束电话交谈，史密斯突然灵感来临，对威兰斯基的电话号码“4937775”产生 了兴趣，总觉得这是个特别的数。可它的特殊之处究竟在哪儿呢？史密斯开始思索考证起来，他先把4937775分解质因 数：4937775=3×5×5×65837，然后再把4937775所有质因数各位上的数字相加得：3＋5＋5＋6＋5＋8＋3＋7=42，接着他又把 4937775各位上的数字相加得：4＋9＋3＋7＋7＋7＋5=42 　　

   秘密终于找到了，原来这两个和相等。这真有意思，难道是巧合么？有没有其他的数也有此特点呢？史密斯不能解答这个问题。不过，他的这一发现引起了许多数学家的浓厚兴趣。数学家们带着这个疑问，纷纷对这一现象进行了研究。

　　结果发现，有许多数具有这样独特的性质，其中最小的数是4。大家不妨检查一 下，4=2×2，2＋2=4。类似有，22=2×11,2＋2=2＋1＋1；27=3×3×3，2＋7=3＋3＋3。随着研究的不断深入，他们发现，在 0至1000之间，共有376个这样的数；并且估计在0至10000之间有3300个这样的数。因为这些有趣数的发现，追本溯源是史密斯的功劳，所以数学 家们把这样的数叫做“史密斯数”。

来自百度百科 http://baike.baidu.com/view/1596160.htm?fr=ala0_1

----------------背景介绍完毕------------------------------

 

这题踩着悬崖边过的,悲剧阿.编译器选择c++ 超时,选择gcc 900+ms,看来gcc产生的代码效率比较高阿.
直接上代码,这题也让我郁闷了下,总是float point error  最后初始化了一下,解决了.

#include <stdio.h> #include <string.h> #include <math.h> #define N 1111111 int a[N+2]; /*--------------计算一个数各位相加的和-----------*/ int sum_num1(int n) { int sum = 0; while (n) { sum += n%10; n /= 10; } return sum; } /*------OVER----------*/ /*---------------计算一个数分解质因数之后各位的和-----------------*/ int sum_num2(int p[],int n) { int r,sum=0,i,nn=n; for (i=1;p[i]<=n/2;i++) { r = n % p[i]; if (!r) { n /= p[i]; sum += sum_num1(p[i]); i--; } if (n<=1) break; } if (nn!=n) sum += sum_num1(n); return sum; } /*--------------OVER----------------*/ /*-----------计算史密斯数----------*/ int smith(int *p,int n) { while ( sum_num1(n) != sum_num2(p,n) ) { n++; } return n; } /*---------------计算史密斯数结束-------------*/ int main(void) { /*----------------筛素数----------------------*/ int p[N+2]; memset(p,2,sizeof(p)); int i,j;int n; memset(a,1,sizeof(a)); a[0] = 0;a[1]=0; for (i = 2 ; i < sqrt(N); i++) if (a[i]) for ( j = 2 ; i*j < N ; j++) { a[i*j] = 0; } j=1; for (i=0;i<N;i++) if (a[i]) { p[j] = i; j++; } /*---------------筛素数结束----------------*/ while (scanf("%d",&n) != EOF && n) { n++; printf("%d/n",smith(p,n)); } return 0; }