この規模の数になると計算が終わらないな。

library(gmp)
m=18200 # 有名人の数(桜を見る会参加人数)
n=1.268e5　# 日本の人口
x=0:n # 感染者数:x, 非感染数:n-x
pmf=1- chooseZ(n-x,m)/chooseZ(n,m) # 1 - (ｍ人全員非感染の確率)
pdf=pmf/sum(pmf) # 確率密度関数化して
(E=sum(x*pdf)) # 期待値を計算
as.numeric(E) #