可以使用Z3py计算Kauffman括号吗？

Question

我正在尝试使用Z3py计算三叶草结的Kauffman支架。到目前为止，我有以下代码：

a, b, c, d, e, f, A, B = Ints('a b c d e f A B')

delta = Function('delta', IntSort(), IntSort(), IntSort())
def X(a,b,c,d):
    return A*delta(a,d)*delta(b,c)+B*delta(a,b)*delta(c,d)
Trefoil = X(a,d,b,e)*X(e,b,f,c)*X(c,f,d,a)  
print simplify(simplify(Trefoil, som= True),mul_to_power=True)

使用此代码我将获得以下输出：

A**3·
delta(a, e)·
delta(b, f)·
delta(c, a)·
delta(d, b)·
delta(e, c)·
delta(f, d) +
A**2·
B·
delta(a, d)·
delta(b, e)·
delta(b, f)·
delta(c, a)·
delta(e, c)·
delta(f, d) +
A**2·
B·
delta(a, e)·
delta(c, a)·
delta(d, b)·
delta(e, b)·
delta(f, c)·
delta(f, d) +
A·
B**2·
delta(a, d)·
delta(b, e)·
delta(c, a)·
delta(e, b)·
delta(f, c)·
delta(f, d) +
A**2·
B·
delta(a, e)·
delta(b, f)·
delta(c, f)·
delta(d, a)·
delta(d, b)·
delta(e, c) +
A·
B**2·
delta(a, d)·
delta(b, e)·
delta(b, f)·
delta(c, f)·
delta(d, a)·
delta(e, c) +
A·
B**2·
delta(a, e)·
delta(c, f)·
delta(d, a)·
delta(d, b)·
delta(e, b)·
delta(f, c) +
B**3·
delta(a, d)·
delta(b, e)·
delta(c, f)·
delta(d, a)·
delta(e, b)·
delta(f, c)

现在我需要应用规则：

delta(a,b)*delta(b,c) = delta(a,c)

并使用此规则简化输出。拜托，你能告诉我怎么做。非常感谢。

Answer 1

您可以使用“替代”功能替换他人的子词。 但是，在您的示例中，没有子项delta（a，b）或delta（a，c） 所以替换将无效。请记住，Z3确实如此 不在运营商的关联性或可交换性下进行替换。 所以替代（x * y * z，（x * z，u））并不简化为u * y。

Answer 2

使用Reduce的解决方案：

如果可以使用Z3PY做同样的事情，请告诉我。