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