OFFSET
1,2
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..6228
EXAMPLE
9801=3^4*11^2 and 1089=3^2*11^2.
PROG
(PARI) is(n)=ispowerful(n) && ispowerful(subst(Polrev(digits(n)), 'x, 10)) \\ Charles R Greathouse IV, Sep 16 2014
(PARI) has(n)=ispowerful(subst(Polrev(digits(n)), 'x, 10))
list(lim)=my(v=List(), t, t2); for(m=1, lim^(1/3), t=m^3; for(n=1, sqrtint(lim\t), if(has(t2=t*n^2), listput(v, t2)))); Set(v) \\ Charles R Greathouse IV, Sep 16 2014
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Giovanni Resta, Jan 28 2006
STATUS
approved