A344263 Numbers m such that 3^(2m+1) - 3^m + 1 is prime.

0, 3, 4, 11, 35, 56, 88, 104, 476, 1367, 1707, 2472, 22232, 25260

Offset: 1

Author

Reza K Ghazi, May 13 2021

Keywords

• nonn
• more

Comments

a(14) > 30000.

a(14) > 50000. - Michael S. Branicky, Aug 12 2024

Crossrefs

Cf. A344170.

Programs

Maple

for m from 0 to 3000 do if isprime(3^(2*m + 1) - 3^m + 1) then print(m); end if; end do;

Mathematica

Do[If[PrimeQ[3^(2 m + 1) - 3^m + 1], Print[m]], {m, 0, 3000}]
Select[Range[0,2500],PrimeQ[3^(2#+1)-3^#+1]&] (* Harvey P. Dale, Mar 01 2023 *)

PARI

for(m=0, 3e3, if(isprime(3^(2*m+1)-3^m+1), print1(m", ")))

SageMath

for m in range(3000):
    if is_prime(3^(2*m + 1) - 3^m + 1):
        print(m)