大整数的阶层

#include"stdio.h"
#include"string.h"
#include<stdlib.h>
int main()
{
char n[10];
int a[20000]={0},b[20000]={0},c[20000]={0},temp;
int i,j,k,len;
         printf("请输入阶乘数:");
scanf("%s",&n);
/*下面以b为乘数,a为被乘数(依次从1到n递增)*/
temp=atoi(n);/*化整*/
b[0]=1;/*初始化b,求阶乘从1开始*/
len=1;/*初始化b的长度*/
for(i=1;i<=temp;i++)/*阶乘从低往高乘*/
{
sprintf(n,"%d",i);/*将阶乘项转换为字符串*/
for(j=strlen(n)-1;j>=0;j--)/*逆向存储阶乘项(串)*/
a[strlen(n)-1-j]=n[j]-'0';/*注意要变成整数才能进行相乘*/
for(j=0;j<len;j++)/*取b的所有位*/
for(k=0;k<strlen(n);k++)/*取阶乘项的每一位*/
  {
c[j+k]+=b[j]*a[k];/*用c来存储每一位相乘的结果*/
c[j+1+k]+=c[j+k]/10;/*若c[j+k]不小于10,则将大于10的部分送至高位*/
c[j+k]%=10;/*送至高位后,c[j+k]应作相应的求余运算。只保留剩余部分*/  }
if(c[j+k-1]!=0)/*判断新的b数的位数长度*/
len=j+k;
  else
  len=j+k-1;
  for(j=0;j<len;j++)/*将b更新*/
b[j]=c[j];
  memset(c,0,sizeof(c));/*将c归零*/
}
for(i=len-1;i>=0;i--)/*逆向输出结果*/
printf("%d",b[i]);
printf("\n");
return 0;
}