%I#22 2021年5月9日11:20:08
%S 16037847649513661464148316217011799223233023492548,
%电话:2567266535363555365338714948504650655655264528353816252,
%电话:6271636965876894699270117118713772107229723573277453819882178315
%N个数,是不同素数的三个立方体之和。
%C这是A024975的子序列。这个序列的奇数项是A138853,偶数项是8+{偶数项A120398}。因此,该序列中的素数A137365与A138853中的素数相同。
%H Chai Wah Wu,n表,n(n)表示n=1..10000(R.J.Mathar的术语1..621)
%H<a href=“/index/Su#ssq”>与多维数据集总和相关的序列索引。
%F A138854={p(i)^3+p(j)^3+p(k)^3;i>j>k>0}=A138853并集{p(ii)^3+p(j
%p为A030078:=进程(n)
%p局部f;
%如果n<8,则为p
%p假;
%p其他
%p f:=系数(n)[2];
%p如果nops(f)=1且op(2,op(1,f))=3,则
%p为真;
%p其他
%p假;
%p end if;
%p end if;
%p结束过程:
%p是A138854:=进程(n)
%p局部i,j,p,q,r,rcub;
%从1do到i的p
%p p:=i素数(i);
%p如果p^3+(p+1)^3+
%p返回false;
%p end if;
%i+1 do中j的p
%pq:=i素数(j);
%p循环:=n-q^3-p^3;
%p如果rcub<=q^3,则
%p断裂;
%p fi;
%p如果是A030078(rcub),则
%p返回true;
%p end if;
%p端do:
%p端do:
%p结束过程:
%从5do到n的p
%p如果是A138854(n),则
%p打印(n);
%p end if;
%截止日期:#R.J.Mathar_,2014年6月9日
%t f[upto_]:=模块[{maxp=PrimePi[Floor[Power[upto,(3)^-1]]},选择[Union[Total/@(Subsets[Prime[Range[maxp]],{3}]^3)],#<=upto&]];f[9000](*哈维·P·戴尔,2011年3月21日*)
%o(PARI)是A138854(n)={如果(n%2,是A13885(n),是A120398(n-8))}
%o表示(n=1,10^4,isA138854(n)&打印1(n“,”))
%Y参考A024975(a^3+b^3+c^3,a>b>c>0),A138853(奇数项),A120398,A137365(A138853/A138854中的素数)。
%K nonn公司
%O 1,1号机组
%A.M.F.Hasler_,2008年4月13日
|