OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from T. D. Noe)
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems: Inversion of Multiplicative Functions (invphi.gp).
FORMULA
a(n) = A228061(n)^2. - Amiram Eldar, Aug 12 2024
MATHEMATICA
nn = 14400; t = Select[Union[DivisorSigma[1, Range[nn]]], IntegerQ[Sqrt[#]] &]; t = Select[t, # <= nn &] (* T. D. Noe, Sep 04 2013 *)
PROG
(PARI) lista(kmax) = for(k = 1, kmax, if(invsigmaNum(k^2) > 0, print1(k^2, ", "))); \\ Amiram Eldar, Aug 12 2024, using Max Alekseyev's invphi.gp
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved