factorial  계승n! 의 값을 다자리수로 산출한다

int factorial(int *data, int n);

data  (출력) n! 의 값이 10진 4자리수씩 들어간 정수의 배열
n     (입력) 계승n의 값

배열 data에 계산값이 들어온 요소의 수(배열의 선두로부터)

int factorial(int *data, int n)
{
    int  i, j;
    int  last;
    int  *dap, *dap2;
    int  t;
    long x;
    
    *data = 1;
    last = 0;
    for (i = 1; i <= n; i++) {
        t = 0;
        for (j = 0, dap = data; j <= last; j++) {
            x = *dap;
            x = x * i + t;
            *dap++ = x % 10000;
            t = x / 10000;
        }
        if (t ! = 0) {
            *dap = t;
            last++;
        }
    }
    
    for (dap = data, dap2 = data + last; dap < dap2; dap++, dap2--) {
        t     = *dap;
        *dap  = *dap2;
        *dap2 = t;
    }
    return last + 1;
}

n! 의 자리수는 스털링의 근사식 $n!\; =\; nne-n(2n\pi )1/2$보다 ,
　　　$n\; log$₁₀(n/e) + (log₁₀(2nπ)/2+1
그리고 요구된다.