cp's OEIS Frontend

A385358 Consecutive states of the linear congruential pseudo-random number generator (741*s + 66037) mod 312500 when started at s=1.

%I A385358 #10 Jul 06 2025 18:21:27
%S A385358 1,66778,173535,217972,20789,158186,94363,301520,54857,90074,248371,
%T A385358 46448,109005,213742,11359,45556,73033,120990,32127,122144,262241,
%U A385358 11618,237475,97512,134929,48426,12203,45960,59897,74714,116611,224788,71445,194282,278999
%N A385358 Consecutive states of the linear congruential pseudo-random number generator (741*s + 66037) mod 312500 when started at s=1.
%C A385358 Periodic with period 312500.
%C A385358 Described in Numerical Recipes as a "quick and dirty" generator.
%H A385358 Sean A. Irvine, <a href="/A385358/b385358.txt">Table of n, a(n) for n = 1..10000</a>
%H A385358 <a href="/index/Ps#PRN">Index entries for sequences related to pseudo-random numbers.</a>
%H A385358 <a href="/index/Rec#order_312500">Index entries for linear recurrences with constant coefficients</a>, order 312500.
%F A385358 a(n) = (741 * a(n-1) + 66037) mod 312500.
%p A385358 a:= proc(n) option remember; `if`(n<2, n,
%p A385358       irem(741*a(n-1)+66037, 312500))
%p A385358     end:
%p A385358 seq(a(n), n=1..44);  # after _Alois P. Heinz_
%t A385358 NestList[Mod[741*# + 66037, 312500] &, 1, 50] (* after _Paolo Xausa_ *)
%Y A385358 Cf. A383127, A385357.
%K A385358 nonn,easy
%O A385358 1,2
%A A385358 _Sean A. Irvine_, Jun 26 2025