>>636
朝飯前にシミュレーションしてみました。

タクシー総数N、観察台数mともに１〜１００の一様分布とし
m <= NならN台からm台を選んで最大数を記録してそれが６０になる場合にNとmを出力させるシミュレーションをしてみました。

sim1 <- function(){
M=m=0 # m:観察台数 M:最大番号 (初期値0)
while(M!=60){　 # M=60になるまで繰り返す
N=sample(1:100,1) # タクシー総数Nを1 ~ 100から選ぶ
m=sample(1:100,1) # 観察する台数mを1 ~ 100から選ぶ
if(m<=N) M=max(sample(1:N,m)) # m<=NならN台からm台選択して最大値をMにいれる
}
return(c(N=N,m=m)) # タクシー総数Nと観察台数mを返す
}
出力はこんな感じ
> replicate(10,sim1())
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
N 70 60 69 60 60 81 60 61 60 60
m 19 22 5 32 50 5 41 7 27 41

1万回繰り返してグラフにすると
re=replicate(1e4,sim1())
BEST::plotPost(re['N',],xlab='総数',cex=1,cex.lab=1,xlim=c(1,100))
BEST::plotPost(re['m',],xlab='観察数',cex=1,cex.lab=1,xlim=c(1,100))
BEST::plotPost(re['N',][re['m',]==5],breaks=20,xlab='観察数=5の時の総台数',cex.lab=1)


https://i.imgur.com/f9cLkZV.png
https://i.imgur.com/FqP1N6H.png

1万回の出力のうち観察台数が５台であったものを選んで
そのときの総台数をグラフにすると
https://i.imgur.com/xX3ZEaJ.png
これは>601に相当する結果。