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

Author

Sean A. Irvine, Jun 17 2025

Keywords

• nonn
• look
• easy

Comments

Periodic with period 134456.

Described in Numerical Recipes as a "quick and dirty" generator.

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.

Crossrefs

Cf. A383127, A384431, A385002, A385003, A385037, A385039, A385078.

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.