cp's OEIS Frontend

A384746 Consecutive states of the linear congruential pseudo-random number generator MCNP from Los Alamos when started at 1.

%I A384746 #15 Jun 11 2025 10:11:58
%S A384746 1,19073486328125,29763723208841,187205367447973,131230026111313,
%T A384746 264374031214925,74735272014937,31978779697717,72377397341089,
%U A384746 127824407320157,39323977335081,168134765887429,73951303845617,27971537168493,266449281326841,41546074810965
%N A384746 Consecutive states of the linear congruential pseudo-random number generator MCNP from Los Alamos when started at 1.
%C A384746 Periodic with period 2^46.
%C A384746 Used for Monte Carlo simulations.
%H A384746 Sean A. Irvine, <a href="/A384746/b384746.txt">Table of n, a(n) for n = 1..10000</a>
%H A384746 Forrest B. Brown and Yasunobu Nagaya, <a href="https://www.osti.gov/servlets/purl/976209">The MCNP5 Random Number Generator</a>, LA-UR-02-3782, American Nuclear Society, Winter Meeting, 2002.
%H A384746 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 A384746 B. D. Ripley, <a href="https://doi.org/10.1016/0377-0427(90)90346-2">Thoughts on pseudorandom number generators</a>, J of Computational and Applied Mathematics, 31, 1 (1990), 153-163.
%H A384746 <a href="/index/Ps#PRN">Index entries for sequences related to pseudo-random numbers</a>.
%F A384746 a(n) = 5^19 * a(n-1) mod 2^48.
%p A384746 a:= proc(n) option remember; `if`(n<2, n,
%p A384746       irem(5^19*a(n-1), 2^48))
%p A384746     end:
%p A384746 seq(a(n), n=1..16);  # _Alois P. Heinz_, Jun 09 2025
%t A384746 NestList[Mod[5^19*#, 2^48] &, 1, 30] (* _Paolo Xausa_, Jun 11 2025 *)
%Y A384746 Cf. A096550-A096561 other pseudo-random number generators.
%K A384746 nonn,easy
%O A384746 1,2
%A A384746 _Sean A. Irvine_, Jun 09 2025