Python 使用scipy.optimize设置优化模型的不连续边界

我有一个优化程序，这是工作良好，但有一些额外的边界条件，我需要实施。目前，边界如下所示：

bnd1 = (50000, 160000)
bnd2 = (357116, 1271242)
bnd3 = (5047, 336340.29)
bnd4 = (11621, 124680.49)
bnd5 = (3655, 521087.8102)
bnd6 = (13145, 347950.8542)

优化应该能够在这些点之间取任何值，但是0也是它应该能够取的值！我是否可以引入某种if/或约束来绕过这个问题？

---

谢谢

请给我们看一些相关代码。现在我们甚至不知道我们在谈论什么样的优化（凸优化、离散优化等）。2^6=64。。。为什么不测试所有的可能性呢？除非一次优化已经花了一个小时。。。