OFFSET
1,1
COMMENTS
Trivially if p is prime then (p+1)*sigma(p) = (p+1)^2.
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..1000
EXAMPLE
1156*sigma(1155) = 1632^2, hence 1155 is in the sequence.
MATHEMATICA
Select[Range[505000], ! PrimeQ[#] && IntegerQ[Sqrt[(# + 1)*DivisorSigma[1, #]]] &] (* Jayanta Basu, Jul 02 2013 *)
PROG
(PARI) isA073586(n)=(!isprime(n))&&issquare((n+1)*sigma(n)) \\ Michael B. Porter, Apr 05 2010
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Aug 28 2002
STATUS
approved