3000!.c

BBS上淘来的一大堆C源代码

                /*这个程序是用字符数组来表示整数,
改变 M 的值可以改变结果位数,为10000时可以,
我算过3000的阶乘,大概9600位,再大了不知会不会内存出错
             用TC2.0运行*/

#include<conio.h>
#include<string.h>
#define M 10000
/********************/
void add(char d1[],char d2[],char d3[])
{
  int x,y,i,j,t;
  char *p;
  for(i=0,t=0;d1[i]!=0&&d2[i]!=0;i++)
  {
    t=d1[i]+d2[i]-2*'0'+t;
    d3[i]=t%10+'0';
    t=t/10;
  }
  p=d1[i]==0?d2:d1;
  for(;p[i]!=0;i++)
  {
    t=p[i]-'0'+t;
    d3[i]=t%10+'0';
    t=t/10;
  }
  if(t==0)d3[i]=0;
  else
  {
     d3[i]=t+'0';
     d3[i+1]=0;
  }
}
/*************/
void input(char d[])
{
  char t[M];
  int i,j;
  gets(t);
  for(i=strlen(t)-1,j=0;i>=0;i--,j++)
  {
     d[j]=t[i];
  }
  d[j]=0;
}
/********************/
void output(char d[])
{
  int i,j;
  if(d[0]==0)putch('0');
  else
  {
    for(i=strlen(d)-1,j=1;i>=0;i--,j++)
    {
       putch(d[i]);
       if(j%1600==0)
       {
      puts("Press any key to continue!\n");
      getch();
       }
    }
  }
}
/******************/
void bitmut(char d1[],int n,char d2[])
{
  int i,t=0;
  if(n>9||n<0){printf("\nerror");exit(0);}
  for(i=0;d1[i]!=0;i++)
  {
     t=(d1[i]-'0')*n+t;
     d2[i]=t%10+'0';
     t=t/10;
  }
  if(t==0)d2[i]=0;
  else
  {
     d2[i]=t+'0';
     d2[i+1]=0;
  }
}
/**********************/
void mut(char d1[],char d2[],char d3[])
{
  int i,j,k,n;
  char t[M],d[M];
  d[0]='0';
  d[1]=0;
  for(i=0;d2[i]!=0;i++)
  {
     n=d2[i]-'0';
     bitmut(d1,n,t);
     for(j=i;j>0;j--)
     {
    for(k=strlen(t)+1;k>0;k--)
           t[k]=t[k-1];
    t[0]='0';
     }

     add(t,d,d3);
     strcpy(d,d3);
  }
}
/**********************/
void discr(char d[])
{
   int i;
   if(d[0]>'0')
     d[0]=d[0]-1;
   else
   {
     if(d[1]!=0)
     {
       discr(d+1);
       d[0]='9';
     }
     else
     goto end;
   }
end:  ;
}
/*****************/
void format(char d[])
{
   int i;
   for(i=strlen(d)-1;d[i]=='0';d[i++]=0);
}
/*********************/
main()
{
  char n[M],r[M],t[M];
  r[0]='1';
  r[1]=0;
  clrscr();
  printf("n=");
  input(n);
  putch('\n');
  output(n);
  puts("!=");
  for(;;)
  {
     strcpy(t,r);
     mut(t,n,r);
     discr(n);
     format(n);
     if(!strcmp(n,"\0"))break;
  }
  output(r);
  getch();
}