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.

A240169 Numbers n such that (6n)^3 is the sum of a twin prime pair.

This page as a plain text file.
%I A240169 #18 Aug 03 2014 18:16:02
%S A240169 1,29,65,81,99,136,165,174,176,191,200,266,295,301,319,346,351,370,
%T A240169 400,411,431,434,436,456,491,494,526,541,599,651,676,714,746,790,924,
%U A240169 956,991,1011,1131,1161,1194,1259,1274,1280,1304,1374,1550,1641,1644,1649,1714,1715,1739,1804,1811,1814,1830,1879,1941,2000
%N A240169 Numbers n such that (6n)^3 is the sum of a twin prime pair.
%C A240169 No terms end with 2, 3, 7, 8. Minimal differences are 1; e,g., a(52) - a(51) = 1715 - 1714. There are no three consecutive terms with common difference 1.
%C A240169 Distribution of last digits for first 61000 terms: W(0..9) = (10190, 10162, 0, 0, 10178, 10222, 10027, 0, 0, 10221).
%C A240169 For "existing" digits distribution is rather uniform.
%H A240169 Jens Kruse Andersen, <a href="/A240169/b240169.txt">Table of n, a(n) for n = 1..10000</a>
%F A240169 a(n) = (1/6)*A245591(n+1)^(1/3).
%e A240169 m = 1: (6m)^3 = 216 = 107 + 109, m = 29: (6m)^3 = 5268024 = 2634011 + 2634013.
%p A240169 select(n -> isprime(108 * n^3 - 1) and isprime(108 * n^3 + 1), [$1..1000]); # _Robert Israel_, Aug 03 2014
%t A240169 Select[Range[1000], PrimeQ[216#^3/2 - 1] && PrimeQ[216#^3/2 + 1] &] (* _Alonso del Arte_, Aug 02 2014 *)
%o A240169 (PARI) N=2*10^3; for(k=1,N,p=216*k^3; if(isprime(p/2-1)&&isprime(1+p/2), print1(k, ", ")))
%Y A240169 Cf. A245591.
%K A240169 nonn
%O A240169 1,2
%A A240169 _Zak Seidov_, Aug 02 2014