我有大约 30 个非负变量和 24 个方程，我想找出每个变量的上限和下限。保证可行的解决方案。

所以对于每个变量，我解决了两个 LP 问题，一个用于最小值，一个用于最大值

$$\begin{array}{rl} \min & c^T x\\ \textrm{s.t.} & Ax = b\\ &x\geq0 \end{array}$$

其中$c=(0, 0, \dots, 0, 1, 0, \dots, 0)^T$和$(0, 0, \dots, 0, -1, 0, \dots, 0)^T$。

因此，我总共有 60 个 LP，它们的约束完全相同，但目标略有不同。我必须为几千天的每一天的数据做同样的事情。

有没有办法以某种方式重用中间结果来加快计算速度？是否可以使用现有的 Python 或 R 包来实现？

编辑：只是为了每天添加，只有 rhs $b$不同，其他一切都保持不变。