|
|
A229627
|
|
a(n) is the smallest prime q such that 2*q^k - 1 is prime for k = 1, 2, ..., n.
|
|
2
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
The prime number 2 in the definition is used because 2 is the only prime p such that p*q^k - 1 can be prime for more than one prime q.
|
|
LINKS
|
|
|
MATHEMATICA
|
a[1]=2; a[n_]:=a[n]=(For[m=PrimePi[a[n-1]], Union[Table[PrimeQ[2 Prime[m]^k-1], {k, n}]]!={True}, m++]; Prime[m])]
|
|
PROG
|
(PARI) a(n)=forprime(m=2, , for(k=1, n, if(!ispseudoprime(2*m^k-1), next(2))); return(m)) \\ Charles R Greathouse IV, Oct 01 2013
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|