OFFSET
1,3
COMMENTS
Sieving can be limited to odd values of k, because 3^k - k + 1 is even when k is even. In fact, if k is even, 3^k - k is odd and the successor is even.
MATHEMATICA
ListA[k_] := Block[{seq = {}, n = 0, i = 0}, While[Length[seq] < k, {n = 3^i - i + 1, If[PrimeQ[n], AppendTo[seq, i]], i += 1}]; seq]
PROG
(Sage)
def list_a(k):
return [i for i in range(k) if (3**i) - i + 1 in Primes()]
(PARI) isok(k) = isprime(3^k - k + 1); \\ Jinyuan Wang, Aug 03 2019
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Giuseppe Bonaccorso, Aug 02 2019
STATUS
approved