OFFSET
1,2
EXAMPLE
a(18) = 1^2+2^2+3^2+6^2+9^2+1^2+8^2 = 196 because divisors of 18 are 1,2,3,6,9,18.
a(31) = 1^2+3^2+1^2 = 11 because divisors of 31 are 1,31.
MATHEMATICA
Table[Total@ Flatten@ Map[#^2 &, IntegerDigits@ Divisors@ n], {n, 61}] (* Michael De Vlieger, Sep 24 2016 *)
PROG
(PARI) ssd(n) = my(d=digits(n)); sum(k=1, #d, d[k]^2);
a(n) = sumdiv(n, d, ssd(d)); \\ Michel Marcus, Sep 22 2016
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Bhushan Bade, Sep 22 2016
EXTENSIONS
More terms from Michel Marcus, Sep 22 2016
STATUS
approved