def theoritical_loss_rate(rho, K):
    """
                        k
               (1 - p) p
      CLR =    ----------
                     k + 1
                1 - p

    """
    return (1-rho) * rho ** K / (1 - rho ** (K+1))