如何使用z3 Solver解决柯克曼的女学生问题?

时间:2019-05-14 11:56:20

标签: python z3 smt

我是SMT问题的新手,并尝试使用Dennis Yurichev编写的示例复制如何使用SAT / SMT中的z3 / MK85解决Kirkman的女学生问题。但是当我尝试获取模型时(我使用的是Z3):

m["%d_%d" % (person , day)]

Python有一些错误:

  

回溯(最近通话最近一次):

     

文件“”,位于

的第1行      

__ getitem___z3_assert中的文件“ D:\ Z3 \ z3-master \ z3-master \ build \ python \ z3 \ z3.py”,行6138(False,“预期为整数,Z3声明或Z3常数”)< / p>      

文件“ D:\ Z3 \ z3-master \ z3-master \ build \ python \ z3 \ z3.py”,_ z3_assert中的第92行,引发Z3Exception(msg)z3.z3types.Z3Exception:整数,Z3声明,或预期的Z3常数

我想知道我是否对模型的表达感到困惑?我很想调试它。 此外,我想知道Z3中是否还有另一种解决方法以及z3和MK85之间的表达差异。

from z3 import *
import itertools
PERSONS , DAYS , GROUPS = 15, 7, 5
s=Solver()
tbl=[[BitVec("%d_%d"%(person , day), 16) for day in range(DAYS)] for person in range(PERSONS)]
for person in range(PERSONS):
   for day in range(DAYS):
      s.add(And(tbl[person][day]>=0, tbl[person][day] < GROUPS))
def only_one_in_pair_can_be_equal(l1, l2):
   assert len(l1)==len(l2)
   expr=[]
   for pair_eq in range(len(l1)):
             tmp=[]
             for i in range(len(l1)):
                     if pair_eq==i:
                             tmp.append(l1[i]==l2[i])
                     else:
                             tmp.append(l1[i]!=l2[i])
             expr.append(And(*tmp))
   s.add(Or(*expr))
for pair in itertools.combinations(range(PERSONS), r=2):
   only_one_in_pair_can_be_equal (tbl[pair[0]], tbl[pair[1]])
assert s.check()
m=s.model()
print("group for each person:")
print("person:"+"".join([chr(ord('A')+i)+" " for i in range(PERSONS)]))
for day in range(DAYS):
   print("day=%d:" % day)
   for person in range(PERSONS):
      print(m["%d_%d" % (person , day)])  #error
   print("")

def persons_in_group(day, group):
     rt=""
     for person in range(PERSONS):
             if m["%d_%d" % (person , day)]==group:#error
                     rt=rt+chr(ord('A')+person)
     return rt
for day in range(DAYS):
     print( "day=%d:" % day)
     for group in range(GROUPS):
             print(persons_in_group(day, group)+" ")
     print(" ")

我希望对其进行调试或以其他方式解决此问题。

1 个答案:

答案 0 :(得分:0)

这里的问题是您正在使用错误的索引访问模型。代替:

m["%d_%d" % (person , day)]

使用:

m[tbl[person][day]]

(您有两个实例,两个都需要。)