OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
FORMULA
a(n) = A046655(n)^2.
EXAMPLE
60516 = 246^2. Sum of square divisors: 60516 + 15129 + 6724 + 1681 + 36 + 9 + 4 + 1 = 84100 = 290^2.
MATHEMATICA
f[p_, e_] := (p^(2*(1 + Floor[e/2])) - 1)/(p^2 - 1); A035316[1] = 1; A035316[n_] := Times @@ f @@@ FactorInteger[n];
Select[Range[200000]^2, IntegerQ[Sqrt[A035316[#]]]&] (* Amiram Eldar, Aug 12 2023 *)
PROG
(PARI) {for(n=1, 10^5, m=n*n; k=sumdiv(m, d, d*issquare(d)); if(issquare(k), print(m)))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Antonio Roldán, Nov 26 2013
STATUS
approved