A383956 Consecutive states of the linear congruential pseudo-random number generator used by BASIC on the Poly-1 computer when started at 1.

1, 7771826, 12906479, 12752200, 14370573, 4177230, 16102619, 5888068, 8967385, 14199722, 1838727, 7559424, 14513509, 9092550, 15771891, 2282364, 11580593, 15929250, 14479391, 2474936, 6872765, 1998142, 6754315, 6251956, 4652937, 6660762, 6157495, 1357168

Offset: 1

Author

Sean A. Irvine, May 15 2025

Keywords

• nonn
• easy

Comments

The sequence is periodic with period 2^24.

In Poly-1 BASIC, random numbers were generated with the RND(m) function. With RND(0) the internal state was returned as a floating-point number state/2^24, otherwise the state was return modulo m.

Links

• Sean A. Irvine, Table of n, a(n) for n = 1..10000
• Andrew Trotman, The Poly Preservation Project
• Wikipedia, Poly-1
• Index entries for sequences related to pseudo-random numbers.

Crossrefs

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

Programs

Maple

a:= proc(n) option remember; `if`(n<2, n,
      irem(4253261*a(n-1)+3518565, 2^24))
    end:
seq(a(n), n=1..28);  # Alois P. Heinz, May 21 2025

Mathematica

NestList[Mod[4253261*# + 3518565, 2^24] &, 1, 50] (* Paolo Xausa, May 22 2025 *)

Formula

a(n) = (4253261 * a(n-1) + 3518565) mod 2^24.