login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A096551 Consecutive internal states of a linear congruential pseudo-random number generator with a parameter proposed by George Marsaglia as a "candidate for the best of all multipliers". 3
1, 69069, 475559465, 2801775573, 1790562961, 3104832285, 4238970681, 2135332261, 381957665, 1744831853, 1303896393, 1945705589, 2707602097, 4198202557, 3820321881, 201201733, 2583294017, 4003049741, 2417848425 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
REFERENCES
D. E. Knuth, The Art of Computer Programming Third Edition. Vol. 2 Seminumerical Algorithms. Chapter 3.3.4 The Spectral Test, Page 108. Addison-Wesley 1997.
Marsaglia, G., The structure of linear congruential sequences, in Applications of Number Theory to Numerical Analysis, (edited by S. K. Zaremba), Academic Press, New York, 249-286, 1972
LINKS
FORMULA
a(1)=1, a(n) = 69069 * a(n-1) mod 2^32. The sequence is periodic with period length 2^30. - corrected by Hugo Pfoertner, Aug 10 2011
MAPLE
a:= proc(n) option remember; `if`(n<2, n,
irem(69069 *a(n-1), 4294967296))
end:
seq(a(n), n=1..30); # Alois P. Heinz, Jun 10 2014
PROG
(PARI) a(n)=lift(Mod(69069, 2^32)^(n-1)) \\ Charles R Greathouse IV, Jan 14 2016
CROSSREFS
Cf. A096550-A096561 for other pseudo-random number generators.
Sequence in context: A199997 A250500 A354194 * A096552 A031663 A105010
KEYWORD
nonn,easy
AUTHOR
Hugo Pfoertner, Jul 18 2004
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 05:48 EDT 2024. Contains 371265 sequences. (Running on oeis4.)