%I#31 2021年5月3日16:38:32

%S 1,4,11,30,66115200302441619829108513951771220026663228，

%电话：38434564535161857143815893491052611934133751489616652，

%电话：1838120370224112462926963294063210134848403776640920441644758751200

%N x^3+y^3+z^3的非负解的个数。

%H David A.Corneth，n的表，A（n）表示n=0..1500</a>

%H David A.Corneth，巴黎程序</a>

%F a（n）=[x^（n^3）]（1/（1-x））*（和{k>=0}x ^（k^3））^3.-_伊利亚·古特科夫斯基，2018年4月20日

%e对于n=1，四个解是{0,0,0}，{0,01}，}0,1,0}和{1,0,0{，因此a（1）=4。

%o（Python）

%o对于范围（99）内的a：

%o n=a*a*a

%o k=0

%o表示范围（99）内的x：

%o s=x*x*x

%o如果s>n：中断

%o表示范围（99）内的y：

%o sy=s+y*y*y

%o如果sy>n：中断

%o对于范围（99）中的z：

%o sz=sy+z*z*z

%o如果sz>n：中断

%o k+=1

%o打印（str（k），end='，'）

%o（PARI）a（n）=n++；p=Pol（（1/（1-x））*和（k=0，n，x^（k^3））^3+O（x^（n^3）））；波尔科夫（p，（n-1）^3）；\\_米歇尔·马库斯，2018年4月21日

%o（PARI）见PARI链接\\_David A.Corneth_，2018年5月22日

%Y参见A224214。

%K nonn公司

%O 0,2

%A _Alex Ratushnyak，2013年4月1日