阶乘计算器改进版

原帖及讨论：http://bbs.bccn.net/thread-118393-1-1.html

只要你有内存和空间,该程序可无限扩展
默认设置最大计算100000!可精确到每一位
这是以前我做的程序的改进版,速度提升了1/3


/*SK-CHINA 阶乘计算器*/
/*本程序可以准确计算十万以内的阶乘*/
/*结果输出为out.txt*/
#include "stdio.h"
long s[500000]={0L}; /*存储数据*/
long h=499999L; /*数据位置*/
int cheng(long num) {
    long t=499999L;
    while(s[h]==0) h++;
    if(num>=100000L) h--;
    if(num>=10000L) h--;
    if(num>=1000L) h--;
    if(num>=100L) h--;
    if(num>=10L) h--;
    h--;
    while(t!=h) {
      s[t]=s[t]*num;
      t--;
    } t=499999L;
    h=h-5L;
    while(t!=h) {
      s[t-1]+=s[t]/10;
      s[t]=s[t]%10;
      t--;
    }
    h=h+5L;
}
int main(void) {
    int i;
    long x;
    FILE *output;
    s[499999]=1L;
    output=fopen("out.txt","w");
    scanf("%ld",&x); /*输入数据*/
    for(i=1;i<=x;i++) cheng(i);
    for(i=0;i<500000L;i++) {
      if(i%500==0) fprintf(output,"/n");
      fprintf(output,"%ld",s[i]);
    }
    return 0;
}

完整版:
速度较上面的稍慢,但较为友好
/*SK-CHINA 阶乘计算器*/
/*本程序可以准确计算十万以内的阶乘*/
/*结果输出为out.txt*/
#include "stdio.h"
long s[500000]={0L}; /*存储数据*/
long h=499999L; /*数据位置*/
int cheng(long num) {
    long t=499999L;
    while(s[h]==0) h++;
    if(num>=100000L) h--;
    if(num>=10000L) h--;
    if(num>=1000L) h--;
    if(num>=100L) h--;
    if(num>=10L) h--;
    h--;
    while(t!=h) {
      s[t]=s[t]*num;
      t--;
    } t=499999L;
    h=h-5L;
    while(t!=h) {
      s[t-1]+=s[t]/10;
      s[t]=s[t]%10;
      t--;
    }
    h=h+5L;
}
int main(void) {
    int i;
    long x;
    FILE *output;
    s[499999]=1L;
    output=fopen("out.txt","w");
    printf("SK-CHINA 阶乘计算器/n请输入数据:");
    scanf("%ld",&x); /*输入数据*/
    for(i=1;i<=x;i++) { cheng(i); printf("/r当前已计算至%ld!",i); }
    printf("/n计算结束,现在程序正在将结果保存到out.txt/n");
    for(i=0;i<500000L;i++) {
      if(i%500==0) fprintf(output,"/n");
      fprintf(output,"%ld",s[i]);
    }
    printf("写入结束!按任意键退出");
    getch();
    return 0;
}