[go: up one dir, main page]

login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A143933
a(n) is the smallest prime x such that x^2-n! is also prime.
3
2, 2, 3, 11, 19, 31, 79, 211, 607, 1931, 6337, 21961, 78919, 295291, 1143563, 4574149, 18859777, 80014843, 348776611, 1559776279, 7147792903, 33526120129, 160785623729, 787685471519, 3938427356629, 20082117944579, 104349745809137, 552166953567737
OFFSET
1,1
COMMENTS
Every prime > 3 in this sequence is bigger than the n-th prime, see comment to A121926. For the smallest number x such that x^2-n! is prime see A143931. For the smallest prime numbers of the form x^2-n! see A143932.
LINKS
MAPLE
f:= proc(n) local p, t;
t:= n!;
p:= floor(sqrt(t));
do
p:= nextprime(p);
if isprime(p^2-t) then return p fi
od
end proc:
map(f, [$1..28]); # Robert Israel, Feb 10 2019
MATHEMATICA
f[n_] := Block[{p = NextPrime[ Sqrt[ n!]]}, While[ !PrimeQ[p^2 - n!], p = NextPrime@ p]; p]; Array[f, 27] (* Robert G. Wilson v, Jan 08 2015 *)
PROG
(PARI) a(n)=my(N=n!, x=sqrtint(N)+1); while(!isprime(x^2-N), x=nextprime(x+1)); x \\ Charles R Greathouse IV, Dec 09 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Sep 05 2008
EXTENSIONS
Corrected by Charles R Greathouse IV, Dec 09 2014
STATUS
approved