In [1]:
fact = 1
for N in range(1,101):
    fact *= N
sum([int(x) for x in str(fact)])


Out[1]:
648

In [2]:
fact


Out[2]:
93326215443944152681699238856266700490715968264381621468592963895217599993229915608941463976156518286253697920827223758251185210916864000000000000000000000000L

In [ ]: