如何在SymPy中判断两个ConditionSet是否不相交？

SymPy里判断两个Interval是否不相交非常直观，is_disjoint方法能直接给出正确结果：

from sympy import Interval

print(Interval(0, 2).is_disjoint(Interval(-2, -1)))  # True
print(Interval(0, 2).is_disjoint(Interval(1, 2)))     # False

但处理ConditionSet时，is_disjoint会返回错误结果——比如两个明显无交集的集合，却返回False：

from sympy import ConditionSet, Symbol, S

x = Symbol('x')
cs1 = ConditionSet(x, x > 1, S.Reals)
cs2 = ConditionSet(x, x < 1, S.Reals)

print(cs1.is_disjoint(cs2))  # 错误输出False

问题根源在于is_disjoint底层是通过判断两个集合的交集是否为EmptySet实现的，但ConditionSet的交集默认不会自动求值，只会返回未化简的Intersection对象，导致判断失效。

以下是几种可行的解决方法：

方法一：手动化简交集后判断 #

直接对两个ConditionSet的交集调用simplify，再判断是否等于EmptySet：

from sympy import ConditionSet, Symbol, S, simplify, EmptySet

x = Symbol('x')
cs1 = ConditionSet(x, x > 1, S.Reals)
cs2 = ConditionSet(x, x < 1, S.Reals)

intersection = simplify(cs1.intersect(cs2))
print(intersection == EmptySet)  # 正确输出True

方法二：直接判断条件的可满足性 #

ConditionSet的本质是满足特定条件的元素集合，因此可以直接判断两个条件的合取是否存在解——如果无解，说明两个集合不相交：

from sympy import ConditionSet, Symbol, S, satisfiable

x = Symbol('x')
cond1 = x > 1
cond2 = x < 1

# satisfiable返回False表示条件无可行解
is_overlap = satisfiable(cond1 & cond2)
print(not is_overlap)  # 正确输出True

方法三：重写ConditionSet的is_disjoint方法 #

如果需要频繁判断ConditionSet是否不相交，可以子类化ConditionSet，重写is_disjoint方法，内部优先判断条件的可满足性：

from sympy import ConditionSet, Symbol, S, satisfiable

class ImprovedConditionSet(ConditionSet):
    def is_disjoint(self, other):
        # 仅处理同变量的ConditionSet情况
        if isinstance(other, ConditionSet) and self.variables == other.variables:
            combined_cond = self.condition & other.condition
            # 合取条件不可满足则集合不相交
            return not satisfiable(combined_cond)
        # 其他情况回退到父类逻辑
        return super().is_disjoint(other)

x = Symbol('x')
cs1 = ImprovedConditionSet(x, x > 1, S.Reals)
cs2 = ImprovedConditionSet(x, x < 1, S.Reals)
print(cs1.is_disjoint(cs2))  # 正确输出True

三种方法各有适用场景：方法一适合快速修改现有代码；方法二更高效，直接针对条件逻辑判断；方法三适合长期频繁使用的场景，封装后无需重复写判断逻辑。

内容的提问来源于stack exchange，提问作者marwan