A383956 - OEIS

%I #21 May 25 2025 20:48:31

%S 1,7771826,12906479,12752200,14370573,4177230,16102619,5888068,

%T 8967385,14199722,1838727,7559424,14513509,9092550,15771891,2282364,

%U 11580593,15929250,14479391,2474936,6872765,1998142,6754315,6251956,4652937,6660762,6157495,1357168

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

%C The sequence is periodic with period 2^24.

%C 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.

%H Sean A. Irvine, <a href="/A383956/b383956.txt">Table of n, a(n) for n = 1..10000</a>

%H Andrew Trotman, <a href="https://www.cs.otago.ac.nz/homepages/andrew/poly/Poly.htm">The Poly Preservation Project</a>

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Poly-1">Poly-1</a>

%H <a href="/index/Ps#PRN">Index entries for sequences related to pseudo-random numbers.</a>

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

%p a:= proc(n) option remember; `if`(n<2, n,

%p       irem(4253261*a(n-1)+3518565, 2^24))

%p     end:

%p seq(a(n), n=1..28);  # _Alois P. Heinz_, May 21 2025

%t NestList[Mod[4253261*# + 3518565, 2^24] &, 1, 50] (* _Paolo Xausa_, May 22 2025 *)

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

%K nonn,easy

%O 1,2

%A _Sean A. Irvine_, May 15 2025