OFFSET
1,1
COMMENTS
Terms of this sequence must be of the form 13(6m)^2, the values for m are listed in A154775. - M. F. Hasler, Jan 15 2009
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = 468 A154775(n)^2 - M. F. Hasler, Jan 15 2009
MATHEMATICA
lst={}; Do[If[PrimeQ[n-1]&&PrimeQ[n+1], s=(n*13)^(1/2); If[Floor[s]==s, AppendTo[lst, n]]], {n, 6, 11!, 6}]; lst
(* ... and/or ... *)
lst={}; Do[If[PrimeQ[n-1]&&PrimeQ[n+1], s=(n/13)^(1/2); If[Floor[s]==s, AppendTo[lst, n]]], {n, 6, 11!, 6}]; lst
Select[13*Range[10^3]^2, And @@ PrimeQ[# + {-1, 1}] &] (* Amiram Eldar, Dec 25 2019 *)
PROG
(PARI) for(i=1, 499, isprime(468*i^2+1) && isprime(468*i^2-1) && print1(468*i^2", ")) \\ M. F. Hasler, Jan 15 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Jan 13 2009
EXTENSIONS
More terms from M. F. Hasler, Jan 15 2009
STATUS
approved