for n in range(100000000): ld = [couples[int(k)] for k in quinary(n)] transitions = [(ld[k][0],ld[k][1],(k+1)%2) for k in range(len(ld))] etat = (3,3,1) for t in transitions: etat = next(etat,t) if etat==(0,0,0) or etat==(0,0,1): print([mot(t) for t in transitions])