在固体力学中,我经常使用Python编写类似以下代码的程序:
for i in range(3):
for j in range(3):
for k in range(3):
for l in range(3):
# do stuff
我经常这样做,开始怀疑是否有更简洁的方法。当前代码的缺点是:如果我符合PEP8
,那么我不能超过每行79个字符的限制,而且左右的空间不多,尤其是如果这又是类的函数中的代码。在固体力学中,我经常使用Python编写类似以下代码的程序:
for i in range(3):
for j in range(3):
for k in range(3):
for l in range(3):
# do stuff
我经常这样做,开始怀疑是否有更简洁的方法。当前代码的缺点是:如果我符合PEP8
,那么我不能超过每行79个字符的限制,而且左右的空间不多,尤其是如果这又是类的函数中的代码。使用嵌套的for循环,您基本上是在尝试创建所谓的输入可迭代对象的(笛卡尔)积,这就是来自itertools
模块的product
函数的作用。
>>> list(product(range(3),repeat=4))
[(0, 0, 0, 0), (0, 0, 0, 1), (0, 0, 0, 2), (0, 0, 1, 0), (0, 0, 1, 1),
(0, 0, 1, 2), (0, 0, 2, 0), (0, 0, 2, 1), (0, 0, 2, 2), (0, 1, 0, 0),
...
而在您的代码中,您可以这样做:
for i,j,k,l in product(range(3),repeat=4):
#do stuff
def product(*args, repeat=1):
# product('ABCD', 'xy') --> Ax Ay Bx By Cx Cy Dx Dy
# product(range(2), repeat=3) --> 000 001 010 011 100 101 110 111
pools = [tuple(pool) for pool in args] * repeat
result = [[]]
for pool in pools:
result = [x+[y] for x in result for y in pool]
for prod in result:
yield tuple(prod)
product(range(3),range(4),['a','b','c'] ,some_other_iterable)
。如果嵌套循环比这更复杂(例如包含一些if / else逻辑),则定义自己的生成函数是下一个最好的选择。在这里,我更喜欢@SergeBallesta版本的生成器,因为它更简单易懂。 - PeterEfor i,j,k,l in product(range(3),repeat=4)
中找到最简单的方式,我还添加了product
函数的等效功能来回答你,你可以查看一下。 - Mazdakproduct()
仍然可以使用。 - PeterEproduct(range(100), repeat=4)
来对比这两个版本。product
花费了我12秒钟,而长版本只需要5秒钟。 - PascalVKooten使用 itertools.product
的想法很好。这里有一种更通用的方法,可以支持不同大小的范围。
from itertools import product
def product_of_ranges(*ns):
for t in product(*map(range, ns)):
yield t
for i, j, k in product_of_ranges(4, 2, 3):
# do stuff
虽然这样会让你需要使用生成器函数,但是它不会更为简洁,但至少你不必被PEP8所困扰:
def tup4(n):
for i in range(n):
for j in range(n):
for k in range(n):
for l in range(n):
yield (i, j, k, l)
for (i, j, k, l) in tup4(3):
# do your stuff
(在Python 2.x中,你应该在生成器函数中使用xrange
而不是range
)。
编辑:
当金字塔的深度已知时,上述方法应该没问题。但是你也可以创建一个通用的生成器,而无需任何外部模块:
def tup(n, m):
""" Generate all different tuples of size n consisting of integers < m """
l = [ 0 for i in range(n)]
def step(i):
if i == n : raise StopIteration()
l[i] += 1
if l[i] == m:
l[i] = 0
step(i+ 1)
while True:
yield tuple(l)
step(0)
for (l, k, j, i) in tup(4, 3):
# do your stuff
我使用了(l, k, j, i)
,因为在上面的生成器中,第一个索引首先变化。
这是等效的:
for c in range(3**4):
i = c // 3**3 % 3
j = c // 3**2 % 3
k = c // 3**1 % 3
l = c // 3**0 % 3
print(i,j,k,l)
如果你经常这样做,考虑使用一个通用的生成器:
def nestedLoop(n, l):
return ((tuple((c//l**x%l for x in range(n-1,-1,-1)))) for c in range(l**n))
for (a,b,c,d) in nestedLoop(4,3):
print(a,b,c,d)
c // 3**3 % 3
、c // 3**2 % 3
等等。 ;-)
- martineauyield
一个元组而不是返回生成器表达式来进一步缩短它,同时通过将reversed(range(n))
更改为range(n)[::-1]
(或range(n,-1,-1)
)来进一步缩短它。 - martineau
numpy.einsum()
,参考使用NumPy实现快速张量旋转。 - jfs