求n!的程序(n=1&&n=1000)

核心提示：这个程序写得比较乱，没好好整理.因为数字大于700的话就要超出屏幕了，求n!的程序(n=1&&n=1000)，所以就用科学表示方法了,其实程序有点浪费空间的,有些地方还可以改进的./*跳蚤侦察2003.9.26*/#define N 3000void PR1(int *num)/*科学方式输出阶层*/{ int i,j

　　这个程序写得比较乱，没好好整理.因为数字大于700的话就要超出屏幕了，所以就用科学表示方法了,其实程序有点浪费空间的,有些地方还可以改进的.
/*跳蚤侦察2003.9.26*/
#define N 3000
void PR1(int *num)/*科学方式输出阶层*/
{
 int i,j;
 for(i=0;;i++)
  if(num[i]!=0)
　 {
　  printf("%d",num[i]);
　  printf(".");
　  for(j=i+1;j<=i+11;j++)
　  printf("%d",num[j]);
　  printf("e+");
　  printf("%d\n",N-i-1);
　  break;
　 }
}
void pr(int *num)/*输入结果*/
{int i,j;
 for(i=0;;i++)
  if(num[i]!=0)
　 {
　 for(j=i;j<N;j++)
　  printf("%d",num[j]);
　 break;
　 }
 printf("\n");
}
fun(int num[],int numx[],int x,int y)/*相乘结果函数*/
{int t[N],i,j,k,tt[N]={0};/*tt数组用来存放竖式的每一行*/
 for(i=0;i<N;i++)/*用来反复乘第一行*/
  {t[i]=num[i];num[i]=0;}
 for(i=3;i>=3-y;i--)/*1位到3位的乘法*/
  {
　 for(j=N-1;j>=N-x;j--)
　  tt[j+i-3]=numx[i]*t[j];/*记录数式*/
　 for(j=N-1;j>=1;j--)/*这里其实可以用个变量来算每一行到第几位,这样可以节约时间*/
  if(tt[j]>=10)
　 {
　  tt[j-1]+=tt[j]/10;
　  tt[j]=tt[j]%10;
　 }
　 for(j=0;j<N;j++)/*竖式相加得出最后乘出来的结果*/ 
　 {num[j]+=tt[j];tt[j]=0;}
  }
 for(i=N-1;i>=N-x;i--)/*结果进位*/
  if(num[i]>=10)
　 {
　  num[i-1]+=num[i]/10;
　  num[i]=num[i]%10;
　 }
}
 main()
{int num[N]={0},numx[4]={0},i,x,y;
 int jx,jy=0;
 clrscr();
 printf("please input n: ");
 scanf("%d",&jx);
 num[N-1]=1;/*一直用来乘的变量初始值*/
 while(1)
 {
 jy++;
 numx[3]++;
 if(numx[3]==10)
 {numx[3]=0;numx[2]++;}
 if(numx[2]==10)/*个位数进位*/
 {numx[2]=0;numx[1]++;}
 if(numx[1]==10)/*十位数进位*/
 {numx[1]=0;numx[0]++;}
 for(i=0;i<N;i++)
  if(num[i]!=0)
　 break;
 x=N-i;/*求结果的位数*/
 for(i=0;i<4;i++)
  if(numx[i]!=0)
　 break;
 y=4-i;/*求阶层的位数*/
 fun(num,numx,x,y);/*传递过去的两个数,x是结果的位数,y是阶层的位数*/
 if(jy==jx)/*阶层大小的判定*/
  break;
 }
 printf("%d!=",jx);/*输出结果*/
 if(jx<=700)
 pr(num);
 else
 pr1(num);
 getch();
}