OFFSET
1,2
COMMENTS
Periodic with period 2^31 (Dyck et al. mistakenly give the period as 2^29).
Proposed by Dyck et al. for FORTRAN 77 on VAX or IBM computers.
REFERENCES
V. A. Dyck, J. D. Lawson, and J. A. Smith, FORTRAN 77: An Introduction to Structured Problem Solving, Reston Pub. Co., 1984 (see p. 467).
LINKS
FORMULA
a(n) = (843314861 * a(n-1) + 453816693) mod 2^31.
MAPLE
a:= proc(n) option remember; `if`(n<2, n,
irem(843314861*a(n-1)+453816693, 2^31))
end:
seq(a(n), n=1..23); # Alois P. Heinz, May 29 2025
MATHEMATICA
NestList[Mod[843314861*# + 453816693, 2^31] &, 1, 50] (* Paolo Xausa, May 30 2025 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Sean A. Irvine, May 29 2025
STATUS
approved