minimize(sum(x[i] for i in range(m)))
    for i in range(m):
        #minimize(x[i])
        st(x[i] <= b[i] - sum(A[i][j]*y[j] for j in range(n)))
        st(x[i] <= -b[i] + sum(A[i][j]*y[j] for j in range(n)))
    solve()