51NOD N!

#include<iostream>
#include<string.h>
#define MAX 40000
using namespace std;


int main()
{
    int n;
int i,j;
    int a[MAX];      //存数运算结果
int p,h,temp;          //p存储当前结果的位数,h为进位
memset(a,0,sizeof(a));
a[0]=1;
p=1;
cin>>n;
for(i=2;i<=n;i++)   //循环与2,3,4.....n相乘
{
 
for(j=1,h=0;j<=p;++j)    //让a[]的每位与i相乘
{
temp=a[j-1]*i+h;
h=temp/10;//进位
a[j-1]=temp%10;//余数存入a[j]
}
while(h>0)         //如果h不为0,需要进一位(h不可能大于10)
{
a[++p-1]=h%10;
h=h/10;
}
}
for(i=p-1;i>=0;--i)//反输出
{
cout<<a[i];
}
cout<<"\n";
 

}

#include <cstdio>
#include <algorithm>
#include <string>
#include <cstring>
#include <vector>
#include <map>
#include<stdio.h>
#include<math.h>
using namespace std;
#define ll long long
#define N 1000005
#define mod 100000000000000
ll ans[10000000];
int main()
{
    ll n;
    scanf("%I64d",&n);
        ll i,j;
        ll l=0;
        ans[0]=1;
        for(i=1; i<=n; i+=1)
        {
            ll c=0;
            for(j=0; j<=l; j+=1)
            {
                ll t=ans[j]*i+c;
                ans[j]=t%mod;
                c=t/mod;
            }
            if(c!=0)
            {
                ans[++l]=c;
            }
        }
        printf("%I64d",ans[l]);
        for(i=l-1; i>=0; i-=1)
        {
            printf("%0.14I64d",ans[i]);
        }
    return 0;
}