RU_N
j = arc4random_uniform(RU_N);
j = (j + 1) % RU_N;
ru_g = pmod(RU_GEN, j, RU_N);
r = permute15(ru_seed ^ pmod(ru_g, ru_seed2 + ru_x, RU_N)) | ru_msb;