模拟公式

我有这样的想法多线程同时算
这边算a=1*......(5000)
      b=5001*....(10000)
        ......
       ......
    n=95000-100000
   
    x=a*b*....*n
不过多线程在c#学过,现在不会了

#include<stdio.h>
#define N 10000
int main()
{
    long x,y,z,a[N],n,c;
    {
        while(scanf("%ld",&n)>0)
        {
            a[0]=1;
            x=0;
            for(y=1;y<=n;y++)
            {
                c=0;
                for(z=0;z<=x;z++)
                {
                a[z]=a[z]*y+c;
                c=a[z]/N;
                a[z]%=N;
                }
                if(c>0)
                {
                x++;
                a[x]=c;
                }
             }
            // printf("%ld",a[x]);
             for(y=x;y>=0;y--)
                 printf("%.4ld",a[y]);  //补0！！
             printf("\n"); //.4=lg(N);
         }
    return 0;
    }
}
贴个精简代码，非原创

我一直不知道这么大的数是怎样表示出来的，照理说用%d打不出这么多位数字吧，但是我的电脑输入9000运行后得到的数值个数是31682.
还有一个很纳闷的问题是，高精度的表示，用duble表示显示小数点后17位之前表示，17到518位都是0，再往后都是乱码。。。