Formule mathématique en python

def fact(a):
    pass

def C(n,p):
    pass

if __name__=='__main__':
    res=0
    for k in range(c,b+1):
        res+=C(b,k)*(1.0/a)**k*(1-1.0/a)**(b-k)

    print res

#!/usr/bin/python
# -*- coding: iso-8859-1 -*-

zmin=int(raw_input("zombies min="))
zmax=int(raw_input("zombies max="))
defville=int(raw_input("défense de la ville="))
c=int(raw_input("citoyens ="))

#retrait de la défense de la ville
zmaxen=zmax-defville
zminen=zmin-defville

#impossible de subir une attaque négative
if zmaxen<0:
	zmaxen=0

if zminen<0:
	zminen=0

#calcule des  chances minimuns et maximums
percentmax=((float(c)-1)/c)**zmaxen
percentmin=((float(c)-1)/c)**zminen

#écriture du résultat
print "pourcentage de survie si l attaque est au plus fort =", percentmax
print "pourcentage de survie si l attaque est au moins fort =", percentmin

def C(n,p) :
    return fact(n)/fact(p)/fact(n-p)