Estimate a log-concave probability mass function from discrete i.i.d. observations
