从题目1000的阶方想到... ...

x2

http://www.physicsforums.com/archive/index.php/t-126744.html
http://en.wikipedia.org/wiki/Factorial

[ 本帖最后由 x2 于 2009-4-25 21:33 编辑 ]

reiase

原帖由 x2 于 2009-4-25 20:43 发表
有没有人能看懂：n! = \Gamma(n+1) = \int_0^\infty t^n e^{-t} dt

latex代码

win_hate

原帖由 x2 于 2009-4-25 20:43 发表
有没有人能看懂：n! = \Gamma(n+1) = \int_0^\infty t^n e^{-t} dt

在 dt 之前插入一个 \, 效果会好一些。见 ``The TeXbook'' 第 18 章 ``Fine Points of Mathematics Typing''.

皇家救星

记得前几天csdn上面也有人在讨论这个

当时有人用sql编程写了一个　耗时不到1s

retuor

1000! 用什么语言算都用不了 1s 吧？

1. 
   
2. print reduce (lambda x, y:x*y, range(1,1001))
   

复制代码

real        0m0.082s
user        0m0.048s
sys        0m0.004s

reiase

原帖由 retuor 于 2009-4-26 11:39 发表
1000! 用什么语言算都用不了 1s 吧？


print reduce (lambda x, y*y, range(1,1001))


real        0m0.082s
user        0m0.048s
sys        0m0.004s

这代码交上去，人家不认

x2

原帖由 retuor 于 2009-4-26 11:39 发表
1000! 用什么语言算都用不了 1s 吧？


print reduce (lambda x, y*y, range(1,1001))


real        0m0.082s
user        0m0.048s
sys        0m0.004s

就是1000!的定义

1. 
   
2. (require-extension (srfi 1))
   
3. (apply * (iota 1000 1))
   

复制代码

x2

原帖由 reiase 于 2009-4-26 19:58 发表


这代码交上去，人家不认

可能题目要求自己实现大数的乘法运算。

action08

就是“自己实现大数的”意思，大家可以参考YZFY的
http://blog.pfan.cn/yzfy/25193.html


我还要继续研究，

比如，1x2x3....x1000会是个多少位的数，如果10000呢？
首先得考虑要选多少B（字节）来来存储它，


同时谢谢8楼的大虾给的参考资料

[ 本帖最后由 action08 于 2009-4-27 09:52 编辑 ]

retuor

python 不认就用 c 写一个好了

1. 
   
2. #include <stdio.h>
   
3. 
   
4. unsigned long b=1000000000UL;
   
5. unsigned long r[286];
   
6. 
   
7. void
   
8. factorial (long n, unsigned long *r)
   
9. {
   
10.     long i, j, d;
    
11.     unsigned long c;
    
12.     unsigned long long x;
    
13.    
    
14. 
    
15.     r[0]=1, d=0;
    
16.     for (i=2; i<=n; i++)
    
17.     {
    
18.         c=0;
    
19.         for (j=0; j<=d; j++)
    
20.         {
    
21.             x=(unsigned long long)r[j] * i+c;
    
22.             r[j]=x%b;
    
23.             c=x/b;
    
24.         }
    
25.         if (c!=0)
    
26.         {
    
27.             d++;
    
28.             r[d]=c;
    
29.         }
    
30.     }
    
31. }
    
32. 
    
33. 
    
34. int
    
35. main()
    
36. {
    
37.     int i;
    
38.     factorial(1000, r);
    
39. 
    
40.     printf ("%ld", r[285]);
    
41.     for (i=284; i>=0; i--)
    
42.         printf ("%09ld", r[i]);
    
43. 
    
44. }
    

复制代码

402387260077093773543702433923003985......410970027753472000000000000000000000......0

real        0m0.019s
user        0m0.012s
sys        0m0.000s