OFFSET
1,2
COMMENTS
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
EXAMPLE
929 - (9*2*9) = 767 (another palindrome). So, 929 is a member of this sequence.
PROG
(Python)
def DP(n):
..p = 1
..for i in str(n):
....p *= int(i)
..return p
def pal(n):
..r = ''
..for i in str(n):
....r = i + r
..return r == str(n)
{print(n, end=', ') for n in range(1, 10**6) if DP(n) and pal(n) and pal(n-DP(n))}
## Simplified by Derek Orr, Apr 05 2015
(PARI) pal(n)=d=digits(n); Vecrev(d)==d
for(n=1, 10^7, d=digits(n); p=prod(i=1, #d, d[i]); if(p&&pal(n)&&pal(n-p), print1(n, ", "))) \\ Derek Orr, Apr 05 2015
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
Derek Orr, Sep 30 2013
EXTENSIONS
More terms from Derek Orr, Apr 05 2015
STATUS
approved