OFFSET
1,1
COMMENTS
EXAMPLE
Prime 17 is in sequence because sigma(2*17-1) = sigma(33) = 48 = 3*(17-1).
MATHEMATICA
Select[Prime@ Range[10^5], DivisorSigma[1, 2 # - 1] == 3 (# - 1) &] (* Michael De Vlieger, Jan 03 2017 *)
PROG
(Magma) [p: p in PrimesUpTo(20000)| SumOfDivisors(2*p-1) eq 3*p-3]
(PARI) forprime(p=1, 10^7, if(sigma(2*p-1)==3*(p-1), print1(p, ", "))) \\ Derek Orr, Sep 25 2014
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Jaroslav Krizek, Sep 24 2014
EXTENSIONS
a(8) from Matthew Campbell, Jan 03 2017
STATUS
approved