OFFSET
1,2
COMMENTS
EXAMPLE
MATHEMATICA
f[n_] := Length@ Select[ 2n - Prime@ Range@ PrimePi@ n, PrimeQ]; fQ[n_] := Mod[ 2n*f[n], DivisorSigma[1, 2n]] == 0; k = 1; lst = {}; While[k < 1250001, If[ fQ@ k, AppendTo[lst, k]; Print@ k]; k++]; lst (* Robert G. Wilson v, Aug 07 2014 *)
PROG
(PARI)
for(n=1, 10^7, my(s); forprime(p=2, n, s+=isprime(2*n-p)); d=divisors(2*n); if(2*n*s%(sum(i=1, #d, d[i]))==0, print1(n, ", "))) \\ Derek Orr, Jul 27 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Ivan N. Ianakiev, Jul 27 2014
EXTENSIONS
a(18)-a(24) from Derek Orr, Jul 27 2014
a(25)-a(48) from Robert G. Wilson v, Aug 07 2014
STATUS
approved