paralell_pairwise_discrete_rcpp function

Wrapper function to discrete-based pairwise two-point estimation in C++