大阶乘大精确值

#include <stdio.h>
#include <string.h>
const int maxn =3000;
int main()
{
    int f[maxn];
    memset(f, 0, sizeof(f));
    int i,j,n;f[0]=1;
    scanf("%d",&n);
    
    for (i=2; i<=n; i++) {
        int c=0;
        for(j=0;j<=maxn;j++)
        {
           int s=f[j]*i+c;
            f[j]=s%10;
            c=s/10;
        }
    }
    for(j=maxn-1;j>=0;j--) if(f[j]) break;
    for (i=j; i>=0; i--) printf("%d",f[i]);
    printf("\n");
}