%I #33 Mar 19 2023 23:20:10
%S 4812191,6353029,8039333,8821867,19876711,60742631,85017061,108879847,
%T 127042367,138853049,170367959,238190951,259108427,414949357,
%U 485941193,512095739,529218559,582868471,623331491,648485381,771656657,1001132351,1098706507,1172752457
%N The smallest of four consecutive primes whose sum is a perfect cube.
%H K. D. Bajpai and Chai Wah Wu, <a href="/A248587/b248587.txt">Table of n, a(n) for n = 1..10000</a> n = 1..42 from K. D. Bajpai.
%e a(2) = 6353029 is prime. Next three primes are 6353033, 6353051 and 6353071. Their sum = 6353029 + 6353033 + 6353051 + 6353071 = 25412184 = 294^3.
%e a(3) = 8039333 is prime. Next three primes are 8039359, 8039363 and 8039377. Their sum = 8039333 + 8039359 + 8039363 + 8039377 = 32157432 = 318^3.
%t t = {}; p = 2; q = 3; r = 5; Do[v = NextPrime[r]; If[IntegerQ[(p + q + r + v)^(1/3)], AppendTo[t, p]; Print[p]]; p = q; q = r; r = v, {5*10^8}]; t
%t Select[Partition[Prime[Range[6*10^7]], 4, 1],IntegerQ[Surd[Total[#], 3]] &] [[All, 1]] (* _Harvey P. Dale_, Oct 07 2016 *)
%o (PARI) lista(nn) = {vp = primes(nn); for (i=1, #vp - 3, if (ispower(vp[i]+vp[i+1]+vp[i+2]+vp[i+3], 3), print1(vp[i], ", ")););} \\ _Michel Marcus_, Oct 24 2014
%o (Python)
%o from sympy import nextprime, prevprime
%o A248587_list = []
%o for i in range(3,10**6):
%o n = i**3
%o p3 = prevprime(n//4)
%o p2, p4 = prevprime(p3), nextprime(p3)
%o p1 = prevprime(p2)
%o q = p1+p2+p3+p4
%o while q <= n:
%o if q == n:
%o A248587_list.append(p1)
%o p1, p2, p3, p4 = p2, p3, p4, nextprime(p4)
%o q = p1+p2+p3+p4 # _Chai Wah Wu_, Dec 31 2015
%Y Cf. A000040 (primes), A000578 (cubes).
%Y Cf. A061308 (two consecutive primes), A210205 (three consecutive primes).
%K nonn
%O 1,1
%A _K. D. Bajpai_, Oct 09 2014
|