趣味问题：100个人围成一圈相互射击，最后存活一人的概率

zgg___

6# zgg___
计算式子算是找到了，可似乎并不能明显的提高速度。
很好奇zgg是怎么算到300多的
wayne 发表于 2012-4-4 12:59

就是用的11层的那段呀，把第07行的m=30换成m=300就可以了，不过答案数比较大，比较费时间了，赫赫。

wayne

噢，我还以为你另外有数值计算的手段呢

wayne

21# mathe
mathe的这个递归式子效率要高一些。
这是按照动态规划的思路写的Mathematica代码：

1. f[0] = 0; f[1] = 1;
   
2. Table[u[n, 2] = 2^(n - 2); tmp = Table[u[n, k] = k^(n - k) (k - 1)^k - Sum[Binomial[k, h] u[n, k - h], {h, k - 2}]; Binomial[n, k] u[n, k], {k, 2, n}]; f[n] = tmp.f /@ Range[n - 2, 0, -1]/(n - 1)^n; {n, f[n]}, {n, 2, 100}]

复制代码