%I#15 2020年12月11日22:57:11

%S 1,1,1,1,2,1,1,3,1,6,5,9,10,25,32,51,9714424446376710622005，

%电话：4177571612101215263530664629114871205337372317718410，

%电话：122632023611124308192730138414615750263820954763120091388286171931627302867194578843590112232587

%N立方体和为立方体的不同正整数的非空集的数目，一个集的最大整数为N。

%e a（13）=10个集合：{13}、{2、3、8、13}，{4、8、11、12、13}、{1、2、6、7、11、13}，{2、5、7、8、12、13%，{3、4、8，10、11、12,13}，3,5,7,8,9,10,11,12,13}。

%o（Python）

%o从functools导入lru_cache

%o def perf_cube（n）：返回圆（n**（1/3））**3==n

%o@lru_cache（最大大小=无）

%o定义b（n，soc，c）：

%o如果n==0：

%o如果perf_cube（soc）：返回1

%o返回0

%o返回b（n-1，soc，c）+b（n-1，soc+n*n*n，c+1）

%o a=λn：b（n-1，n*n*n，1）

%o打印（[a（n）代表范围（1，21）中的n）]#_Michael S.Branicky_，2020年12月10日

%Y参考A000578、A126111、A339612。

%K nonn公司

%O 1,5型

%A _Ilya Gutkovskiy_，2020年12月10日

%E a（24）-a（41）摘自_Michael S.Branicky_，2020年12月10日

%E a（42）-a（45）摘自_Alois P.Heinz，2020年12月11日