OFFSET
1,1
COMMENTS
If it exists, a(17) > 120000. - Michael S. Branicky, Aug 19 2024
EXAMPLE
5 is in the sequence because 2^sigma(5)-5 = 2^6-5 = 59 is prime.
MATHEMATICA
a[n_] := Select[Range@ n, PrimeQ[2^DivisorSigma[1, #] - #] &]; a[20000]
DeleteCases[ParallelTable[If[PrimeQ[2^DivisorSigma[1, k]-k], k, n], {k, 1, 10^4}], n]
PROG
(Magma) [n: n in[1..10000] | IsPrime((2^SumOfDivisors(n)) - n)];
CROSSREFS
KEYWORD
nonn,more
AUTHOR
J.W.L. (Jan) Eerland, Jan 02 2024
EXTENSIONS
a(16) from J.W.L. (Jan) Eerland, Jan 25 2024
STATUS
approved