我们假设x,y,z必须是正整数。否则就有无穷多的解。
从x和z计算y
这里有一个方法应该比你的方法快500倍,因为它不在y上迭代:for z in range(500):
for x in range(500):
fifty_y = 1950 - (x * 80) - (z * 65)
if fifty_y >= 0 and (fifty_y % 50) == 0:
y = fifty_y // 50
print("x " + str(x) + " z " + str(z) + " y " + str(y))
通过迭代x、y和z,你基本上是在黑暗中拍摄,希望它落在1950上。
但是你知道50 * y = 1950 - x * 80 - z * 65,所以你可以直接从x和z计算y。
50 * y应该是正的,如果y是整数,它应该可以被50整除。
限制x和z
如果我们希望y是阳性的话,range(500)对于x和z来说太大了。
range(1950 // 65 + 1)应该足以满足z。
知道z,range((1950 - 65 * z)// 80 + 1)就足够了。
作为奖励,我们确信50 * y是阳性的,并且我们可以删除一个测试:for z in range(1950 // 65 + 1):
for x in range((1950 - 65 * z) // 80 + 1):
fifty_y = 1950 - (x * 80) - (z * 65)
if (fifty_y % 50) == 0:
y = fifty_y // 50
print("x " + str(x) + " z " + str(z) + " y " + str(y))
与沃尔夫拉姆阿尔法
通过在wolfram alpha中键入方程式,我们得到:Integer solution: y = 13 n + x, z = -10 n - 2 x + 30, n element Z
太好了!对于任何x,我们只需要选择n,这样y和z都是阳性的。不再需要if。此代码重复42次以显示42个解决方案:for x in range(1950 // 80 + 1):
for n in range(- (x // 13), (30 - 2 * x) // 10 + 1):
z = -10 * n - 2 * x + 30
y = 13 * n + x
print("x " + str(x) + " z " + str(z) + " y " + str(y))
它几乎可以排成一行:print [(x, 13 * n + x, 30 - 10 * n - 2 * x) for x in range(25) for n in range(-(x // 13), (30 - 2 * x) // 10 + 1)]
这段代码比原来的代码快300万倍:)