A385341 Consecutive states of the linear congruential pseudo-random number generator (1021*s + 25673) mod 121500 when started at s=1.

1, 26694, 64247, 11860, 106233, 111566, 89059, 72912, 110825, 61498, 121131, 13424, 2077, 80790, 13763, 105196, 24789, 63242, 79255, 26028, 113261, 118654, 35907, 115220, 53293, 5826, 20519, 77572, 8685, 23558, 21391, 117384, 75737, 79150, 40323, 6956, 80749

Offset: 1

Author

Sean A. Irvine, Jun 25 2025

Keywords

• nonn
• easy

Comments

Periodic with period 121500.

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

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 121500.

Crossrefs

Cf. A383126, A385340.

Programs

Maple

a:= proc(n) option remember; `if`(n<2, n,
      irem(1021*a(n-1)+25673, 121500))
    end:
seq(a(n), n=1..44);  # after Alois P. Heinz

Mathematica

NestList[Mod[1021*# + 25673, 121500] &, 1, 50] (* after Paolo Xausa *)

Formula

a(n) = (1021 * a(n-1) + 25673) mod 121500.