Add efficient range function for integral types

[Shihab-Shahriar]

Code mostly taken from here, credits to Lemire.

For range N and floating point types, we use the basic rand() * N formula in the uniform method.

For integral type, uniform by default uses biased version. Users can use non-biased verison by directly invoking range.