https://twitter.com/kfuku0502/status/1234498687608512512

fukushima.R

library(tidyverse)
library(furrr)
library(future)
future::plan("multicore")

L = 1000
M = 100
N = 200
w = rgamma(L, 1)
prob = w / sum(w)

simulate = function(L, M, N, prob) {
  stones1 = sample(seq_len(L), M, replace = FALSE, prob = prob)
  stones2 = sample(seq_len(L), N, replace = FALSE, prob = prob)
  sum(stones1 %in% stones2)
}

x = furrr::future_map_int(rep(L, 10000), simulate, M = M, N = N, prob = prob)
tibble::tibble(x) %>%
  ggplot(aes(x)) +
  geom_bar()