OFFSET
1,1
COMMENTS
If k is a term, then so is 1000*k.
Contains 27 * (10^(3*k+2)+1)^3 for k >= 0.
EXAMPLE
a(3) = 216648648216 = 6006^3 is a term because it is the concatenation of 6^3 = 216, 4^3 = 64, 2^3 = 8, 4^3 = 64, 2^3 = 9 and 6^3 = 216.
MAPLE
g:= proc(x) local j, y;
for j from 1 to ilog10(x) do
y:= x mod 10^j;
if y < 10^(j-1) then next fi;
if surd(y, 3)::integer and h((x-y)/10^j) then return true fi;
od;
false
end proc:
h:= proc(x) option remember; local j, y;
if surd(x, 3)::integer then return true fi;
for j from 1 to ilog10(x) do
y:= x mod 10^j;
if y < 10^(j-1) then next fi;
if surd(y, 3)::integer and procname((x-y)/10^j) then return true fi;
od;
false
end proc:
select(g, [seq(x^3, x = 1 .. 10^7)]);
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Robert Israel, Apr 16 2025
STATUS
approved
