A383126 Consecutive internal states of the linear congruential pseudo-random number generator (281*s + 28411) mod 134456 when started at 1.
1, 28692, 23503, 44410, 3213, 124528, 62219, 32670, 65673, 62052, 120199, 55874, 132109, 41184, 37899, 56006, 34745, 110924, 4263, 16210, 11917, 15688, 134147, 76038, 16585, 117292, 45743, 108874, 100493, 31184, 51475, 106094, 126049, 86252, 63143, 23402
Offset: 1
References
- William H. Press, Saul A. Teukolsky, William T. Vetterling, and Brian P. Flannery, Numerical Recipes in C (2nd ed), Cambridge University Press, 1999 (see p. 285).
Links
- Sean A. Irvine, Table of n, a(n) for n = 1..10000
- W. E. Sharp and Carter Bays, A review of portable random number generators, Computers and Geosciences, 18, 1 (1982), 79-87.
- Index entries for sequences related to pseudo-random numbers.
- Index entries for linear recurrences with constant coefficients, order 134456.
Programs
-
Maple
a:= proc(n) option remember; `if`(n<2, n, irem(281*a(n-1)+28411, 134456)) end: seq(a(n), n=1..45); # after Alois P. Heinz
-
Mathematica
NestList[Mod[281*# + 28411, 134456] &, 1, 50] (* Paolo Xausa, Jun 18 2025 *)
Formula
a(n) = (281*a(n-1) + 28411) mod 134456.
Comments