login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A328445
a(n) is the smallest prime p such that n = Omega(p^n - 2) = Omega(p^n + 2) where Omega = A001222.
0
5, 11, 127, 401, 1487, 1153, 6199, 10301, 22193, 72277, 1301423
OFFSET
1,1
EXAMPLE
5 is a term of a(1) because 1 = Omega(5) = Omega(7),
11 is a term of a(2) because 2 = Omega(119) = Omega(123).
MATHEMATICA
a[n_] := Module[{p = 2}, While[PrimeOmega[p^n - 2] != n || PrimeOmega[p^n + 2] != n, p = NextPrime[p]]; p]; Array[a, 10] (* Amiram Eldar, Oct 15 2019 *)
PROG
(PARI) a(n) = {my(p=3); while (! ((bigomega(p^n-2) == n) && (bigomega(p^n+2) == n)), p = nextprime(p+1)); p; } \\ Michel Marcus, Oct 17 2019
CROSSREFS
Cf. A001222 (bigomega), A328359.
Sequence in context: A248253 A262309 A266421 * A094108 A222674 A197538
KEYWORD
nonn,more
AUTHOR
EXTENSIONS
a(11) from Daniel Suteu and Giovanni Resta, Nov 07 2019
STATUS
approved