Comment filtrer des éléments de produit cartésien après des conditions de commande spécifiques

Je dois générer plusieurs réactions avec différentes variables. Ils ont 3 éléments. Appelons-les B code>, s code> et h code>. Et ils commencent tous avec b1 code>. s code> peut être ajouté à l'élément s'il y a au moins un B code>. Donc, il peut être b1s1 code> ou b2s2 code> ou b2s1 code> etc ... mais pas b1s2 code>. Il en va de même pour h code>. b1s1h1 code> ou b2s2h1 code> ou b4s1h1 code> mais jamais b2s2h3 code>. La variation finale serait b5s5h5 code>. J'ai essayé avec itheroTools.Product code>. Mais je ne sais pas comment vous débarrasser des éléments qui ne correspondent pas à mon état et comment ajouter l'élément suivant. Voici mon code:

import itertools

a = list(itertools.product([1, 2, 3, 4], repeat=4))

#print (a)

met = open('random_dat.dat', 'w')
met.write('Reactions')
met.write('\n')
for i in range(1,256): 
    met.write('\n')
    met.write('%s: B%sS%sH%s -> B%sS%sH%s' %(i, a[i][3], a[i][2], a[i][1], a[i][3], a[i][2], a[i][1]))
    met.write('\n')
met.close()