OFFSET
1,1
COMMENTS
EXAMPLE
52 is in the sequence since the divisors of 52 are {1, 2, 4, 13, 26 & 52} so d(52) = 98 and 2*52 - 98 = 6, a perfect number.
MATHEMATICA
fQ[n_] := PerfectNumberQ[2n - DivisorSigma[1, n]]; Select[ Range@ 2500000, fQ]
PROG
(PARI) for(n=1, 25*10^5, d=2*n-sigma(n); d>0&&sigma(d)==2*d&&print1(n ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v and Waldemar Puszkarz, Apr 01 2018
STATUS
approved