A384315 Consecutive internal states of the linear congruential pseudo-random number generator for Demos started at 1.

1, 8192, 9317, 9225317, 19706942, 64858329, 25218022, 53630558, 40796927, 52681924, 53134651, 4299603, 62185148, 971592, 41535118, 60983366, 19606857, 50156573, 32119735, 27545333, 62690922, 51199833, 56863186, 18164438, 43380397, 13011312, 34587268, 44612022

Offset

1,2

Comments

Periodic with period 2^26-9318 = 67099546.

Originally defined by Downham and Roberts in 1967, subsequently used in the Demos modelling system in 2003.

Links

Sean A. Irvine, Table of n, a(n) for n = 1..10000

Graham Birtwistle, Demos: a system for Discrete Event Modelling on Simula, Jan 2003.

D. Y. Downham and F. D. K. Roberts, Multiplicative congruential pseudo-random number generators, The Computer J., 10, 1 (1967), 74-77

Index entries for sequences related to pseudo-random numbers.

Formula

a(n) = 2^13 * a(n-1) mod (2^26-9317).

Maple

a:= proc(n) option remember; `if`(n<2, n,
      irem(2^13*a(n-1), 2^26-9317))
    end:
seq(a(n), n=1..28);  # Alois P. Heinz, May 26 2025

Mathematica

NestList[Mod[2^13*#, 2^26 - 9317] &, 1, 50] (* Paolo Xausa, May 26 2025 *)

Crossrefs

Cf. A096550-A096561 other pseudo-random number generators.

Sequence in context: A289477 A222527 A035908 * A069274 A220585 A305756

Adjacent sequences: A384312 A384313 A384314 * A384316 A384317 A384318

Keyword

nonn,easy

Author

Sean A. Irvine, May 25 2025

Status

approved