OFFSET
1,2
COMMENTS
Periodic with period 145800.
Described in Numerical Recipes as a "quick and dirty" generator.
LINKS
FORMULA
a(n) = (3661 * a(n-1) + 30809) mod 145800.
MAPLE
a:= proc(n) option remember; `if`(n<2, n,
irem(3661*a(n-1)+30809, 145800))
end:
seq(a(n), n=1..44); # after Alois P. Heinz
MATHEMATICA
NestList[Mod[3661*# + 30809, 145800] &, 1, 50] (* after Paolo Xausa *)
CROSSREFS
KEYWORD
AUTHOR
Sean A. Irvine, Jun 26 2025
STATUS
approved
