acm-大数阶乘

#include <iostream>  
#define N 18000        
#include<cstring>  
using namespace std;  
int a[N];     
int main( )  
{  
 int n,i,j,c,s;  
 cin>>n;  
 memset(a,0,sizeof(a));      
 a[0]=1;                           
 for(j=2;j<=n;j++)   
 {  
  c=0;                  
  for(i=0;i<N;i++)  
  {  
   s=a[i]*j+c;  
   a[i]=s%10;  
  c=s/10;  
  }  
 }  
 for(i=N-1;i>=0;i--)  
 if(a[i])break;           
 for(j=i;j>=0;j--)    
 cout<<a[j];  
   
 cout<<endl;  
 return 0;  
}