A385465 Consecutive internal states of the linear congruential pseudo-random number generator (4096*s + 150889) mod 714025 when started at 1.

1, 154985, 201224, 379543, 326592, 502896, 50780, 364494, 92038, 133337, 70116, 307975, 648339, 288458, 677507, 518411, 51995, 342959, 423778, 150802, 204256, 660190, 276454, 62823, 424897, 450076, 49635, 672749, 308318, 625217, 546071, 531405, 437569, 230763

Offset: 1

Author

Sean A. Irvine, Jun 29 2025

Keywords

• nonn
• easy

Comments

Periodic with period 714025.

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
• Index entries for sequences related to pseudo-random numbers.
• Index entries for linear recurrences with constant coefficients, order 714025.

Crossrefs

Cf. A385358, A385460.

Programs

Maple

a:= proc(n) option remember; `if`(n<2, n,
      irem(4096*a(n-1)+150889, 714025))
    end:
seq(a(n), n=1..45);  # after Alois P. Heinz

Mathematica

NestList[Mod[4096*# + 150889, 714025] &, 1, 50]

Formula

a(n) = (4096*a(n-1) + 150889) mod 714025.