>>175
可読性を度外視して1行に纏めた。自分でも解読するのに難渋しそう。


* 血液型頻度比はA:O:B:AB＝4:3:2:1, n 人に m 種類の血液型が存在する確率 *)
solve[n_,m_] := Total[Binomial[n,#[[1]]] (1/10)^#[[1]] Binomial[n-#[[1]],#[[2]]] (2/10)^#[[2]] Binomial[n-#[[1]]-#[[2]],#[[3]]] (3/10)^#[[3]] (4/10)^#[[4]]& /@ (Select[Select[Partition[Flatten@Table[{x1,x2,x3,n-x1-x2-x3},{x1,Range[0,n]},{x2,Range[0,n]},{x3,Range[0,n]}],4],#[[4]]>=0&],Total@Boole@Table[#[[i]]>0,{i,Range[4]}] == m &])]

実行結果

In[1]:= solve[n_,m_] := Total[Binomial[n,#[[1]]] (1/10)^#[[1]] Binomial[n-#[[1]],#[[2]]] (2/10)^#[[2]] Binomial[n-#[[1]]-#[[2]],#[[3]]] (3/10)^#[[3]] (4/10)^#[[4]]& /@ (Select[Select[Partition[Flatten@Table[{x1,x2,x3,n-x1-x2-x3},{x1,Range[0,n]},{x2,Range[0,n]},{x3,Range[0,n]}],4],#[[4]]>=0&],Total@Boole@Table[#[[i]]>0,{i,Range[4]}] == m &])]

In[2]:= solve[40,3]

4660364912515525380760230058213856169
Out[2]= ---------------------------------------
312500000000000000000000000000000000000

尿瓶チンパフェチの嫁なしPhimoseくん（別名、自演認定厨）はWolfram言語はわからないみたいだな。
初学の俺と達人のコードが区別できないみたいだし。
嫁なしPhimoseくんが自演認定厨であることは度々、識者より指摘されている。
俺の想定解よりエレガントな解が投稿されていることを理解できないようだ。