A385365 Consecutive states of the linear congruential pseudo-random number generator (3661*s + 30809) mod 145800 when started at s=1.

1, 34470, 108479, 13228, 52917, 137546, 139315, 54624, 117473, 135262, 88191, 96860, 49669, 56418, 124307, 76936, 7905, 102614, 119863, 137052, 80381, 81250, 55059, 106208, 9697, 102126, 82895, 99604, 35253, 59042, 107971, 48840, 83249, 83398, 45687, 58316

Offset: 1

Author

Sean A. Irvine, Jun 26 2025

Keywords

• nonn
• look
• easy

Comments

Periodic with period 145800.

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

Crossrefs

Cf. A383127, A385361.

Programs

Maple

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

Mathematica

NestList[Mod[3661*# + 30809, 145800] &, 1, 50] (* after Paolo Xausa *)

Formula

a(n) = (3661 * a(n-1) + 30809) mod 145800.