我试图用Z3解决以下问题:
拥有四个运算符(+
,-
,*
,/
),决定哪些运算符可以替换为以下表达式,以使其成立:< / p>
(((((1 <op1> 2) <op2> 3) <op3> 4) <op4> 5) <op5> 6) = 35
我需要打印所有有效答案。这是Z3语言的示例程序:
(declare-datatypes () ((operator (Plus) (Minus) (Mult) (Divid))))
(define-fun myf ((x Int) (z operator) (y Int)) Int
(ite (= z Plus) (+ x y)
(ite (= z Minus) (- x y)
(ite (= z Mult) (* x y)
(div x y)))))
(assert
(exists ((b1 operator) (b2 operator) (b3 operator) (b4 operator) (b5 operator))
(= (myf(myf(myf(myf(myf 1 b1 2) b2 3) b3 4) b4 5) b5 6) 35)
)
)
(check-sat)
(get-model)
(exit)
它不会打印所有解决方案。我知道要打印我需要在每次迭代like that answer recommends之后向求解器添加约束的所有解决方案,但我无法用C#来解决这个问题。
这是我目前的代码(C#7):
using (var context = new Context())
{
var @operator = context.MkEnumSort("operator", "Plus", "Minus", "Mult", "Div");
var plus = @operator.Consts[0];
var minus = @operator.Consts[1];
var mult = @operator.Consts[2];
IntExpr myf(IntExpr x, Expr z, IntExpr y) =>
(IntExpr)context.MkITE(context.MkEq(z, plus), context.MkAdd(x, y),
context.MkITE(context.MkEq(z, minus), context.MkSub(x, y),
context.MkITE(context.MkEq(z, mult), context.MkMul(x, y),
context.MkDiv(x, y))));
var solver = context.MkSolver();
var b1 = context.MkConst("b1", @operator);
var b2 = context.MkConst("b2", @operator);
var b3 = context.MkConst("b3", @operator);
var b4 = context.MkConst("b4", @operator);
var b5 = context.MkConst("b5", @operator);
solver.Assert(
context.MkExists(
new[] { b1, b2, b3, b4, b5 },
context.MkEq(
myf(
myf(
myf(
myf(
myf(
context.MkInt(1),
b1,
context.MkInt(2)),
b2,
context.MkInt(3)),
b3,
context.MkInt(4)),
b4,
context.MkInt(5)),
b5,
context.MkInt(6)),
context.MkInt(35))));
while (Status.SATISFIABLE == solver.Check())
{
var operators = new[] { b1, b2, b3, b4, b5 };
var model = solver.Model;
var values = operators.Select(o => model.Eval(o, true)); // That doesn't return the right values
Console.WriteLine(model);
Console.WriteLine(string.Join(" ", values));
solver.Add(context.MkOr(
operators.Select(o => context.MkNot(context.MkEq(o, model.Eval(o, true)))))); // That's supposed to work, but it doesn't
}
我在访问exists
范围内的变量时遇到问题:似乎model.Eval(b1, true)
返回一些值,但不值解算器决定在当前迭代中使用。甚至询问这些常量的值会进一步污染求解器范围(例如,我可以在模型输出中一直看到这些常量):
(define-fun b3!2 () operator
Mult)
(define-fun b2!3 () operator
Mult)
(define-fun b5!0 () operator
Minus)
(define-fun b1!4 () operator
Plus)
(define-fun b4!1 () operator
Plus)
# ^ ^ ^ ^ these seems like the proper values
(define-fun b1 () operator
Minus)
(define-fun b2 () operator
Minus)
(define-fun b3 () operator
Minus)
(define-fun b4 () operator
Minus)
(define-fun b5 () operator
Minus)
# ^ ^ ^ ^ and I don't know why it prints these
如何修复程序以添加对正确值的约束并且不会污染范围?
答案 0 :(得分:4)
不幸的是我无法访问C#来试试这个,但我很好奇为什么你需要打电话给mkExist
?由于您已经通过调用b1
创建了变量b6
.. MkConst
,因此您只需使用它们来断言您的约束,然后在执行“全部”时反驳模型-sat“循环,没有对mkExist
的任何调用以及没有在循环中调用new
来创建新的。{/ p>