%I#43 2023年10月27日08:19:44

%S 0,1,2,3,8,9,10,16,17,24,27,28,29,35,36,43,54,55,62,64,65,66,72,73,80，

%电话81,91,92,99118125126128129133134136141152153155160，

%U 179189190192197216217218224225232243244250251253

%N个数字，最多是3个正立方的总和。

%C或：3个（不一定是不同的）非负立方体之和的数字_R.J.Mathar，2015年9月9日

%C Deshouillers、Hennecart和Landreau推测，该序列的密度为0.0999425…=lim_K Sum_{K=1..K}-exp（C*rho（K，K）/K^2）/K，其中C=-gamma（4/3）^3/6=-0.1186788…，K在A003418（或等价地，A051451）中取递增值，rho（k0，K）是三元组数1<=k1，k2，k3<=K，这样k0=k1^3+k2^3+k3^3 mod K.-Charles R Greathouse IV_，2016年9月16日

%H Alois P.Heinz，n的表格，n=1..100000的a（n）（来自T.D.Noe的前1000个术语）

%H Jean-Marc Deshouillers、François Hennecart和Bernard Landreau，<a href=“http://link.springer.com/chapter/10.1007/11792086_11“>关于三个立方体总和的密度，ANTS-VI（2006），第141-155页。

%H<a href=“/index/Su#ssq”>与多维数据集和相关的序列的索引项</a>

%p是A004825:=进程（n）

%p局部x，y，zc；

%从0到do的x的p

%p如果3*x^3>n，则

%p返回false；

%p end if；

%x do中y的p

%p如果x^3+2*y^3>n那么

%p断裂；

%p其他

%pzc：=n-x^3-y^3；

%p如果zc>=y^3且isA000578（zc），则

%p返回true；

%p end if；

%p end if；

%p端do：

%p端do：

%p端程序：

%p A004825:=程序（n）

%p选项记住；

%p局部a；

%如果n=1，则为p

%p 0；

%p其他

%p表示来自procname（n-1）+1 do的a

%p如果是A004825（a），则

%p返回a；

%p end if；

%p端do：

%p end if；

%p端程序：

%p序列（A004825（n），n=1..100）；#_R.J.Mathar，2015年9月9日

%p#第二个Maple程序：

%pb:=proc（n，i，t）选项记忆；n=0或i>0且t>0

%p和（b（n，i-1，t）或i^3<=n和b（n-i^3，i，t-1））

%p端：

%p a：=proc（n）选项记住；局部k；

%p表示从1+`if`（n=1，-1，a（n-1））开始的k

%p而不是b（k，iroot（k，3），3）做od；k个

%p端：

%p序列（a（n），n=1..100）；#_Alois P.Heinz，2016年9月16日

%t q=7；imax=q^3；选择[Union[Flatten[Table[x^3+y^3+z^3，{x，0，q}，{y，x，q}，{z，y，q}]]，#<=imax&]（*_Vladimir Joseph Stephan Orlovsky_，2011年4月20日*）

%o（PARI）列表（lim）=我的（v=列表（），k，t）；对于（x=0，sqrtnint（lim=1,3），对于（y=0，min（sqrtinint（lim-x^3,3），x），k=x^3+y^3；对于（z=0，min（sqrtnint（lim-k，3），y），列表输入（v，k+z^3））；集（v）\\_Charles R Greathouse IV_，2015年9月14日

%Y A003072是一个子序列。

%Y参考A004999。

%A336820的Y列k=3。

%K nonn公司

%氧1,3

%A _N.J.A.斯隆_