OFFSET
0,20
COMMENTS
Palindromes themselves are not 'Reverse and Add!'ed, so they yield a zero!
Numbers n that may have a(n) = -1 (i.e., potential Lychrel numbers) appear in A023108. - Michael De Vlieger, Jan 11 2018
REFERENCES
D. Wells, The Penguin Dictionary of Curious and Interesting Numbers Penguin Books, 1987, pp. 142-143.
LINKS
Kerry Mitchell, Table of n, a(n) for n = 0..195
P. De Geest, Some thematic websources
Jason Doucette, World Records
S. K. Eddins, The Palindromic Order Of A Number [archived page]
Kerry Mitchell, Table of n, a(n) for n = 0..10000 (The -1 entries are only conjectural)
Kerry Mitchell, Table of n, a(n) for n = 0..100000 (The -1 entries are only conjectural)
I. Peter, Search for the biggest numeric palindrome [archived page]
T. Trotter, Jr., Palindrome Power [archived page]
Eric Weisstein's World of Mathematics, 196-Algorithm.
EXAMPLE
19 -> 19+91 = 110 -> 110+011 = 121 = palindrome, took 2 steps, so a(19)=2.
n = 89 needs 24 steps to end up with the palindrome 8813200023188. See A240510. - Wolfdieter Lang, Jan 12 2018
MATHEMATICA
rev[n_]:=FromDigits[Reverse[IntegerDigits[n]]]; radd[n_]:=n+rev[n];
pal[n_]:=If[n==rev[n], True, False];
raddN[n_]:=Length[NestWhileList[radd[#]&, n, pal[#]==False&]]-1;
raddN/@Range[0, 195] (* Ivan N. Ianakiev, Aug 31 2015 *)
With[{nn = 10^3}, Array[-1 + Length@ NestWhileList[# + IntegerReverse@ # &, #, !PalindromeQ@ # &, 1, nn] /. k_ /; k == nn -> -1 &, 200]] (* Michael De Vlieger, Jan 11 2018 *)
PROG
(PARI) rev(n)={d=digits(n); p=""; for(i=1, #d, p=concat(Str(d[i]), p)); return(eval(p))}
a(n)=if(n==rev(n), return(0)); for(k=1, 10^3, i=n+rev(n); if(rev(i)==i, return(k)); n=i)
n=0; while(n<100, print1(a(n), ", "); n++) \\ Derek Orr, Jul 28 2014
(PARI) A033665(n, LIM=333)={-!for(i=0, LIM, my(r=A004086(n)); n==r&&return(i); n+=r)} \\ with {A004086(n)=fromdigits(Vecrev(digits(n)))}. The second optional arg is a search limit that could be taken smaller up to very large n, e.g., 99 for n < 10^9, 200 for n < 10^14, 250 for n < 10^18: see A065199 for the records and A065198 for the n's. - M. F. Hasler, Apr 13 2019, edited Feb 16 2020
(Python) A033665 = lambda n, LIM=333: next((i for i in range(LIM) if is_A002113(n) or not(n := A004086(n)+n)), -1) # The second, optional argument is a search limit, see above. - M. F. Hasler, May 23 2024
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
EXTENSIONS
More terms from Patrick De Geest, Jun 15 1998
I truncated the b-file at n=195, since the value of a(196) is not presently known (cf. A006960). The old b-files are now a-files. - N. J. A. Sloane, May 09 2015
STATUS
approved