implemented exact p-value calculation + parallelization + zeta random variate for the discrete case