projecteuler/problem34.py

#!/usr/bin/env python
from itertools import combinations
from math import factorial

def digits(n):
    return map(int, str(n))

def facsum(n):
    return sum(facs[d] for d in digits(n))

facs = map(factorial, range(10))
s = 0

for n in xrange(10, 99999):
    if n == facsum(n):
        s += n

print s