三千大数阶乘

1. #include   <stdio.h>  
2.  
3. #include   <conio.h>  
4.  
5. #define   STRMAX   2000  
6.  
7. #define   P   100000  
8.  
9. #define   N   3000  
10.   
11. int   main()  
12.   
13. {  
14.   
15. long   t[STRMAX]={0},a,b,c,g=0,x,e,p=P;  
16.   
17. t[STRMAX-1]=1;  
18.   
19. for   (e=1;e <=N;e++)  
20.   
21. for(x=STRMAX-1;x>=0;c=t[x],c*=e,a=c/p,b=c%p,t[x]=b+g,g=a,x--);  
22.   
23. while(t[++x]==0);  
24.   
25. printf( "%ld ",t[x]);  
26.   
27. while(++x <STRMAX)  
28.   
29. printf( "%.5ld ",t[x]);  
30.   
31. getch();return   0;  
32.   
33. }