OFFSET
1,1
COMMENTS
Inspired by Zak Seidov's post to the SeqFan list, cf. link: This yields A154673 as 252 a(n)^2. Indeed, if N/7 is a square, then N=7m^2 and this can't be the average of a twin prime pair unless m=6a.
LINKS
Zak Seidov, "A154676", Jan 15 2009
FORMULA
a(n) = sqrt(A154673(n)/252).
MATHEMATICA
Select[Range[600], And@@PrimeQ[252#^2+{1, -1}]&] (* Harvey P. Dale, Dec 13 2012 *)
PROG
(PARI) for(i=1, 999, isprime(252*i^2+1) & isprime(252*i^2-1) & print1(i", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Jan 15 2009
STATUS
approved