線形計画問題の解法の一つである「シンプレックス法」をPythonで実装してみました。 def LP(c, A, b, xb_idx): xn_idx = np.arange(A.shape[1]) xn_idx = np.delete(xn_idx, xb_idx) while True: B = A[:, xb_idx] N = A[:, xn_idx] B_inv = np.linalg.inv(B…