你这个问题有问题
丢10000次，如何做到每个面各出现一次
应该是每个面至少出现一次吧

无论你想表达的意思是什么都不为 0。
就算是荒诞的每一面各出现一次，一共 6 次都有可能，因为骰子有可能出现各种意外比如摔的粉碎，摩擦炸裂，甚至太阳爆炸

下个r跑个几百万次，几行代码就搞定了。

先考虑十面骰掷100次的简化版问题
P（至少全掷出一次）=1-P（至少有一种没掷出来）
记P(n)为刚好有n种没掷出来的情况数，
则P（至少有一种没掷出来）
=ΣP（n）÷10^100
P（1）=C(1,10)×1^100=10
记Q（1）=1^100,即只出现一种的情况数
P（2）
=C(2,10)×[2^100-C(1,2)×1^100]
=C(2,10)×[2^100-C(1,2)×Q（1）]
记Q（2）=2^100-C(1,2)×Q（1），即只出现两种的情况数.
P（3）=C(3,10)×[3^100-C(2,3)×Q（2）-C(1,3)×Q（1）]
记Q（3）=3^100-C(2,3)×Q（2）-C(1,3)×Q（1）,即只出现三种的情况数.
观察可得，Q（n）
=n^100-C(n-1,n)×Q（n-1）-C(n-2,n)×Q（n-2）-…-
C(1,n)×Q（1）
P（n）=C(n,10)×Q（n）
已知Q（1）和Q（2），运用递归可以求出Q（n），则P（1）…P（10）均可以求得。写个程序很快就可以出结果。
至于手算我就不算了，化简这块掌握的不好。

C10000，100 ×A100，100

还是1/100

容斥原理

@哦吼吼看的哦～～
信誓旦旦地说百面骰不在Persi Diaconis列出的30类fair dice之中的某人

看到"two infinite classes of prism and anti-prisms"了吗？棱柱与反棱柱，作为两类按顶点对称的多面体，按下一段所述取对偶，得到的双棱锥与偏方面体（trapezohedron，又称双反角锥、鸢形多面体），是fair polyhedra种的两类。取底面为正50边形的棱柱/反棱柱得到的就是100面的双棱锥/双反角锥。
哪怕你没看过他的论文原文，只看过Numberphile上对他的采访的那两期节目，在列出30类Fair Dice的前后也明确地举了双五棱锥做例子并且提到了它是如何推广到任意正多边形的。

我来个随机过程的解法吧。这个可以看作一个随机游走，每个节点的取值是0到100的整数，代表当前已经出现过的面数，每个节点i有(100-i)/100的概率跳到下一个节点，i/100的概率留在当前节点，重复10000次相当于求转移矩阵的10000次方，再跟初始状态相乘。几年前学的了，应该没错。

P≈100