4行代码计算9999的阶乘

本文来源于csdn2001年的一个帖子，题目是有谁可以用四行代码求出1000000的阶乘！！！，见
http://topic.csdn.net/t/20001226/11/52263.html，我在该帖中给出了代码。但严格的将，这个代码的行数并未控制在4行以内。这次给出一个只用4行代码计算10000以内的阶乘的完整程序，并给出代码优化（压缩)过程）

   下面看看这个程序的最初形式

 

1. //第一版，正常的格式 
2. #define N 10    //计算10的阶乘，最多只能算到9999 
3. long a[8916]={1,0},n,i,c,len;
4. int main() 
5. {
6.     n=N;
7.     for ( len=1;n>1; n--)
8.     {
9.         for (c=0,i=0; i<len;i++)
10.         {
11.             long p= a[i]*n+c;
12.             a[i]= p % 10000;
13.             c=p / 10000;
14.         }
15.         a[i]= c;
16.         if (c>0)
17.             len++;
18.     }
19. 
    
20.     printf("%d",a[len-1]);
21. 
    
22.     for( len--;len>=0;len--)
23.         printf("%04d",a[len]);
24.     return 0;
25. }

代码的首次改进，压缩代码行数

1. // 第二版： 将打印语句缩写到1行,同时化简计算部分的内循环
2. #define N 100   //计算9999的阶乘，最多只能算到9999
3. long a[8916]={1,0},n,i,c,len;
4. int main() 
5. {
6.     n=N;
7.     for ( len=1;n>1; n--)
8.     {
9.         for (c=0,i=0; i<len;i++)
10.         {
11.             a[i]= ( c+= a[i]*n ) % 10000; c/=10000;
12.         }
13.         a[i]= c;
14.         if (c>0)
15.         len++;
16.     }
17. 
    
18.     for( len--,printf("%d",a[len--]);len>=0; len--) printf("%04d",a[len]);
19.         
20.     return 0;
21. }

进一步改进

1. // 版本3：将外循环的最后2条语句改为一个 "?" 语句 将
2. #define N 10    //计算9999的阶乘，最多只能算到9999
3. long a[8916]={1,0},n,i,c,len;
4. int main() 
5. {
6.     n=N;
7.     for ( len=1;n>1; n--)
8.     {
9.         for (c=0,i=0; i<len;i++ )
10.         {
11.             a[i]= ( c+= a[i]*n ) % 10000; c/=10000;
12.         }
13.         
14.         ((a[i]=c)>0)?len++:0;
15. 
    
16.     }   
17.     for( len--,printf("%d",a[len--]);len>=0; len--) printf("%04d",a[len]);
18.         
19.     return 0;
20. }

再次改进

1. // 版本4
2. // 将内循环中的两条语句变成一条语句，用逗号分隔,去掉{}
3. // 将外循环最后一行 和 外循环 for语句第三部分 "n--"合并,
4. // 进行格式重排，将main函数和变量定义语句放入一行
5. // 将将 len从1开始改为从0开始
6. 
   
7. #include <stdio.h>
8. #define N 10 //计算N的阶乘，修改N的定义可计算10000以内任意数的阶乘    
9. long a[N]={1,0},n,i,c,len; void main() {
10. for ( n=N,len=0;n>1;((a[i]=c)>0) ? (len++,n--) : n--)
11.     for (c=0,i=0; i<=len;i++ )  a[i]= ( c+= a[i]*n ) % 10000, c/=10000;
12. for( printf("%d",a[len--]);len>=0; len--) printf("%04d",a[len]); }

版本5，进一步优化，不算include, #define 及回车换行符，只有170个byte. 主要改进有：
1. 删除多余的空格
2. "a[8916]={1,0}" 改为 "a[N]={1}",节省5个字符
3. 变量赋值 n=N 放到变量声明语句，减少1个字符
4. len 改为单变量m, 减少了10个以上的字符
5. 全局变量自动初始化为0，len=0,这句去掉
6. "((a[i]=c)>0)?(m++,n--):n--" 改为 "(a[i]=c)>0?m++:0,n--"，节省了6个字符
7. 最后一个for 语句中，由"printf("%d",a[m--]);m>=0;m--" 调整为 "printf("%d",a[m]);--m>=0;", 节省了4个字符
8. 虽然说long 改为int 可节约一个字符，但这样的话，在16为编译器的运行结果将出错，失去了通用性

1. //版本5
2. #include <stdio.h> 
3. #define N 10 //计算N的阶乘，修改N的定义可计算10000以内任意数的阶乘     
4. long a[N]={1},n=N,i,c,m;void main(){
5. for(;n>1;(a[i]=c)>0?m++:0,n--)
6. for(c=0,i=0;i<=m;i++)a[i]=(c+=a[i]*n)%10000,c/=10000;
7. for(printf("%d",a[m]);--m>=0;)printf("%04d",a[m]);}