pub fn pflip(weights: &[f64], n: usize, rng: &mut impl Rng) -> Vec<usize>
Draw n indices in proportion to their weights
n
weights