A unified approach to minimum variance unbiased estimation of a probability function belonging to an exponential family