A088753 Numbers k such that the Reverse and Add! trajectory of k (presumably) does not reach a palindrome (with the exception of k itself) and does not join the trajectory of any term m < k.
196, 879, 1997, 7059, 9999, 10553, 10563, 10577, 10583, 10585, 10638, 10663, 10668, 10697, 10715, 10728, 10735, 10746, 10748, 10783, 10785, 10787, 10788, 10877, 10883, 10963, 10965, 10969, 10977, 10983, 10985, 12797, 12898, 13097, 13197, 13694, 14096, 14698, 15297, 15597, 18598, 18798
Offset: 1
Examples
From _M. F. Hasler_, Apr 13 2019: (Start) All numbers < 196 quickly reach a palindrome under iterations of the reverse-and-add function A056964, cf. A033665. a(1) = 196 is the smallest integer which appears to never reach a palindrome (checked up to 10^9 iterations!). Next, A056964(196) = 196 + 691 = 887 is in the orbit of 196 and will therefore never reach a palindrome if 196 does not. However, we do not list this term in this sequence because it is in the orbit of the smaller term 196. Similarly, 295 + 592 = 887 = A056964(196). Therefore, 295 will also never reach a palindrome if 196 (and therefore 887) doesn't. But again we will not list this number, because its orbit merges into that of the smaller term 196. The next number which appears to be a Lychrel and has an orbit (conjectured to be) disjoint with that of 196 is 897 = a(2). (End)
Links
- M. F. Hasler, Table of n, a(n) for n = 1..74 (all terms up to 10^5), Apr 13 2019.
- W. VanLandingham, 196 and Other Lychrel Numbers [Copy on web.archive.org as of 05/2018, original site p196.org seems no longer available. - _M. F. Hasler_, Apr 13 2019]
Crossrefs
Programs
-
Mathematica
limit = 10^3; (* Assumes that there is no palindrome if none is found before "limit" iterations *) utraj = {}; Select[Range[0, 20000], (np = # + IntegerReverse[#]; x = NestWhileList[ # + IntegerReverse[#] &, np, ! PalindromeQ[#] &, 1, limit]; If[Length[x] >= limit && Intersection[x, utraj] == {}, utraj = Union[utraj, x, {np}]; True, utraj = Union[utraj, x, {np}]]) &] (* Robert Price, Oct 16 2019 *)
-
PARI
A088753_upto(LIM=2e4,M=1e199)={my(U=[],a=List());for(n=1,LIM, my(t=n); while( t
A002113(t=A056964(t)) && next(2)); setsearch(U,t) && next; U=setunion(U,[t]); print1(n","); listput(a,n)); Set(a)} \\ M. F. Hasler, Apr 13 2019
Extensions
Edited by M. F. Hasler, Apr 13 2019
Comments