pub fn random_permutation(n: usize, seed: Option<u64>) -> Vec<usize>
Generate a random permutation of indices