def func_mv_ppf(x, variant_params, p):
    """Function CDF of max of pareto random variables for root-finding."""
    cdf = 1.0
    for a, b in variant_params:
        cdf *= 1 - (b / x) ** a
    return cdf - p
