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