OFFSET
0,2
COMMENTS
All numbers in this sequence, except for a(0), are even.
LINKS
Marius A. Burtea, Table of n, a(n) for n = 0..1000
EXAMPLE
1*3^1+1 = 4 is not prime. 2*3^1+1 = 7 is prime. Thus, a(1) = 2.
1*3^3+1 = 28 is not prime. 2*3^3+1 = 57 is not prime. 3*3^3+1 = 82 is not prime. 4*3^3+1 = 109 is prime. Thus, a(3) = 4.
PROG
(Python)
import sympy
from sympy import isprime
def Pow3(n):
..for k in range(10**4):
....if isprime(k*(3**n)+1):
......return n
n = 1
while n < 100:
..print(Pow3(n))
..n += 1
(PARI)
for(n=0, 100, k=0; while(!isprime(k*3^n+1), k++); print1(k, ", ")) \\ Colin Barker, Mar 24 2014
(Magma) sol:=[]; m:=1; for n in [0..73] do k:=0; while not IsPrime(k*3^n+1) do k:=k+1; end while; sol[m]:=k; m:=m+1; end for; sol; // Marius A. Burtea, Jun 05 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Mar 23 2014
STATUS
approved