大数阶乘

#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
unsigned short  s[20000];//最多可以容纳5000阶乘；且必须为int型，当为short型时，数组元素大小溢出。因为c为int型，当不段的
int main()
{
    int m;
    while(~scanf("%d",&m)) 
    {
        memset(s,0,sizeof(s));
        s[0]=1;
        int c,l=0;
        for(int i=1;i<=m;i++)//m的阶乘
        {
            c=0;
            for(int j=0;j<=20000;j++)//乘以i之后结果
            {
                s[j]=s[j]*i+c;
                c=s[j]/10;
                s[j]%=10;
            }
        }
        //printf("%d",s[l]);
        int i;
        for(i=20000;i>=0;i--)
            if(s[i])
            break;
        for(;i>=0;i--)
            printf("%d",s[i]);
        printf("\n");
    }
    return 0;
}