projecteuler/55.py

#!/usr/bin/env python

def rev(n):
    r = 0

    while n > 0:
        r = r * 10 + n % 10
        n /= 10

    return r

def is_lychrel(n):
    niter = 0

    while niter < 50:
        n += rev(n)

        if rev(n) == n:
            return False

        niter += 1

    return True

print sum(int(is_lychrel(n)) for n in xrange(1, 10000))