/*
 *
 *
 *
 */
#include <stdio.h>

void ladd(short *, short *, short *);
void lsub(short *, short *, short *);
void lmul(short *, short, short *);
void printresult(short *);

#define L 100000
#define L2 ((L+3)/4)
int ppt;

int main()
{
  int target;
  printf("enter...>");
  scanf("%d",&target);
  printf("\n");
  static short s[L2];
  short k;
  for (k=0;k<L2;k++)
    s[k]=0;
  s[L2-1]=1;
  for(k=1;k<=target;k++){
    lmul(s,k,s);
    if (k>=target){
    printf("%2d!=",k);
    printresult(s);
    }
  }
}

void lmul(short a[], short b, short c[])
{
  short i; long d, cy=0;
  for (i=L2-1;i>=0;i--){
    d=a[i];
    c[i]=(d*b+cy)%10000;
    cy=(d*b+cy)/10000;
  }
}

void printresult(short c[])
{
  short i;
  int sinri=0,ketasuu=4;
  for (i=0;i<L2;i++){
    if (c[i]!=0)
      sinri=1;
    if (sinri==1){
      printf("%04d",c[i]);
      ketasuu=ketasuu+4;
    }
  }
    printf("\n\n%d",ketasuu);
    printf("\n");
}

void ladd(short a[], short b[], short c[])
{
  short i,cy=0;
  for(i=L2-1;i>=0;i--){
    c[i]=a[i]+b[i]+cy;
    if(c[i]<10000)
      cy=0;
    else {
      c[i]=c[i]-10000;
      cy=1;
    }
  }
}

void lsub(short a[], short b[], short c[])
{
  short i,brrw=0;
  for(i=L2-1;i>=0;i--){
    c[i]=a[i]-b[i]-brrw;
    if(c[i]>=0)
      brrw=1;
    else {
      c[i]=c[i]+10000;
      brrw=1;
    }
  }
}