我想我已经想出了一个解决方案,但如果这个错误,请纠正我。我也很想看看是否有更优雅的解决方案。
首先,我想出了组合的总数。所有没有替换的组合都等于 n!/r!(nr)! 并且替换等于 (m+s-1)!/s!(m-1)! 其中 m 和 n 是可供选择的项目数,r 和 s 是您实际选择的项目数。因为我知道每个组合中我想要的总项目(我们称之为上限),我找到了无替换类型(n=0)的 0 和替换类型的“上限”(m=3)的组合数并将这些数字相乘。然后,将无替换类型 (n=1) 的 1 的组合数乘以替换类型 (m=2) 的组合“cap-1”。这样做直到您最终添加了无替换类型(n = 3)乘以替换类型(m = 0)的“上限”的组合(感谢@André Nicolas)。
import itertools
from math import factorial as fact
norep = ['A','B','C']
rep = ['1','2','3']
cap = 3 #length of combinations, e.g. cap=3, combo1=123,combo2=A12,etc
combos = 0
for i in range(cap+1):
combnorep = fact(len(norep))/(fact(cap-i)*fact(len(norep)-(cap-i)))
combrep = fact(len(rep)+i-1)/(fact(i)*fact(len(rep)-1))
combos = combos + combnorep*combrep
print combos
对于此示例,组合数为 38。接下来,我想打印所有组合。为此,我确定了所有替换的组合、所有无替换的组合以及两者的任意组合,例如 n=0,m=3;n=1,m=2;等。这就是我想出的:
for i in range(cap+1):
norepcomb = [j for j in itertools.combinations(norep,i)]
repcomb = [k for k in itertools.combinations_with_replacement(rep,cap-i)]
for l in itertools.product(norepcomb,repcomb):
print list(itertools.chain.from_iterable(l))
要包括none
,我只会将none
替换组合包括在我的列表中。我希望对此有任何反馈,特别是如果有更好的解决方案,或者如果它不像我认为的那样有效。谢谢!