%I #20 Oct 27 2023 21:39:37
%S 64,1728,4096,17576,21952,46656,110592,195112,287496,314432,405224,
%T 474552,1061208,1191016,1404928,1601613,1906624,2000376,2146689,
%U 2197000,3241792,3511808,4913000,5268024,6229504,6751269,6859000,7077888,11239424,20346417,21952000
%N Cubes which are arithmetic mean of two consecutive primes.
%H Alois P. Heinz, <a href="/A234240/b234240.txt">Table of n, a(n) for n = 1..1000</a> (terms n = 1..57 from K. D. Bajpai)
%F a(n) = A075191(n)^3.
%e 64 is in the sequence because cube 64 = 4^3 = (61+67)/2 is arithmetic mean of two consecutive primes.
%e 1728 is in the sequence because 1728 = 12^3 = (1723+1733)/2.
%p a:= proc(n) option remember; local k, kk, p, q;
%p for k from 1 +`if`(n=1, 1, iroot(a(n-1), 3))
%p do kk:= k^3; p, q:= prevprime(kk), nextprime(kk);
%p if (p+q)/2=kk then return kk fi
%p od
%p end:
%p seq(a(n), n=1..60); # _Alois P. Heinz_, Dec 21 2013
%t Select[Mean/@Partition[Prime[Range[1500000]],2,1],IntegerQ[Surd[#,3]]&] (* _Harvey P. Dale_, Oct 08 2014 *)
%t Select[Range[300]^3,#==Mean[{NextPrime[#],NextPrime[#,-1]}]&] (* _Harvey P. Dale_, Sep 02 2015 *)
%o (PARI) is(n)=nextprime(n)+precprime(n)==2*n && ispower(n,3)
%o for(n=8,1e4,if(is(n^3), print1(n^3", "))) \\ _Charles R Greathouse IV_, Aug 25 2014
%Y Cf. A000578, A062703, A069495, A075191.
%K nonn
%O 1,1
%A _K. D. Bajpai_, Dec 21 2013
|