OFFSET
1,1
EXAMPLE
Sum of square divisor of 15376, 15376+3844+961+16+4+1=20202, a palindrome with five digits.
MATHEMATICA
ssdpQ[n_]:=Module[{s=Total[Select[Divisors[n], IntegerQ[Sqrt[#]]&]]}, IntegerLength[ s]>1 &&s==IntegerReverse[s]]; Select[Range[ 900000], ssdpQ] (* The program uses the IntegerReverse function from Mathematica version 10 *) (* Harvey P. Dale, Jan 31 2016 *)
PROG
(PARI)
reverse(n)=concat(Vecrev(Str(n)))
palind(n)=(Str(n)==reverse(n)&&n>10)
{for(n=1, 10^6, k=sumdiv(n, d, d*issquare(d)); if(palind(k), print(n)))}
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Antonio Roldán, Dec 02 2013
STATUS
approved