import math


def approximate_product(pair_one, pair_two):
    """
    Calculates the product of two estimates and approximates its margin of error.

    Follows the U.S. Census Bureau's `official guidelines`_.

    Args:
        pair_one (list): a two-item sequence with a U.S. Census
            bureau estimate and its margin of error.

        pair_two (list): a two-item sequence with a U.S. Census
            bureau estimate and its margin of error.

    Returns:
        A two-item sequence containing the estimate and its approximate margin of error.

        (61393366, 202289)

    Examples:
        >>> approximate_product((74506512, 228238), (0.824, 0.001))
        (61393366, 202289)

    ... _official guidelines:
        https://www.documentcloud.org/documents/6177941-Acs-General-Handbook-2018-ch08.html#document/p8
    """
    # Pull out the values
    estimate_one, moe_one = pair_one
    estimate_two, moe_two = pair_two

    # Approximate the product
    product_estimate = estimate_one * estimate_two

    # Approximate the margin of error
    squared_product_moe = (estimate_one**2 * moe_two**2) + (estimate_two**2 * moe_one**2)
    product_moe = math.sqrt(squared_product_moe)

    # Return the results
    return product_estimate, product_moe