OFFSET
1,2
EXAMPLE
17 is in the sequence because 4^sigma(17) - 17 = 4^18 - 17 = 68719476719 is prime.
MATHEMATICA
a[n_] := Select[Range@ n, PrimeQ[4^DivisorSigma[1, #] - #] &]; a[20000]
DeleteCases[ParallelTable[If[PrimeQ[4^DivisorSigma[1, k]-k], k, n], {k, 1, 10^4}], n]
PROG
(Magma) [n: n in[1..10000] | IsPrime((4^SumOfDivisors(n)) - n)];
CROSSREFS
KEYWORD
nonn,more
AUTHOR
J.W.L. (Jan) Eerland, Nov 11 2024
EXTENSIONS
a(9) from Michael S. Branicky, Nov 11 2024
STATUS
approved