OFFSET
1,1
COMMENTS
a(18) > 2*10^5. - Robert Price, Oct 18 2013
EXAMPLE
3^3 - 22 = 5 (prime), so 3 is in the sequence.
MATHEMATICA
Do[If[PrimeQ[3^n - 22], Print[n]], {n, 10000}]
PROG
(PARI) is(n)=isprime(3^n-22) \\ Charles R Greathouse IV, Feb 17 2017
(Python)
from sympy import isprime
def ok(n): return isprime(3**n - 22)
print([m for m in range(700) if ok(m)]) # Michael S. Branicky, Mar 04 2021
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Nicolas M. Perrault, Nov 10 2012
EXTENSIONS
a(17) from Robert Price, Oct 18 2013
STATUS
approved