１０－2８位水仙花－－国家级一流本科课程《C语言程序设计》实践教学正式上线头歌教学平台！欢迎大家指导！https://forge.educoder.net/paths/uiqwncjx

#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组成Ｎ位串的所有情况，要求，后面数字不小于前者
	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位自幂数的程序。

1. 超大整数的存储与运算。
2. 如何减少搜索空间。

• 一个数的 21 次幂极为庞大，即使 unsigned long 也显得无能为力；可以借助数组存储超大整数，用逐位运算实现，加、乘等运算。
• 即使是 C 语言这样运行效率极高的语言，穷举判断也很有可能超出 3 分钟的时限。
• 为了进一步减少程序运行时间，我们可以预存 0 到 9 的 21 次幂的结果，之后采用查表的方法求各位数的 N 次幂之和。
• 一种减少搜索空间的办法是，用0-9的Ｎ的重复组合，并进行排序，即后面的数字不小于前面的数字，如000，001，002，。。。，009，011，

• 运行结果示例：国家级一流本科课程《C语言程序设计》实践教学正式上线头歌教学平台！欢迎大家指导！https://forge.educoder.net/paths/uiqwncjx