import functools
import operator


def combine_probability_matrices(matrices):
    """given a sequence of probability matrices, combine them into a
    single matrix with sum 1.0 and return it"""

    distribution = functools.reduce(operator.mul, matrices)
    # normalize
    return distribution / distribution.sum()