cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A154774 Numbers n such that 9900n^2 is the average of a twin prime pair.

This page as a plain text file.
%I A154774 #9 Jan 19 2019 13:39:26
%S A154774 10,14,15,25,60,74,76,87,127,129,130,140,171,196,207,224,259,263,302,
%T A154774 314,315,319,333,337,377,389,451,454,470,491,508,518,549,568,574,589,
%U A154774 592,624,629,636,690,696,729,748,753,769,770,771,781,791,802,833,844
%N A154774 Numbers n such that 9900n^2 is the average of a twin prime pair.
%C A154774 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).
%H A154774 Zak Seidov, <a href="http://zak08.livejournal.com/4070.html">"A154676"</a>, Jan 15 2009
%F A154774 a(n) = sqrt(A154674(n)/9900).
%t A154774 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 *)
%o A154774 (PARI) for(i=1,999, isprime(9900*i^2+1) & isprime(9900*i^2-1) & print1(i","))
%Y A154774 Cf. A037073, A154331, A154772-A154773.
%K A154774 nonn
%O A154774 1,1
%A A154774 _M. F. Hasler_, Jan 15 2009