1 #include<stdio.h>
2 #include<time.h>
3 #include<string.h>
4 //#include<conio.h>
5 //#include<math.h>
6 //有的注释 是用来测试数据
7 #define N 21
8 int a[10][N+1];
9 void fun()
10 {
11 int i,temp,j,k,t;
12 memset(a,0,sizeof(a));
13 for(i=1;i<10;i++)
14 {
15 a[i][0]=1;
16 for(j=0;j<N;j++)
17 {
18 temp=0;
19 for(k=0;k<N+1;k++)
20 {
21 t=a[i][k]*i+temp;
22 temp=t/10;
23 a[i][k]=t%10;
24 }
25 }
26 }
27 for(i=0;i<10;i++)
28 {
29 for(j=N-1;j>=0;j--)
30 if(a[i][j]!=0)break;
31
32 a[i][N]=j+1; //记录从第几个开始之后全部为0
33 }
34 //可以在此处输出结果以供检查错误
35 for(i=0;i<10;i++)
36 {
37 printf("i=%d\n",i);
38 for(j=0;j<N+1;j++)
39 printf("%3d",a[i][j]);
40 printf("\n");
41 }
42 }
43 void fun1()
44 {
45 int b[N];
46 int c[N]={0},d[N]={0},e[10]={0};//c[N]存b[N]中每个数字的21次方之和
47 int x,i,j,k,temp,flag; //d[N]存c[N]从大到小排好序的值
48 int f[100][21]; //用来存放满足条件的数 然后从小到大输出
49 for(i=0;i<9;i++) //b[N]存从99……99到100……00之间的数 初始化为9个9后面全是8,因为10个9的21次方之和将超过21位数
50 b[i]=9;
51 for(i=9;i<N;i++)
52 b[i]=8;
53 k=0;//先将满足条件的数组个数清为0
54 while(b[0]!=0) //e[N]存中间排序的中间量 选择的排序方法比较特别（已知最大值的排序）
55 {
56 flag=1;
57
58 //求21个数的21次方之和
59 for(i=0;i<N;i++)
60 {
61 temp=0;
62 for(j=0;j<N;j++)
63 {
64 temp=temp/10+a[b[i]][j]+c[j];
65 c[j]=temp%10;
66 if(j==N-1&&temp>9) //
67 {
68 flag=0;
69 //printf("超过21位数\n");
70 break;
71 }
72 if(a[b[i]][j]==0&&j>=a[b[i]][N])break;//
73 }
74 if(flag==0)break;// 这些都是为了节省时间 进行程序优化
75
76 }
77 //getch();
78 if(c[N-1]==0)flag=0;
79 if(flag)
80 {
81 //将c[N]排序 注意排序方法 复杂度较低
82 x=0;
83 for(i=0;i<N;i++)
84 e[c[i]]++; //分别有几个0到9 存到e[N]中 曾将c[i]错写成b[i]~~
85 for(i=10-1;i>=0;i--)
86 {
87 for(j=0;j<e[i];j++)
88 d[x++]=i;
89 }
90
91 //比较二者是否相等
92 for(i=0;i<N;i++)
93 {
94 if(b[i]!=d[i])
95 {
96 flag=0;
97 break;
98 }
99 }
100 }//if
101 //如果标记falg=1则输出结果
102 if(flag)
103 {
104 /*printf("输出结果为:");
105 for(i=N-1;i>=0;i--)
106 printf("%d",c[i]); 将c[i]曾错写成b[i]
107 printf("\n");若这样输出则为从大到小输出*/
108 j=0;
109 for(i=N-1;i>=0;i--)
110 f[k][j++]=c[i];
111 k++;
112 }
113 //找下一个a[N]
114 for(i=N-1;i>=0;i--)
115 if(b[i]!=0) //此处曾经出错 将b[i]写成a[i]
116 break;
117 temp=b[i];
118 // printf("temp=%d ",temp);
119 while(i<N)
120 {
121 b[i++]=temp-1; //保证去掉重复的数 并按从大到小顺序 进行处理
122 }
123 //printf("b[2]=%d b[3]=%d\n",b[2],b[3]);
124
125 memset(c,0,sizeof(c));//全部置0 因为后面还要用
126 memset(d,0,sizeof(d));
127 memset(e,0,sizeof(e));
128 }
129 //从小到大输出
130 for(i=k-1;i>=0;i--)
131 {
132 for(j=0;j<N;j++)
133 printf("%d",f[i][j]);
134 printf("\n");
135 }
136 }
137 int main()
138 {
139 fun();
140 // fun1();
141
142 printf("\n程序运行了%.2lf秒\n",(double)clock()/CLOCKS_PER_SEC);//只是为了测试一下程序运行的时间
143 return 0;
144 }
