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.
MATHEMATICA
bpQ[n_]:=DigitCount[n, 10, 0]==0&&AllTrue[{n, n-Times@@IntegerDigits[n]}, PalindromeQ]; Select[Range[16*10^5], bpQ] (* Harvey P. Dale, Nov 11 2024 *)
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
