OFFSET
0,1
COMMENTS
a(7) > 46472.
a(7) > 2*10^7. - Michael S. Branicky, Nov 09 2024
MATHEMATICA
Table[m = 3;
k = 0; Monitor[
Parallelize[
While[True,
If[And[PrimeQ[m^((m + 1)*m^n) - k],
PrimeQ[((m^((m + 1)*m^n) - k) - 1)/2]], Break[]]; k++]; k],
k], {n, 0, 5}]
PROG
(PARI) a(n) = {my(k=0); while (!(isprime(p=3^(4*3^n) - k) && isprime((p-1)/2)), k++); k; }
(Python)
from sympy import isprime, prevprime
def A(n):
m = 3**(4*3**n)
p = prevprime(m)
while not isprime((p-1)//2):
p = prevprime(p)
return m-p #
CROSSREFS
KEYWORD
nonn,more
AUTHOR
J.W.L. (Jan) Eerland, Oct 10 2024
STATUS
approved