Speed up and implement both version for discrete distributions

1. Change call signature to have array of arrays and optional keyword to
   specify which version of the k-sample AD test should be computed.

2. Get rid of all inner loops and list comprehensions by using
   np.searchsorted.

3. Both version given by Scholz & Stephens for discrete sample can be
   computed now.