palidiv.py

f = lambda m, *a: sum([10**i*(a[i] - m*a[-i-1]) for i in range(len(a))])
sp = lambda x: [int(y) for y in str(x)]
dom = lambda k: xrange(10**k, 5*10**(k))
X = dom(5)

for m in xrange(2,10):
  for x in X:
    if f(m, *sp(x)) == 0: print x