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: A083863 A057626 A063968 this_sequence A063831 A064023 A024350

Adjacent sequences: A064498 A064499 A064500 this_sequence 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