OFFSET
1,1
COMMENTS
The associated (p-1)^3/8+(p+1)^2/4 are in A163424.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
EXAMPLE
p=3 is in the sequence because (3-1)^3/8+(3+1)^2/4=1+4=5 is also prime.
p=5 is in the sequence because (5-1)^3/8+(5+1)^2/4=17 is also prime.
MATHEMATICA
f[n_]:=((p-1)/2)^3+((p+1)/2)^2; lst={}; Do[p=Prime[n]; If[PrimeQ[f[p]], AppendTo[lst, p]], {n, 7!}]; lst
Select[Prime[Range[1500]], PrimeQ[(# - 1)^3 / 8 + (# + 1)^2 / 4]&] (* Vincenzo Librandi, Apr 08 2013 *)
PROG
(Magma) [p: p in PrimesInInterval(3, 2000) | IsPrime((p-1)^3 div 8 + (p+1)^2 div 4)]; // Vincenzo Librandi, Apr 08 2013
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vladimir Joseph Stephan Orlovsky, Jul 27 2009
EXTENSIONS
Comment turned into examples by R. J. Mathar, Sep 02 2009
STATUS
approved