Recaman numbers
John Layman
layman at calvin.math.vt.edu
Thu Oct 4 17:48:31 CEST 2001
Following the recent flurry of notes on the Recaman numbers (A005132),
I finally couldn't resist the temptation to get involved. I have completed
the calculation of R(n) to beyond n=2 billion and a calculation proceeding
as I write will push this to 15 billion. The much sought-after R=1355 does
not occur by n=2 billion, the seven smallest values that have not occurred
by then being 1355, 1643, 1814, 2405, 2406, 2520, and 2991.
An interesting aspect of the Recamans that has not been previously noted is
the fact that the successive differences of the integer part of the ratio
R(n)/n is almost always negative, with only 36 occurrences where
FractPart(R(n+1)/(n+1))-FractPart(R(n)/n) is positive up to n=2*10^9, these
being at n=2,4,7,12,22,40,77,135,249,454,845,1521,2753,5046,9318,17224,31222,
57072,99742,181694,328256,589933,1034839,1788538,3225919,5784586,10212211,
18399785,32148795,58056876,101769230,173395920,302890749,511561221,904036925,
1610187039. The ratio of successive numbers in this sequence is usually about
1.7 to 1.8.
Neil, my calculations use routines which I (pardon the expression!) "hacked", that allow 15 boolean,
or bit, values to be set or read per 2-byte INTEGER*2 Fortran data value. This allows a 1-dim
array of 10^9 INTEGER*2 values to suffice for calculating R(n) up to n=15*10^9. My server Fortran
permits array indexes only up to 2^30 or just over 1 billion.
More information about the SeqFan
mailing list