18/11/06 17:20:32.86 jkx5i2bQ.net
ド底辺シリツ医大卒の裏口バカのカウントが終わりそうもないので、
シミュレーションのコードも書いてみた。
# simulation
mhs = function(x){ # maximum head sequence
y=paste(x,collapse='')
str="1"
if(!grepl(str,y)) return(0)
else{
while(grepl(str,y)){
str=paste0(str,"1")
}
return(nchar(str)-1)
}
}
(x=sample(0:1,100,rep=T)) ; mhs(x)
sim <- function(n){
mhs(sample(0:1,n,rep=T))
}
re=replicate(1e3,sim(10000))
hist(re,col='lightblue',main='',ann=F,freq=F)
summary(re)
which.max(tabulate(re))
max(tabulate(re))/length(re)