底辺私立医大を卒業した医者って頭悪いよね？ Part30

350 ：卵の名無しさん (ﾜｯﾁｮｲ 0324-dI8p [149.50.210.17 [上級国民]])：2024/03/12(火) 07:01:18.48 ID:zL8P6Hg/0.net

calc=function(n,frac=FALSE){
n2v=function(n){
nn=1:(n-5)
re=NULL
for(x1 in nn){
for(x2 in nn){
for(x3 in nn){
for(x4 in nn){
for(x5 in nn){
x6=n-x1-x2-x3-x4-x5
if(x6 %in% nn) re=rbind(re,c(x1,x2,x3,x4,x5,x6))
}}}}}
re
}
re=n2v(n)
nu=0
for(i in 1:nrow(re)){
v=re[i,]
cs=c(0,cumsum(v))
nu=nu+prod(choose(n-cs,v)[1:6])
# choose(n,v[1])*choose(n-v[1],v[2])*choose(n-v[1]-v[2],v[3])*
#choose(n-v[1]-v[2]-v[3],v[4])*choose(n-v[1]-v[2]-v[3]-v[4],v[5])*
#choose(n-v[1]-v[2]-v[3]-v[4]-v[5],v[6])
}
nu
de=6^n
# cat(nu,'/',de,'\n')
if(frac) print(MASS::fractions(nu/de))
nu/de
}

pn=diff(c(0,sapply(6:20,calc)))
plot(6:20,pn)