OFFSET
1,2
COMMENTS
Periodic with period 2^47.
LINKS
Sean A. Irvine, Table of n, a(n) for n = 1..10000
Heidi G. Kuehn, A 48-bit pseudo-random generator, C ACM, 4, 8 (1961), 350-352.
FORMULA
a(n) = (513*a(n-1) + 29741096258473) mod 2^47.
MAPLE
a:= proc(n) option remember; `if`(n<2, n,
irem(513*a(n-1)+29741096258473, 2^47))
end:
seq(a(n), n=1..17); # Alois P. Heinz, Jun 09 2025
MATHEMATICA
NestList[Mod[513*# + 29741096258473, 2^47] &, 1, 20] (* Paolo Xausa, Jun 11 2025 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Sean A. Irvine, Jun 08 2025
STATUS
approved
