def invbinom(alpha,n,p):
    k = 0
    S = 0
    while S<alpha:
        S += binomFdp(n,p,k)
        k += 1
    return [0,k-1]