(33/16660)*(28561/31104) * log (28561/31104)+(1-(33/16660)*(28561/31104))log(172416709/172388736)