1,2
(PARI) {usigma(n, s=1, fac, i)= fac=factor(n); for(i=1, matsize(fac)[1], s=s*(1+fac[i, 1]^fac[i, 2]); ); return(s); } for(n=1, 10^6, if(issquare(usigma(n^3)), print1(n, " ")))
Cf. A034448.
Sequence in context: A201310 A063968 A172136 * A063831 A064023 A179959
Adjacent sequences: A064498 A064499 A064500 * A064502 A064503 A064504
nonn
Jason Earls (zevi_35711(AT)yahoo.com), Oct 05 2001
More terms from Naohiro Nomoto (n_nomoto(AT)yabumi.com), Oct 22 2001
