def homo_line(a, b):
    """
    Return the homogenous equation of a line passing through a and b,
    i.e. [ax, ay, 1] cross [bx, by, 1]
    """
    return (a[1] - b[1], b[0] - a[0], a[0] * b[1] - a[1] * b[0])