Chinaunix首页 | 论坛 | 博客
  • 博客访问: 103543364
  • 博文数量: 19283
  • 博客积分: 9968
  • 博客等级: 上将
  • 技术积分: 196062
  • 用 户 组: 普通用户
  • 注册时间: 2007-02-07 14:28
文章分类

全部博文(19283)

文章存档

2011年(1)

2009年(125)

2008年(19094)

2007年(63)

分类: C/C++

2008-05-25 00:07:29

帖及讨论:

只要你有内存和空间,该程序可无限扩展
默认设置最大计算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;

阅读(1359) | 评论(0) | 转发(0) |
给主人留下些什么吧!~~