我想你明白你的需求是一个非常大的列表。即使我们只考虑2个变量之间的d-分离（由一组节点约束）。
无论如何，使用pyAgrum（https://agrum.org）可以很容易地做到这一点。

import itertools
import pyAgrum as gum

# create a BN
bn=gum.fastBN("A->B<-C->D->E->F;B->E<-G");

# print the indepency model by testing d-separations

# how to iterate for each subset of an interable
def powerset(iterable):
    """
    powerset([1,2,3]) --> () (1,) (2,) (3,) (1,2) (1,3) (2,3) (1,2,3)
    """
    xs = list(iterable)
    # note we return an iterator rather than a list
    return itertools.chain.from_iterable(itertools.combinations(xs,n) for n in range(len(xs)+1))

# testing every d-separation
for i in bn.names():
    for j in bn.names()-{i}:
        for k in powerset(bn.names()-{i,j}):
            if bn.isIndependent(i,j,k):
                print(f"{i} indep {j} given {k}")

结果（在笔记本中）：

我做了一个软件包，准确地实现了你想要的功能。真巧啊！

pip install cdpi

from cdpi import pattern
ptn = pattern()
ptn.add_edges([
    ('A', 'B'),
    ('C', 'B'),
    ('B', 'D'),
    ('E', 'D'),
    ('F', 'D'),
    ('D', 'G')
])
ptn.draw()

[out]图形可视化也支持！但10声誉限制...😂我不能上传图像...

ds = ptn.get_all_d_separation() # ds is 2-d dictionary
ds['F']['A'] #example

[出]

{(),
 ('B',),
 ('B', 'C'),
 ('B', 'D'),
 ('B', 'D', 'C'),
 ('B', 'E'),
 ('B', 'E', 'C'),
 ('B', 'E', 'D'),
 ('B', 'E', 'D', 'C'),
 ('B', 'E', 'D', 'G'),
 ('B', 'E', 'G'),
 ('B', 'E', 'G', 'C'),
 ('B', 'G'),
 ('B', 'G', 'C'),
 ('C',),
 ('E',),
 ('E', 'C'),
 ('G', 'B', 'D'),
 ('G', 'B', 'D', 'C'),
 ('G', 'E', 'C', 'B', 'D')}

如果要查看整个代码或其工作原理，请参见here的d-separation part