C语言求21位水仙花数

/*（编程题）花朵数一个N位的十进制正整数，如果它的每个位上的数字的N次方的和等于这个数本身，则称其为花朵数。例如：当N=3时，153就满足条件，因为1^3+5^3+3^3=153，这样的数字也被称为水仙花数（其中，“^”表示乘方，5^3表示5的3次方，也就是立方）。当N=4时，1634满足条件，因为1^4+6^4+3^4+4^4=1634。当N=5时，92727满足条件。实际上，对N的每个取值，可能有多个数字满足条件。程序的任务是：求N=21时，所有满足条件的花朵数。注意：这个整数有21位，它的各个位数字的21次方之和正好等于这个数本身。如果满足条件的数字不只有一个，请从小到大输出所有符合条件的数字，每个数字占一行。因为这个数字很大，请注意解法时间上的可行性。要求程序在3分钟内运行完毕。*//*本程序比较容易理解。执行时间40秒左右（P6100处理器，2G内存，WIN7系统）。思路：任何一个水仙花数，比如153=1^3+5^3+3^3，里面包括1个1,1个5,1个3。如果是513,315,351这样的数字……，虽然这几个不是水仙花数，但这几个数字中1,3,5出现的次数和水仙花数是一样的……那么513,315,351这些数字，每位上数字的和应该就是水仙花数。如果是从1000……-9999……的枚举法，我想问题肯定大了。而如果采用上面的方法，这样的话，采用对数字出现次数进行枚举会让程序加快很多。*/#includestdio.h#includetime.hintnum_time[10]={0};//num_time[10]用来得到0-9数字出现的个数intsum[10][21]={(0,0)},sumc[10][21]={(0,0)};//sum数组用来计算出0-9的21次方。sumc用来计算数字出现次数*它的21次方intcount=0;//计算得到几个水仙花数了，得到2个以后就结束程序voidinitProgramm()//这个函数用来获得0-9的21次方，存在sum数组里{inti,j,k;for(i=0;i10;i++)sum[i][20]=i;for(k=2;k10;k++){for(i=1;i21;i++){for(j=20;j=0;j--){sum[k][j]*=k;}for(j=20;j0;j--){sum[k][j-1]+=sum[k][j]/10;sum[k][j]%=10;}}}}voidcheck(inti0,inti1,inti2,inti3,inti4,inti5,inti6,inti7,inti8,inti9)//检测数字是不是水仙花数{inti,j;intgetAdd[21]={0};for(j=0;j21;j++)//把数字出现次数*它的21次方{sumc[1][j]=sum[1][j]*i1;sumc[2][j]=sum[2][j]*i2;sumc[3][j]=sum[3][j]*i3;sumc[4][j]=sum[4][j]*i4;sumc[5][j]=sum[5][j]*i5;sumc[6][j]=sum[6][j]*i6;sumc[7][j]=sum[7][j]*i7;sumc[8][j]=sum[8][j]*i8;sumc[9][j]=sum[9][j]*i9;}for(i=0;i10;i++)//进位for(j=20;j0;j--){sumc[i][j-1]+=sumc[i][j]/10;sumc[i][j]%=10;}for(i=0;i10;i++)//得到一个数每位的21次方的和，就是把sumc叠加起来for(j=20;j=0;j--)getAdd[j]+=sumc[i][j];for(i=20;i0;i--)//进位{getAdd[i-1]+=getAdd[i]/10;getAdd[i]%=10;}intj1=0,j2=0,j3=0,j4=0,j5=0,j6=0,j7=0,j8=0,j9=0,j0=0;for(i=20;i=0;i--)//用来判断和里面每个数字出现的次数{switch(getAdd[i]){case0:j0++;break;case1:j1++;break;case2:j2++;break;case3:j3++;break;case4:j4++;break;case5:j5++;break;case6:j6++;break;case7:j7++;break;case8:j8++;break;case9:j9++;break;}}/*如果一个数字，和里0-9出现的次数与这个数字里0-9出现的次数相同，那么和就应该是水仙花数（第一位数字不能为0）*/if((i0==j0)&&(i1==j1)&&(i2==j2)&&(i3==j3)&&(i4==j4)&&(i5==j5)&&(i6==j6)&&(i7==j7)&&(i8==j8)&&(i9==j9)&&(getAdd[0]!=0)){printf(\n);count++;for(i=0;i21;i++)printf(%d,getAdd[i]);printf(\n);}}voidmain(){intt,t1,t2;initProgramm();t1=time(NULL);inti0,i1,i2,i3,i4,i5,i6,i7,i8,i9;for(i9=0;i910;i9++){for(i0=1;i022;i0++){if(count==2)//出现2个水仙花数以后breakbreak;if(i9+i0==21)//几个数字的出现次数和为21以后就break,因为后面的数字出现次数和一定大于21，就超过了21位{check(i0,0,0,0,0,0,0,0,0,i9);break;}for(i2=0;i222;i2++){if(count==2)break;if(i9+i0+i2==21){check(i0,0,i2,0,0,0,0,0,0,i9);break;}for(i3=0;i322;i3++){if(count==2)break;if(i9+i0+i2+i3==21){check(i0,0,i2,i3,0,0,0,0,0,i9);break;}for(i4=0;i422;i4++){if(count==2)break;if(i9+i0+i2+i3+i4==21){check(i0,0,i2,i3,i4,0,0,0,0,i9);break;}for(i5=0;i522;i5++){if(count==2)break;if(i9+i0+i2+i3+i4+i5==21){check(i0,0,i2,i3,i4,i5,0,0,0,i9);break;}for(i6=0;i622;i6++){if(count==2)break;if(i9+i0+i2+i3+i4+i5+i6==21){check(i0,0,i2,i3,i4,i5,i6,0,0,i9);break;}for(i7=0;i722;i7++){if(count==2)break;if(i9+i0+i2+i3+i4+i5+i6+i7==21){check(i0,0,i2,i3,i4,i5,i6,i7,0,i9);break;}for(i8=0;i822;i8++){if(count==2)break;if(i9+i0+i2+i3+i4+i5+i6+i7+i8==21){check(i0,0,i2,i3,i4,i5,i6,i7,i8,i9);break;}for(i1=0;i122;i1++){if(count==2)break;if(i9+i0+i2+i3+i4+i5+i6+i7+i8+i1==21){check(i0,i1,i2,i3,i4,i5,i6,i7,i8,i9);break;}}}}}}}}}}}t2=time(NULL);t=t2-t1;printf(\n%ds\n,t);}