OFFSET
1,2
COMMENTS
Periodic with period 2^46.
LINKS
Sean A. Irvine, Table of n, a(n) for n = 1..10000
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.
FORMULA
a(n) = 43272750451645 * a(n-1) mod 2^48.
MAPLE
a:= proc(n) option remember; `if`(n<2, n,
irem(43272750451645*a(n-1), 2^48))
end:
seq(a(n), n=1..16); # Alois P. Heinz, Jun 09 2025
MATHEMATICA
NestList[Mod[43272750451645*#, 2^48] &, 1, 30] (* Paolo Xausa, Jun 11 2025 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Sean A. Irvine, Jun 09 2025
STATUS
approved
