python-3.x - 求解变量取二进制值的线性方程组

Question

我正在尝试求解变量取二进制值的线性方程组

import numpy as np
import numpy.linalg as LA
import scipy.optimize as optimize

A = np.array([[1/5, 1/2, 1/3], [1/4, 1/5, 1/6], [1/6, 2/9, 3/10]])
b = np.array([1, 1, 1])
x = LA.solve(A, b)

我的矩阵 A 大约有 400 列和 3000 行。由于篇幅限制，只发布了几列。我的输出向量始终为 1。我所有的变量只能采用二进制值 - 0 或 1