ru_n
j = arc4random() % p->ru_n;
j = (j + 1) % p->ru_n;
p->ru_g = pmod(p->ru_gen, j, p->ru_n);
return (p->ru_seed ^ pmod(p->ru_g, p->ru_seed2 + p->ru_x, p->ru_n)) |