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

1, 27714, 92447, 20500, 70173, 266, 82579, 48912, 103565, 113338, 12531, 86444, 39877, 9630, 119003, 32296, 88809, 6842, 17695, 55668, 41561, 44674, 80307, 28760, 40333, 89826, 17039, 53272, 11325, 54998, 10591, 14904, 69737, 82390, 27663, 109856, 74269, 98202

Offset: 1

Author

Sean A. Irvine, Jun 26 2025

Keywords

• nonn
• look
• 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, A385341.

Programs

Magma

[n le 1 select 1 else (2041*Self(n-1)+ 25673) mod 121500: n in [1..40]]; // Vincenzo Librandi, Jun 27 2025

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 *)

Formula

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