OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1800 from Vincenzo Librandi)
FORMULA
A093616(a(n)) = n.
MATHEMATICA
A093616[n_] := For[k = 1, True, k++, If[DivisorSigma[0, k*n] == DivisorSigma[0, n^2], Return[k]]]; Select[Range[100], A093616[#] == # &] (* Jean-François Alcover, Aug 14 2014 *)
PROG
(PARI) is(n) = {my(k = 1, d = numdiv(n^2)); while(k < n && numdiv(k*n) != d, k++); k == n; } \\ Amiram Eldar, Apr 15 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Apr 06 2004
EXTENSIONS
Corrected by Franklin T. Adams-Watters, Dec 19 2006
STATUS
approved