寻找所有不等式解的优化方法

Question

我正在寻找不平等解决方案的数量 c > (a+(b^2)−1)/(a−1)

受约束2<=a<=A, 1<=b<=B, 1<=c<=C。

到目前为止，我使用的方法是对 a 使用带外部循环的嵌套循环，对b使用内部循环。由于A可以和10^9一样大，因此我正在尝试寻找优化方法的方法。

任何建议将不胜感激。

Answer 1

使用numpy进行矢量化操作。您可以执行以下操作->

import numpy as np

a = np.arange(2,A)
b = np.arange(1,B)
a1 = a.repeat(len(b), axis=0)
b1 = b.repeat(len(a), axis=0)

rhs = (a1 + (b1**2) -1) / (a1 - 1)