[公告]挑战C板块的 阶乘算法

[QUOTE]#include <stdio.h>
#include <math.h>
#include <malloc.h>
#include <time.h>
#include <stdlib.h>
#define CARRY 100000
clock_t start,end ;
int GetFactorialMemSize(int n);
long*InitiFactorialMem(int Size);
int Factorial(register long*FactorialMem,int n);
void PrintFactorial(long*FactorialMem,int Start);
int main(void)
{
int N,FactorialStart,FactorialMemSize ;
long*PFactorial ;
char Key='y' ;
while(Key=='y'||Key=='Y')
{
printf("N=?");
while(scanf("%d",&N)!=1)getchar();
FactorialMemSize=GetFactorialMemSize(N)/5+1 ;
PFactorial=InitiFactorialMem(FactorialMemSize);
start=clock();
FactorialStart=Factorial(PFactorial,N);
end=clock();
printf("\n计算%d!花费: %f 秒\n按回车打印....\n\n",N,(float)(end-start)/CLK_TCK);
getchar();
getchar();
start=clock();
PrintFactorial(PFactorial,FactorialStart);
end=clock();
printf("\n\n输出%d!花费: %f 秒\n",N,(float)(end-start)/CLK_TCK);
free(PFactorial);
printf("继续? y 或 Y 退出? 任意键\n");
Key=getchar();
}
return 0 ;
}
/*-----------------------------------------------------------------*/
int GetFactorialMemSize(int n)
{
double sum=1.0 ;
int i ;
for(i=1;i<=n;i++)
sum+=log10((long double)i);
return(int)sum ;
}
/*-------------------------------------------------------------------*/
long*InitiFactorialMem(int Size)
{
long*FactorialMem=(long*)calloc(Size,sizeof(long));
if(!FactorialMem)
{
printf(" Allocation Error !\n");
exit(1);
}
return FactorialMem ;
}
/*------------------------------------------------------------------*/
int Factorial(register long*FactorialMem,int n)
{
int NotZero=0,End=0,i,j ;
long int CARRY_number ;
register long int tem ;
FactorialMem[0]=1 ;
for(i=1;i<=n;i++)
{
CARRY_number=0 ;
for(j=NotZero;j<=End;j++)
{
tem=FactorialMem[j]*i+CARRY_number ;
CARRY_number=tem/CARRY ;
FactorialMem[j]=tem%CARRY ;
}
while(!FactorialMem[NotZero])NotZero++;
if(CARRY_number)FactorialMem[++End]=CARRY_number ;
}
return End ;
}
/*-------------------------------------------------------------------*/
void PrintFactorial(long*FactorialMem,int Start)
{
printf("%ld",FactorialMem[Start--]);
while(Start>=0)printf("%.5ld",FactorialMem[Start--]);
}[/QUOTE]

环境：PII 233 winme
我的


[QUOTE]N=?10000
计算10000!花费: 4.500000 秒
按回车打印....

............（略去）
0000000000000000000000000000000000
0000000000000000000000000000000000
0000000000000000000000000000000000
输出10000!花费: 6.590000 秒
继续? y 或 Y 退出? 任意键[/QUOTE]

你的，计算部分就要

23秒多

纠正一下

哈哈，程序严重有问题。
kai 的10000！的前几位280202.........
knocker 的10000！的前几位4393181........
网上公布的测试过的 10000！
2846259680
9170545189 0641321211 9868890148 0514017027 9923079417
9994274411 3400037644 4377299078 6757784775 8158840621
.............
这程序你们都有底吗，到底谁对了！

没弄明白！

进一步的研究发现，knocker的算法有抄袭乌鸦丘比特的嫌疑，只是对他的算法的无关紧要部分进行了美化
加入了华丽的内存分配，但是主体部分是显然一样的！

不要斗，to kai，有没有看过《C++性能优化***》（忘了书名），反正是蓝色封面很薄的一本书，里面说到，C++比起C语言有一定的起伏，就是，有些命令比C快，有些比C慢，而且慢N倍，其中的绞绞者就是类，用类的话，就算全部inline也不及C快，所以没得比。

这种变量单一的算法题体现不了C++的优点，要用我以前问过那种求路径最长的旅游那种N个变量N个函数N种的关系的，用C就会做到头晕。

哈哈哈哈，学过计算机组成原理后就知道，如果想比快，没有任何东西比精练的汇编来得凶悍了~~

To: kai 版主的二次方

假如你记性不好，我可以帮你回忆一下你设立这个贴子的目的：你是为了打击C版的学弟学妹们，你认为他们写的程序计算50次1000!竟然要7秒多，真是烂得可以。所以，你还提出一个口号：挑战速度的极限。
现在，你用什么通用性来搪塞我，只是为不给我3000分找理由！而且我们也说定比的是10000!的速度，不是10001!也不是9999!所以，你的理由不成立。
不给我3000分，你就对不行C版的学弟学妹！你是我的偶像，不给我3000分，你在我心中
的光辉形像就大打折扣了。^_^

To: 神vLinux飘飘 前C高手版版主

你说的是否正确，我愚味不知道。你就用汇编写，我用C，我们俩个再比。可是，你会汇编么？哈哈，哈哈，要么你用JAVA，JAVA你会的。哈哈，哈哈,B4你这个PT。

To: 激情依旧 现任版主

C++方便是因为别人把你认为烦的事都处理好了。如果你用VB，BCB写程序还会更方便。

草，你去死吧。。。我在看dbx协议没空理你。

kai居然屈服，钱是事小，给老k的烂舌说服了事大。

不愤中......