18/11/13 15:39:05.37 UTrDDuDZ.net
シミュレーション結果と一致したので答には自信があったのだが、数学板でのレスと一致したので確信がもてた。
Rのスクリプトはこれ
Nmin=g+s-x
sub <- function(){
n=sample(Nmin:Nmax,1)
a=sample(n,g)
b=sample(n,s)
if(sum(a %in% b)==x) return(n)
}
sim <- function(){
r=sub()
while(is.null(r)) r=sub()
return(r)
}
re=replicate(1e5,sim())
summary(re) ; hist(re,freq=F,col="lightblue",main='',xlab='N')