Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #14 Sep 18 2024 14:15:31
%S 17,59,137,149,1277,1619,2729,4127,6359,11717,13397,17489,20807,23027,
%T 28307,29669,51869,61379,69827,73847,78539,92567,97499,115769,138569,
%U 166667,176609,193379,200927,258317,304979,430007,500909,518759,536189
%N Smaller member of a twin prime pair with a triangular sum.
%C It is easily shown that all terms end in 7 or 9.
%H Harvey P. Dale, <a href="/A086816/b086816.txt">Table of n, a(n) for n = 1..300</a>
%F a(n) = 6*A166119(n) - 1. - _Hugo Pfoertner_, Sep 18 2024
%e 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.
%t Transpose[Select[Partition[Prime[Range[45000]],2,1],#[[2]]-#[[1]]==2 && IntegerQ[(Sqrt[1+8Total[#]]-1)/2]&]][[1]] (* _Harvey P. Dale_, Dec 15 2011 *)
%o (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
%Y Cf. A000217, A001359, A166119.
%K nonn
%O 1,1
%A _Jason Earls_, Aug 06 2003
%E Offset changed by _Andrew Howroyd_, Sep 18 2024