OFFSET
1,2
COMMENTS
Periodic with period 2^46.
Used for Monte Carlo simulations.
LINKS
Sean A. Irvine, Table of n, a(n) for n = 1..10000
Forrest B. Brown and Yasunobu Nagaya, The MCNP5 Random Number Generator, LA-UR-02-3782, American Nuclear Society, Winter Meeting, 2002.
George S. Fishman, Multiplicative Congruential Random Number Generators with Modulus 2^beta: An Exhaustive Analysis for beta = 32 and a Partial Analysis for beta = 48, Math. Comp., 54, 189 (1990), 331-344.
B. D. Ripley, Thoughts on pseudorandom number generators, J of Computational and Applied Mathematics, 31, 1 (1990), 153-163.
FORMULA
a(n) = 5^19 * a(n-1) mod 2^48.
MAPLE
a:= proc(n) option remember; `if`(n<2, n,
irem(5^19*a(n-1), 2^48))
end:
seq(a(n), n=1..16); # Alois P. Heinz, Jun 09 2025
MATHEMATICA
NestList[Mod[5^19*#, 2^48] &, 1, 30] (* Paolo Xausa, Jun 11 2025 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Sean A. Irvine, Jun 09 2025
STATUS
approved