A linear-time algorithm to sample the dual-birth model