>>516の続き

変数を動かして総和をとるシグマ記号を
x, a, b, c の4重の繰り返しに書き直すと
P(C_m<=N)
=∑[x=0, min(n, N-3m)](
∑[a=m, x+m](
∑[b=m, x+(a-m)+m](
∑[c=m, x+(a-m)+(b-m)+m](
  ((a+b+c)!n!/(a!b!c!x!(n-x)!))
  *((9^x)/((3^(a+b+c))(10^n))
))))

これをC言語など任意のプログラムで表せば
中央値、達成率90%・95%・99%に
対応する回数の計算や
ヒストグラムの描画などができるはず