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