如何求一个大数的N！

#include<iostream>
using namespace std;
#define I 10000000000 
__int64 a[35670]; 
void fac(__int64 n,__int64 *a)
{ 
    __int64 i,j, k,t;
    k=a[0]; 
    for(i=2;i<=n;++i)
        a[i]=0; 
    for(i=2;i<=n;i++)
    { 
        for(j=1;j<=k;++j)
        {
            a[j]*=i;
        }
        a[0]=0;
        for(t=k;t>=1;--t)
        {
            if(a[t]>=I)
            {
                a[t-1]+=a[t]/I;
                a[t]%=I;
            }
        }
        if(a[0])
        {
            k++;
            for(t=k;t>=1;--t)
                a[t]=a[t-1];