OFFSET
1,1
COMMENTS
It is easily shown that all terms end in 7 or 9.
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..300
FORMULA
a(n) = 6*A166119(n) - 1. - Hugo Pfoertner, Sep 18 2024
EXAMPLE
59 is a term because it is the smaller member of the twin prime pair (59,61) and 59 + 61 = 120 = 15*(15+1)/2.
MATHEMATICA
Transpose[Select[Partition[Prime[Range[45000]], 2, 1], #[[2]]-#[[1]]==2 && IntegerQ[(Sqrt[1+8Total[#]]-1)/2]&]][[1]] (* Harvey P. Dale, Dec 15 2011 *)
PROG
(PARI) lista(n)={my(L=List(), k=0); while(#L<n, k++; my(t=k*(k+1)/4); if(!frac(t)&&isprime(t-1)&&isprime(t+1), listput(L, t-1))); Vec(L)} \\ Andrew Howroyd, Sep 18 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Aug 06 2003
EXTENSIONS
Offset changed by Andrew Howroyd, Sep 18 2024
STATUS
approved
