OFFSET
1,2
EXAMPLE
7776*(7*7*7*6) = 1600030008 = 252^3. Thus, 7776 is a member of this sequence.
PROG
(Python)
def DP(n):
..p = 1
..for i in str(n):
....p *= int(i)
..return p
def a(n):
..k = 0
..while k < n:
....if k**3 == n*DP(n):
......return n
....if k**3 > n*DP(n):
......return 0
....k += 1
n = 1
while n < 10**6:
..if a(n):
....print(n, end=', ')
..n += 1
# Simplified by Derek Orr, Mar 22 2015
(PARI) for(n=1, 10^7, d=digits(n); p=prod(i=1, #d, d[i]); if(p&&ispower(n*p, 3), print1(n, ", "))) \\ Derek Orr, Mar 22 2015
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Derek Orr, Sep 25 2013
EXTENSIONS
Corrected and extended by Derek Orr, Mar 22 2015
STATUS
approved