%I#35 2015年2月8日08:49:07

%S 10,11,29,34,99100101106110114119120124125274275276，

%电话：279281290296299314315313193203203243253293303353340，

%电话：966970975976979986990909969100010011004105100610101020102110241025103410491051

%N数字N，这样可以通过插入内部（但不一定是连续的）数字从N中获得N^3。

%C这是A046851，因为立方体A000578是方形A000290。A086458的子集（但请注意，即104^3=1124864以与104相同的数字开始和结束，但缺少内部“0”）。如果我们需要插入的数字来填充相邻的位置，则会产生另一个序列，例如不包含106。

%e 34^3=39304（在“34”中插入“930”）。

%e 106^3=1191016（在“106”中插入“191”和“1”）。

%p A000030:=proc（n），如果n=0，则为0；else op（-1，转换（n，base，10））；结束条件：；结束过程：

%p A010879:=程序（n）n模块10；结束过程：

%p为A086458：=程序（n）A000030（n）=A000030；结束过程：

%p subsI：=从1到nops（L）的i的proc（c，L）do，如果op（i，L）=c，则返回i；结束条件：；结束do；返回-1；结束过程：

%p isSubS:=proc（Sub，Sup）如果nops（Sub）=1，则如果subsI（op（1，Sub），Sup）>0，则返回true；否则返回false；结束条件：；elif nops（Sub）=0，然后返回true；否则f：=subsI（op（1，Sub），Sup）；如果f<0，则返回false；else进程名（subsop（1=NULL，Sub），[op（f+1..nops（Sup），Sup）]）；结束条件：；结束条件：；结束过程：

%p是A179856：=proc（n），如果是A086458（n）则dgsn：=convert（n，base，10）；dgsn：=op（2..nops（dgsn）-1，dgsn；dgsn3:=转换（n^3，基数，10）；dgsn3:=op（2..nops（dgsn3）-1，dgsn2）；是SubS（[dgsn]，[dgsn3]）；否则为假；结束条件：；结束过程：

%p表示n从10到1400，如果是A179856（n），则打印f（“%d，”，n）；结束条件：；结束日期：#R.J.Mathar，2011年1月30日

%Y参考A000578、A046851、A086458。

%K nonn，基础，简单

%O 1,1号机组

%《乔纳森·沃斯邮报》，2011年1月28日