OFFSET
1,2
COMMENTS
Periodic with period 2^31-2.
Also used by SAS and GPSS/PC.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000
Visual Numerics, IMSL Fortran Numerical Library: User's Guide, Stat Library, Version 7.0 (see p. 1430).
FORMULA
a(n) = 397204094 * a(n-1) mod (2^31-1).
MAPLE
a:= proc(n) option remember; `if`(n<2, n,
irem(397204094*a(n-1), 2^31-1))
end:
seq(a(n), n=1..23); # Alois P. Heinz, May 29 2025
MATHEMATICA
NestList[Mod[397204094*#, 2^31 - 1] &, 1, 50] (* Paolo Xausa, May 30 2025 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Sean A. Irvine, May 27 2025
STATUS
approved
