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
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)
Comments