这个题该用啥函数啊

n = 1000000; MC = {0, 0, 0, 0, 0}; Do[list = RandomInteger[5, 5];
MC[[Max[Count[list, #] & /@ Range[0, 5]]]] += 1, n];
MC/n // N

f[x_] := x + 1; N1 = 100000; Table[
Max[Nest[MapAt[f, #, RandomChoice[Range[6]]] &, Table[0, {i, 6}],
5]], {i, N1}];
Count[%, #] & /@ Range[1, 5]/N1 // N[#, 20] & //
DecimalForm[#, {2, 6}] &

Tally[#][[{5, 1}, 2]]/6^5 &@(Max /@
Nest[Flatten[Table[MapAt[# + 1 &, x,
#] & /@ Range[6], {x, #
}],
1] &, {Table[0, 6]}, 5])

m = 5;
n = 6;
caseList = IntegerPartitions[m, n];
f[case_] := Module[{},
Times @@((Times @@
Binomial[FoldList[
#1 - #
2 &,
#1, Most@#
2], #2]) & @璐村惂鐢ㄦ埛_053SRMK馃惥{{m,
case}, {n, Tally[case][[All, 2]]}})]
GroupBy[{Max@
#, f[#
]} & /@ caseList, First, Total[#[[All, 2]]] &]/n^m
由挡板法可知，m个球放入n个杯子最终可能的状态数是C(m+n-1, n-1)，其量级远远小于穷举的n^m，因此可以先用IntegerPartitions列出所有可能的状态数，再针对每种状态数计算对应的分布数，最后统计即可。这样速度比较快，即使m取到50也能1秒之内出结果

敲了半天代码回帖的时候才发现楼主设置了“关注的人才能回复”，考虑到这个回复还算有点意义，我迫不得已点了个关注发表了回复，但非常不爽。
我认为发帖者设置“关注的人才能回复”这种行为一方面阻碍了贴吧正常的交流，另一方面有很大的骗关注嫌疑，应当予以严厉打击，希望@xzcyr 吧主考虑在吧规中添加相应条款
我将在本条回复发布完后取关楼主，且不再回复任何设置了“关注的人才能回复”的帖子。

F = FactorialPower;
B = Binomial;
P[1] = F[6, 5]/6^5
P[2] = (B[5, 4] B[6, 2] B[4, 2] B[4, 1] + B[5, 2] B[6, 1] F[5, 3])/6^5
P[3] = B[5, 3] B[6, 1] 5^2/6^5
P[4] = B[5, 4] B[6, 1] B[5, 1]/6^5
P[5] = B[5, 5] B[6, 1]/6^5