近似计算方法,同时因为这个数比较小,N相应要取得比较大,用了并行计算的包
library(future.apply)
plan(multisession)#To go back to sequential processing, use plan(sequential)
N = 10000000
tmp = future_replicate(N,length(intersect(sample(300,110),sample(300,110))))
mean(tmp/110 >= 0.5)
# 0.0002286