8嵌套循环可能组合

我已经建立一个程序，发现在这个等式中8个不同的常数：使用蛮力（A-H < = 5）8嵌套循环可能组合

a*40 +b*0 +c*3 +d*10 +e*10 +f*0 +g*9 +h*7.5 =292(+-5) 
a*4 +b*7 +c*5 +d*3 +e*0 +f*0 +g*7 +h*0 =63(+-5) 
a*0 +b*6 +c*3 +d*0 +e*0 +f*5 +g*7 +h*0 =85(+-5) 
a*175 +b*50 +c*50 +d*75 +e*75 +f*50 +g*110 +h*50 =635(+-5)

。但它需要很长时间（我知道，我知道你不需要说）我怎么能加快这个过程？

基本上这是我的代码。在现实中，我的程序有4个：

chofound=[] 
konstanta=[5,5,5,5,5,5,5,5] 
## konstanta=[10,0,5,8,2,0,4, 
for h in range(0,5): 
    for g in range(0,5): 
     for f in range(0,5): 
      for e in range(0,5): 
       for d in range(0,5): 
        for c in range(0,5): 
         for b in range(0,5): 
          for a in range(0,5): 
           hasil= a*konstanta[0]+\ 
             b*konstanta[1]+\ 
             c*konstanta[2]+\ 
             d*konstanta[3]+\ 
             e*konstanta[4]+\ 
             f*konstanta[5]+\ 
             g*konstanta[6]+\ 
             h*konstanta[7] 

           if (hasil>=(292-5) and hasil <=(292+5)): 

            asd=[a,b,c,d,e,f,g,h] 
            print ("found with config: {}".format(asd)) 
            chofound.append(asd) 


return chofound

有没有任何有效的方法来真正知道没有暴力的a-h？或者使我的代码有效运行的任何算法？

来源

2017-03-22 aji

我有什么你正在尝试做的，但猜'itertools'的关键是这样做没有真正的想法不太详细。对于现代计算机来说，5 ** 8次迭代也是微不足道的。无论你在做什么，循环本身不是瓶颈。也许这就是所有的印刷。 –

使用来自'numpy'的'array'类将允许您将其重写为矢量操作，并且应该比for循环快得多。 – Craig

292（+ - 5）是否意味着您的解决方案具有+/- 5的容差？ – Crispin

我认为这会做非常快，给大家很好的配置，以及：

import numpy as np 
from itertools import product 

konstanta = np.array([5,5,5,5,5,5,5,5]) 

configs = np.array(list(product(range(5),repeat=8))) # big array: iterating over rows is equivalent to your 8 nested for loops 
hasil = (konstanta*configs).sum(axis=1) # element wise multiplication followed by sum over rows 
good_configs = configs[(hasil>=0) & (hasil<=10)] # keep only rows where `hasil` is in desired range

来源

2017-03-22 03:35:23 Julien

这工作，但是因为我找不到任何匹配的解决方案的所有他们。我可以将“范围（5）”的步进更改为0.1吗？ – aji

您确定可以使用'np.arange（0,5,0.1）'。 – Julien

8嵌套循环可能组合

回答

相关问题