OFFSET
1,1
COMMENTS
This is a "subtract with carry" generator and is used to define the ranlux24 (A221561) generator.
Period is (2^572 - 2^236)/3. - Charles R Greathouse IV, Oct 28 2022
LINKS
Eric M. Schmidt, Table of n, a(n) for n = 1..1000
George Marsaglia and Arif Zaman, A New Class of Random Number Generators, Annals of Applied Probability, Volume 1, Number 3 (1991), 462-480.
Alexei Sibidanov, A revision of the subtract-with-borrow random number generators, arXiv preprint (2017). arXiv:1705.03123 [physics.comp-ph]
PROG
(C++)
#include <iostream>
#include <random>
void A221559(int max)
{
std::ranlux24_base gen;
for (int i = 1; i <= max; ++i)
std::cout << i << ' ' << gen() << '\n';
}
CROSSREFS
KEYWORD
nonn
AUTHOR
Eric M. Schmidt, Jan 20 2013
STATUS
approved