cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A384696 Consecutive states of the linear congruential pseudo-random number generator Cray RANF when started at 1.

This page as a plain text file.
%I A384696 #19 Jun 12 2025 21:54:59
%S A384696 1,44485709377909,232253848878969,94800993741645,243522309605169,
%T A384696 20783065360997,154093299791145,161954398135485,183663036741473,
%U A384696 207319719370837,142356556532697,278312552510253,242082341486737,37630394630981,176334633251721,233894773868189
%N A384696 Consecutive states of the linear congruential pseudo-random number generator Cray RANF when started at 1.
%C A384696 Periodic with period 2^46.
%C A384696 Also implemented in the GNU Scientific Library as gsl_rng_ranf.
%C A384696 Also implemented in PASCLIB for Pascal programs on the CDC Cyber computers.
%H A384696 Sean A. Irvine, <a href="/A384696/b384696.txt">Table of n, a(n) for n = 1..10000</a>
%H A384696 George S. Fishman, <a href="https://doi.org/10.2307/2008698">Multiplicative Congruential Random Number Generators with Modulus 2^beta: An Exhaustive Analysis for beta = 32 and a Partial Analysis for beta = 48</a>, Math. Comp., 54, 189 (1990), 331-344.
%H A384696 Mark Galassi, Jim Davies, James Theiler, Brian Gough, Gerard Jungman, Michael Booth, and Fabrice Rossi, <a href="https://www.gnu.org/software/gsl/doc/latex/gsl-ref.pdf">GNU Scientific Library: Reference Manual</a>, 2003 (see p. 181).
%H A384696 B. D. Ripley, <a href="https://doi.org/10.2307/1402590">Computer Generation of Random Variables: A Tutorial</a>, International Statistical Review, 51 (1983), 301-309.
%H A384696 <a href="/index/Ps#PRN">Index entries for sequences related to pseudo-random numbers</a>.
%F A384696 a(n) = 44485709377909 * a(n-1) mod 2^48.
%p A384696 a:= proc(n) option remember; `if`(n<2, n,
%p A384696       irem(44485709377909*a(n-1), 2^48))
%p A384696     end:
%p A384696 seq(a(n), n=1..16);  # _Alois P. Heinz_, Jun 09 2025
%t A384696 NestList[Mod[44485709377909*#, 2^48] &, 1, 30] (* _Paolo Xausa_, Jun 11 2025 *)
%Y A384696 Cf. A384746, A384775, A384776, A384778, A384779, A384780.
%Y A384696 Cf. A384546, A384547, A384548, A384549, A384550, A384551, A384552
%Y A384696 Cf. A096550-A096561 other pseudo-random number generators.
%K A384696 nonn,easy
%O A384696 1,2
%A A384696 _Sean A. Irvine_, Jun 07 2025