OFFSET
1,1
COMMENTS
Inspired by Zak Seidov's post to the SeqFan list, cf. link: This yields A154674 as 9900 a(n)^2. Indeed, if N/11 is a square, then N=11 m^2 and this can't be the average of a twin prime pair unless m=30a (considering N+1 mod 2,3,5 and N-1 mod 5).
LINKS
Zak Seidov, "A154676", Jan 15 2009
FORMULA
a(n) = sqrt(A154674(n)/9900).
MATHEMATICA
Select[Range[1000], AllTrue[9900#^2+{1, -1}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Jan 19 2019 *)
PROG
(PARI) for(i=1, 999, isprime(9900*i^2+1) & isprime(9900*i^2-1) & print1(i", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Jan 15 2009
STATUS
approved