OFFSET
1,1
COMMENTS
For any term a(n), all numbers a(n)*10^k, k >= 0, are also in the sequence. More interestingly, all numbers N == 471 (mod 1000) are in the sequence, since 471^3*3 == 333 (mod 1000). - M. F. Hasler, Jul 16 2024
REFERENCES
C. A. Pickover, Keys to Infinity. New York: Wiley, p. 7, 1995.
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
Giovanni Resta, super-d numbers, personal web site "Numbers Aplenty", 2013.
Eric Weisstein's World of Mathematics, Super-d Number.
EXAMPLE
n = 1752 so 3 * 1752^3 = 161'333'13024.
MATHEMATICA
Select[Range[5000], MemberQ[Partition[IntegerDigits[3#^3], 3, 1], {3, 3, 3}]&] (* Harvey P. Dale, Feb 01 2013 *)
PROG
(PARI) select( {is_A014569(n, d=3, m=10^d, r=m\9*d)=n=d*n^d; until(r>n\=10, n%m==r && return(1))}, [0..4999]) \\ Using the (optional) 2nd arg d=2..9 allows to compute the sequences A032743-A032749. - M. F. Hasler, Jul 16 2024
(Python) is_A014569=lambda n, d=3: str(d)*d in str(d*n**d) # M. F. Hasler, Jul 16 2024
CROSSREFS
KEYWORD
nonn,base
AUTHOR
EXTENSIONS
Corrected and extended by Patrick De Geest, May 15 1998
Offset changed to 1 by M. F. Hasler, Jul 16 2024
STATUS
approved