fatorial_100.py

from unittest import TestCase

import math

class FatorialDeN(TestCase):
    def test_soma_do_resultado_do_fatorial_de_10(self):
        self.assertEquals(27, soma_do_fatorial_de_n(n=10))

    def test_soma_do_resultado_do_fatorial_de_100(self):
        self.assertEquals(648, soma_do_fatorial_de_n(n=100))

def soma_do_fatorial_de_n(n):
    fatorial_n = math.factorial(n)
    resultado = 0

    for numero in  str(fatorial_n):
        resultado += int(numero)
    return resultado

Sugiro a seguinte função pra calcular o fatorial:

def fatorial(n):
    return reduce(lambda x, y: x*y, range(1,n+1))

Acho que ela vai propiciar uma performance bem maneira. Se combinar com o list comprehension que o @rennerocha sugeriu então, ficará ótimo!

Douglas,
A biblioteca padrão já possui a função math.factorial. Com certeza a performance dela já foi otimizada, então não vejo válido usar uma função lambda neste caso. Isso sim eu acho que reduz a legibilidade :-)