1,2

A093616(a(n)) = n; complement of A093617.

Vincenzo Librandi, Table of n, a(n) for n = 1..1800

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 *)

Cf. A000005, A000290, A048691.

Sequence in context: A090322 A076084 A151764 * A317257 A242031 A335376

Adjacent sequences: A093615 A093616 A093617 * A093619 A093620 A093621

nonn

Reinhard Zumkeller, Apr 06 2004

Corrected by Franklin T. Adams-Watters, Dec 19 2006

approved