%I#26 2023年1月9日07:41:25
%S 0,0,00,0,1,0,0,0,0',0,0,
%T 0,0,1,0,0,0,0,00,0'0,0',0,0,
%U 0,0,00,0,1,0,0,0,0',0,0
%N将N划分为3个不同的正立方体的分区数。
%换句话说,方程n=x^3+y^3+z^3的解的个数,其中x>y>z>0_Antti Karttunen,2017年8月29日
%H Antti Karttunen,n的表,n=0..17073的a(n)</a>
%H<a href=“/index/Su#ssq”>与多维数据集和相关的序列的索引项</a>
%F a(n)=A025465(n)-A025468(n).-_Antti Karttunen,2017年8月29日
%e摘自2017年8月29日的《安蒂·卡图宁》:(开始)
%e对于n=36,有一个解:36=27+8+1,因此a(36)=1。
%e对于n=1009,有两个解:1009=10^3+2^3+1^3=9^3+6^3+4^3,因此a(1009)=2。这也是序列达到大于1的值的第一个点。
%e(结束)
%p A025469:=程序(n)
%p局部a,x,y,zcu;
%p a:=0;
%从1do到x的p
%p如果3*x^3>n那么
%p返回a;
%p end if;
%x+1 do中y的p
%p如果x^3+2*y^3>n那么
%p断裂;
%p end if;
%pzcu:=n-x^3-y^3;
%p如果zcu>y^3且为A000578(zcu),则
%pa:=a+1;
%p end if;
%p端do:
%p端do:
%p端程序:
%p序列(A025469(n),n=1..80);#_R.J.Mathar,2018年6月15日
%t表[Count[Integer Partitions[n,{3}],_?(和[UnsameQ@@#,AllTrue[#,IntegerQ[#^(1/3)]&]]&)],{n,105}](*_Michael De Vlieger_,2017年8月29日*)
%o(PARI)A025469(n)={my(s=0_Antti Karttunen,2017年8月29日
%Y参考A025465(不一定不同)、A025468、A025419(贪婪逆)。
%Y参见A024975(非零项位置)、A024974(项位置>1)、A025399-A025402。
%K nonn公司
%O 01010号
%A·热心的W·威尔逊_
|