#include<stdio.h>
#include<time.h>
#include<string.h>
//=======begin=======
#define MD 30
static int a[10][MD][MD];//0-9的10-29次方的值表
static x[MD],y[MD],count;
char res[10][MD];
int listBD(int a[])//显示大数,a[MD-1]记录位数
{
int i;
for(i=a[MD-1]-1;i>=0;i--) printf("%d",a[i]);printf("\n");
return 0;
}
int mulBD(int a[],int x)
{
int i;
for(i=0;i<a[MD-1];i++) a[i]=a[i]*x;
for(i=0;i<a[MD-1];i++)
if(a[i]>=10)
{
a[i+1]=a[i+1]+a[i]/10;
a[i]=a[i]%10;
}
if(a[a[MD-1]]>0) a[MD-1]++;
return 0;
}
int atob(int a[],int b[])
{
int i;
for(i=0;i<MD;i++) b[i]=a[i];
return 0;
}
int addb(int a[],int b[])
{
int i,m;
if(a[MD-1]>b[MD-1]) m=a[MD-1];
else m=b[MD-1];
a[MD-1]=m;
for(i=0;i<m;i++) a[i]=a[i]+b[i];
for(i=0;i<m;i++)
if(a[i]>=10)
{
a[i+1]=a[i+1]+a[i]/10;
a[i]=a[i]%10;
}
if(a[m]>0) a[MD-1]++;
return 0;
}
int init(int n)
{//0-9的n次方的值表
int i,j;
for(i=0;i<10;i++) {a[i][1][0]=i;a[i][1][MD-1]=1;}
for(i=0;i<10;i++)
for(j=2;j<=n;j++)
{
atob(a[i][j-1],a[i][j]);
mulBD(a[i][j],i);
}
return 0;
}
int list(int N)
{//枚举出0-9组成N位串的所有情况,要求,后面数字不小于前者
int a[MD]={0};
int i,j,p,n=0;
a[0]=0;
j=0;
while(j<N)
{//
if(j==N-1)//最后一位
{
//for(i=0;i<N;i++)printf("%d",a[i]);printf("\n");
//判断是否为自幂数;
judge(a,N);//if(p) res[count++]=p;
if(a[j]<9) {a[j]++;continue;}
else
while(j>0)
{
j--;
if(a[j]<9){a[j]++;break;}
}
}
if(j<N) {j++;a[j]=a[j-1];}
if(a[0]==9)break;
}
return 0;
}
int sort(int s[],int N)
{
int i,j,x;
for (i=0;i<N-1;i++)
for(j=0;j<N-1-i;j++)
{
if(s[j]>s[j+1])
{
x=s[j+1];s[j+1]=s[j];s[j]=x;
}
}
return 0;
}
int judge(int b[],int N)
{
int ss[MD],x,j=0;
int i,sum[MD];
for(i=0;i<MD;i++) sum[i]=0;
for(i=0;i<N;i++) addb(sum,a[b[i]][N]);
if(sum[MD-1]!=N)return 0;
atob(sum,ss);
sort(sum,N);
x=1;//假定是,比较两个串是否一致
for(i=0;i<N;i++) if(b[i]!=sum[i]) x=0;
if(x)
{for(i=0;i<N;i++)
res[count][i]=ss[N-i-1]+'0';
res[count][N]='\0';
count++;
} //listBD(ss);
return 0;
}
int main(void)
{
int n=10,i,j,start;
scanf("%d",&n);
if(n<10||n>28)
{
printf("error!");return 0;
}
start=clock();
init(n);
//for(i=2;i<10;i++)listBD(a[i][n]);
list(n);
for (i=0;i<count-1;i++)
for(j=0;j<count-1-i;j++)
{
if(strcmp(res[j],res[j+1])>0)
{
strcpy(res[count],res[j+1]);strcpy(res[j+1],res[j]);strcpy(res[j],res[count]);
}
}
for(i=0;i<count;i++)puts(res[i]);
printf("用时:%d ms\n",clock()-start);
//========end========
return 0;
}
编写一个能计算10到21位自幂数的程序。
- 超大整数的存储与运算。
- 如何减少搜索空间。
- 一个数的 21 次幂极为庞大,即使 unsigned long 也显得无能为力;可以借助数组存储超大整数,用逐位运算实现,加、乘等运算。
- 即使是 C 语言这样运行效率极高的语言,穷举判断也很有可能超出 3 分钟的时限。
- 为了进一步减少程序运行时间,我们可以预存 0 到 9 的 21 次幂的结果,之后采用查表的方法求各位数的 N 次幂之和。
- 一种减少搜索空间的办法是,用0-9的N的重复组合,并进行排序,即后面的数字不小于前面的数字,如000,001,002,。。。,009,011,
-
运行结果示例:国家级一流本科课程《C语言程序设计》实践教学正式上线头歌教学平台!欢迎大家指导!https://forge.educoder.net/paths/uiqwncjx