%I #20 Apr 04 2023 07:46:27
%S 0,1,4,5,6,9,10,11,12,21,24,25,29,32,33,34,39,40,44,49,50,51,54,55,56,
%T 59,60,61,64,65,66,67,71,72,73,75,76,88,90,97,99,100,101,102,106,109,
%U 110,111,112,114,115,116,119,120,124,125,129,137,151,153,176
%N Digits of n appear in n^3.
%H Robert Israel, <a href="/A029776/b029776.txt">Table of n, a(n) for n = 1..10000</a>
%e 66 is in the list because 66^3 = 287496 and 6 appears in 287496.
%e 153 is in the list because 153^3 = 3581577 and 1, 5, 3 appear in 3581577.
%p filter:= n -> convert(convert(n,base,10),set) subset convert(convert(n^3,base,10),set):select(filter, [$0..200]); # _Robert Israel_, Mar 18 2020
%o (Magma) [n: n in [0..200] | Intseq(n) subset Intseq(n^3)]; // _Bruno Berselli_, Aug 01 2013
%o (PARI) isok(m) = my(d=Set(digits(m)), ddd=Set(digits(m^3))); setintersect(d, ddd) == d; \\ _Michel Marcus_, Mar 18 2020
%o (Python)
%o from itertools import count, islice
%o def A029776_gen(startvalue=0): # generator of terms >= startvalue
%o return filter(lambda n:set(str(n)) <= set(str(n**3)), count(max(startvalue,0)))
%o A029776_list = list(islice(A029776_gen(),20)) # _Chai Wah Wu_, Apr 03 2023
%Y Contains A119735.
%K nonn,base
%O 1,3
%A _Patrick De Geest_
%E Offset changed to 1 by _Robert Israel_, Mar 18 2020
|