pub fn get_multibranch_close_score(
    seq: &[usize],
    pos_pair_close: &(usize, usize)
) -> f32