OFFSET
1,2
COMMENTS
Periodic with period 23^4.
LINKS
Sean A. Irvine, Table of n, a(n) for n = 1..10000
J. Eichenauer-Herrmann and H. Grothe, Upper bounds for the Beyer ratios of linear congruential generators, J Comp. Appl. Math., 31 1 (1990), 73-80.
W. E. Sharp and Carter Bays, A review of portable random number generators, Computers and Geosciences, 18, 1 (1982), 79-87.
FORMULA
a(n) = (7200*a(n-1) + 1) mod 23^4.
MAPLE
a:= proc(n) option remember; `if`(n<2, n,
irem(7200*a(n-1)+1, 279841))
end:
seq(a(n), n=1..45); # after Alois P. Heinz
MATHEMATICA
NestList[Mod[7200*# + 1, 23^4] &, 1, 50] (* Paolo Xausa, Jun 18 2025 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Sean A. Irvine, Jun 13 2025
STATUS
approved
