def symm_area(col, n):
    """
    returns n + (n - 1) + ... + (n - col + 1)
    i.e., the number of matrix elements below and including the diagonal and
    from column 0 to column `col`
    """
    return col * (2 * n - col + 1) // 2