z3给出了一个令人惊讶的答案

Question

我试过这个：

import z3
n,x,y,z = z3.Ints('n x y z')
s = z3.Solver()
s.add(4/n==1/x+1/y+1/z)
s.add(x>0)
s.add(n>0)
s.add(y>0)
s.add(z>0)
s.check()
s.model()

我得到：

[x = 1, n = 2, z = 3, y = 1, div0 = [(1, 1) → 1, (4, 2) → 2, else → 0], mod0 = [(1, 3) → 1, else → 0]]

但是4/2不等于1/1+1/1+1/3

我做错了什么？

Answer 1

您已将 n, x, y, z 声明为整数。所以除法是作为一个整数完成的，给你 1/1 = 1 和 1/3 = 0；因此满足等式 2=2.

显而易见的事情是对这个问题使用 Real 值。将声明更改为：

n,x,y,z = z3.Reals('n x y z')

产生：

[z = 1, y = 1, n = 4/3, x = 1]

这确实满足了你提出的方程式。

以防您确实希望 n, x, y, z 为整数；那么你应该在除法之前将它们转换为实数，如下所示：

import z3
n,x,y,z = z3.Ints('n x y z')
s = z3.Solver()
s.add(4/z3.ToReal(n)==1/z3.ToReal(x)+1/z3.ToReal(y)+1/z3.ToReal(z))
s.add(x>0)
s.add(n>0)
s.add(y>0)
s.add(z>0)
print(s.check())
print(s.model())

这会打印：

sat
[n = 4, x = 3, z = 3, y = 3]

再次满足您的限制。

z3给出了一个令人惊讶的答案

z3 giving a surprising answer

python

z3

z3py