OFFSET
1,2
COMMENTS
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..91
EXAMPLE
42624 is divisible by 4*2*6*2*4 and by 4+2+6+2+4.
PROG
(Python)
from operator import mul
from functools import reduce
from gmpy2 import t_mod, mpz
A117228 = sorted([mpz(n) for n in (str(x)+str(x)[::-1] for x in range(1, 10**8))
if not (n.count('0') or t_mod(mpz(n), sum((mpz(d) for d in n)))
or t_mod(mpz(n), reduce(mul, (mpz(d) for d in n))))]+
[mpz(n) for n in (str(x)+str(x)[-2::-1] for x in range(10**8))
if not (n.count('0') or t_mod(mpz(n), sum((mpz(d) for d in n)))
or t_mod(mpz(n), reduce(mul, (mpz(d) for d in n))))])
# Chai Wah Wu, Aug 25 2014
(PARI)
rev(n)=r=""; d=digits(n); for(i=1, #d, r=concat(Str(d[i]), r)); eval(r)
for(n=1, 10^7, d=digits(n); if(rev(n)==n, p=prod(i=1, #d, d[i]); if(p&&n%p==0&&n%sumdigits(n)==0, print1(n, ", ")))) \\ Derek Orr, Aug 25 2014
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Giovanni Resta, Apr 22 2006
EXTENSIONS
More terms from Chai Wah Wu, Aug 22 2014
STATUS
approved