OFFSET
1,1
COMMENTS
This sequence is believed to be finite. a(9) > 10^5, if it exists.
EXAMPLE
62 is in this sequence because 62 - (6+2)^2 = -2. Then -2 + (2)^2 = 2.
PROG
(PARI)
a(n)=k=n; c=1; v=List(); listput(v, k); while(c, if(k>=0, k-=sumdigits(k)^2; c+=1; if(k==2||k==3||k==0||k==6||k==9, return(k)); if(vecsearch(Vec(v), k), return(sumdigits(abs(k)))); listput(v, k)); if(k<0, k+=sumdigits(-k)^2; c+=1; if(k==2||k==3||k==0||k==6||k==9, return(k)); if(vecsearch(Vec(v), k), return(sumdigits(abs(k)))); listput(v, k)); c+=1)
for(n=1, 10^5, if(a(n)==2, print1(n, ", ")))
CROSSREFS
KEYWORD
nonn,base,fini,full
AUTHOR
Derek Orr, Jul 24 2017
STATUS
approved