cp's OEIS Frontend

A240859 Cubes k^3 such that k^3 + (k+1)^3 is semiprime.

%I A240859 #12 Dec 15 2015 11:55:48
%S A240859 1,8,27,125,216,512,2744,3375,8000,9261,35937,68921,125000,157464,
%T A240859 328509,421875,474552,704969,729000,970299,1157625,1367631,1685159,
%U A240859 2248091,2628072,2803221,3581577,3723875,4741632,5177717,5451776,6751269,7301384,9129329,9938375
%N A240859 Cubes k^3 such that k^3 + (k+1)^3 is semiprime.
%H A240859 K. D. Bajpai, <a href="/A240859/b240859.txt">Table of n, a(n) for n = 1..3689</a>
%e A240859 a(2) = 8 = 2^3: 2^3 + 3^3 = 35 = 5*7 which is a semiprime.
%e A240859 a(4) = 125 = 5^3: 5^3 + 6^3 = 341 = 11*31 which is a semiprime.
%p A240859 with(numtheory):KD:= proc() local a,b; a:=n^3+(n+1)^3;b:=bigomega(a); if b=2 then RETURN (n^3); fi; end: seq(KD(), n=1..500);
%t A240859 Transpose[Select[Partition[Range[250]^3,2,1],PrimeOmega[Total[#]] == 2&]] [[1]] (* _Harvey P. Dale_, Dec 15 2015 *)
%o A240859 (PARI) forprime(p=3,1e4,if(isprime((p^2+3)/4),print1((p^3 - 3*p^2 + 3*p - 1)/8", "))) \\ _Charles R Greathouse IV_, Aug 25 2014
%Y A240859 Cf. A001358, A005898.
%K A240859 nonn,easy
%O A240859 1,2
%A A240859 _K. D. Bajpai_, Apr 13 2014