1067 lines
35 KiB
Python
1067 lines
35 KiB
Python
# Copyright 2010-2017 Google
|
|
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
# you may not use this file except in compliance with the License.
|
|
# You may obtain a copy of the License at
|
|
#
|
|
# http://www.apache.org/licenses/LICENSE-2.0
|
|
#
|
|
# Unless required by applicable law or agreed to in writing, software
|
|
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
# See the License for the specific language governing permissions and
|
|
# limitations under the License.
|
|
"""Propose a natural linear API on top of cp_model_pb2 python proto.
|
|
|
|
This file implements a easy to use API on top of the cp_model_pb2 API
|
|
defined in ../ .
|
|
"""
|
|
|
|
from __future__ import absolute_import
|
|
from __future__ import division
|
|
from __future__ import print_function
|
|
|
|
import collections
|
|
import numbers
|
|
from six import iteritems
|
|
|
|
from ortools.sat import cp_model_pb2
|
|
from ortools.sat import sat_parameters_pb2
|
|
from ortools.sat import pywrapsat
|
|
|
|
# The classes below allow linear expressions to be expressed naturally with the
|
|
# usual arithmetic operators +-*/ and with constant numbers, which makes the
|
|
# python API very intuitive. See cp_model_test.py for examples.
|
|
|
|
INT_MIN = -9223372036854775808 # hardcoded to be platform independent.
|
|
INT_MAX = 9223372036854775807
|
|
|
|
|
|
# Cp Solver status (exported to avoid importing cp_model_cp2).
|
|
UNKNOWN = cp_model_pb2.UNKNOWN
|
|
MODEL_INVALID = cp_model_pb2.MODEL_INVALID
|
|
MODEL_SAT = cp_model_pb2.MODEL_SAT
|
|
MODEL_UNSAT = cp_model_pb2.MODEL_UNSAT
|
|
OPTIMAL = cp_model_pb2.OPTIMAL
|
|
|
|
|
|
def AssertIsInt64(x):
|
|
"""Asserts that x is integer and x is in [min_int_64, max_int_64]."""
|
|
if not isinstance(x, numbers.Integral):
|
|
raise TypeError('Not an integer: %s' % x)
|
|
if x < INT_MIN or x > INT_MAX:
|
|
raise OverflowError('Does not fit in an int64: %s' % x)
|
|
|
|
|
|
def AssertIsBoolean(x):
|
|
"""Asserts that x is 0 or 1."""
|
|
if not isinstance(x, numbers.Integral) or x < 0 or x > 1:
|
|
raise TypeError('Not an boolean: %s' % x)
|
|
|
|
|
|
def CapInt64(v):
|
|
if v > INT_MAX:
|
|
return INT_MAX
|
|
if v < INT_MIN:
|
|
return INT_MIN
|
|
return v
|
|
|
|
|
|
def CapSub(x, y):
|
|
"""Saturated arithmetics. Returns x - y truncated to the int64 range."""
|
|
if not isinstance(x, numbers.Integral):
|
|
raise TypeError('Not integral: ' + str(x))
|
|
if not isinstance(y, numbers.Integral):
|
|
raise TypeError('Not integral: ' + str(y))
|
|
AssertIsInt64(x)
|
|
AssertIsInt64(y)
|
|
if y == 0:
|
|
return x
|
|
if x == y:
|
|
if x == INT_MAX or x == INT_MIN:
|
|
raise OverflowError(
|
|
'Integer NaN: subtracting INT_MAX or INT_MIN to itself')
|
|
return 0
|
|
if x == INT_MAX or x == INT_MIN:
|
|
return x
|
|
if y == INT_MAX:
|
|
return INT_MIN
|
|
if y == INT_MIN:
|
|
return INT_MAX
|
|
return CapInt64(x - y)
|
|
|
|
|
|
def DisplayBounds(bounds):
|
|
"""Displays a flattened list of intervals."""
|
|
out = ''
|
|
for i in range(0, len(bounds), 2):
|
|
if i != 0:
|
|
out += ', '
|
|
if bounds[i] == bounds[i + 1]:
|
|
out += str(bounds[i])
|
|
else:
|
|
out += str(bounds[i]) + '..' + str(bounds[i + 1])
|
|
return out
|
|
|
|
|
|
def ShortName(model, i):
|
|
"""Returns a short name of an integer variable, or its negation."""
|
|
if i < 0:
|
|
return 'Not(%s)' % ShortName(model, -i - 1)
|
|
v = model.variables[i]
|
|
if v.name:
|
|
return v.name
|
|
elif len(v.domain) == 2 and v.domain[0] == v.domain[1]:
|
|
return str(v.domain[0])
|
|
else:
|
|
return '[%s]' % DisplayBounds(v.domain)
|
|
|
|
|
|
class IntegerExpression(object):
|
|
"""Holds an integer expression."""
|
|
|
|
def GetVarValueMap(self):
|
|
"""Scan the expression, and return a list of (var_coef_map, constant)."""
|
|
coeffs = collections.defaultdict(int)
|
|
constant = 0
|
|
to_process = [(self, 1)]
|
|
while to_process: # Flatten to avoid recursion.
|
|
expr, coef = to_process.pop()
|
|
if isinstance(expr, _ProductCst):
|
|
to_process.append((expr.Expression(), coef * expr.Coefficient()))
|
|
elif isinstance(expr, _SumArray):
|
|
for e in expr.Array():
|
|
to_process.append((e, coef))
|
|
constant += expr.Constant() * coef
|
|
elif isinstance(expr, IntVar):
|
|
coeffs[expr] += coef
|
|
elif isinstance(expr, _NotBooleanVariable):
|
|
raise TypeError('Cannot interpret literals in a integer expression.')
|
|
else:
|
|
raise TypeError('Unrecognized integer expression: ' + str(expr))
|
|
|
|
return coeffs, constant
|
|
|
|
def __hash__(self):
|
|
return object.__hash__(self)
|
|
|
|
def __add__(self, expr):
|
|
return _SumArray([self, expr])
|
|
|
|
def __radd__(self, arg):
|
|
return _SumArray([self, arg])
|
|
|
|
def __sub__(self, expr):
|
|
return _SumArray([self, -expr])
|
|
|
|
def __rsub__(self, arg):
|
|
return _SumArray([-self, arg])
|
|
|
|
def __mul__(self, arg):
|
|
if isinstance(arg, numbers.Integral):
|
|
if arg == 1:
|
|
return self
|
|
AssertIsInt64(arg)
|
|
return _ProductCst(self, arg)
|
|
elif isinstance(arg, IntegerExpression):
|
|
return _Product(self, arg)
|
|
else:
|
|
raise TypeError('Not an integer expression: ' + str(arg))
|
|
|
|
def __rmul__(self, arg):
|
|
AssertIsInt64(arg)
|
|
if arg == 1:
|
|
return self
|
|
return _ProductCst(self, arg)
|
|
|
|
def __div__(self, _):
|
|
raise NotImplementedError('IntegerExpression.__div__')
|
|
|
|
def __truediv__(self, _):
|
|
raise NotImplementedError('IntegerExpression.__truediv__')
|
|
|
|
def __mod__(self, _):
|
|
raise NotImplementedError('IntegerExpression.__mod__')
|
|
|
|
def __neg__(self):
|
|
return _ProductCst(self, -1)
|
|
|
|
def __eq__(self, arg):
|
|
if isinstance(arg, numbers.Integral):
|
|
AssertIsInt64(arg)
|
|
return BoundIntegerExpression(self, [arg, arg])
|
|
else:
|
|
return BoundIntegerExpression(self - arg, [0, 0])
|
|
|
|
def __ge__(self, arg):
|
|
if isinstance(arg, numbers.Integral):
|
|
AssertIsInt64(arg)
|
|
return BoundIntegerExpression(self, [arg, INT_MAX])
|
|
else:
|
|
return BoundIntegerExpression(self - arg, [0, INT_MAX])
|
|
|
|
def __le__(self, arg):
|
|
if isinstance(arg, numbers.Integral):
|
|
AssertIsInt64(arg)
|
|
return BoundIntegerExpression(self, [INT_MIN, arg])
|
|
else:
|
|
return BoundIntegerExpression(self - arg, [INT_MIN, 0])
|
|
|
|
def __lt__(self, arg):
|
|
if isinstance(arg, numbers.Integral):
|
|
AssertIsInt64(arg)
|
|
if arg == INT_MIN:
|
|
raise ArithmeticError('< INT_MIN is not supported')
|
|
return BoundIntegerExpression(self, [INT_MIN, CapInt64(arg - 1)])
|
|
else:
|
|
return BoundIntegerExpression(self - arg, [INT_MIN, -1])
|
|
|
|
def __gt__(self, arg):
|
|
if isinstance(arg, numbers.Integral):
|
|
AssertIsInt64(arg)
|
|
if arg == INT_MAX:
|
|
raise ArithmeticError('> INT_MAX is not supported')
|
|
return BoundIntegerExpression(self, [CapInt64(arg + 1), INT_MAX])
|
|
else:
|
|
return BoundIntegerExpression(self - arg, [1, INT_MAX])
|
|
|
|
def __ne__(self, arg):
|
|
if isinstance(arg, numbers.Integral):
|
|
AssertIsInt64(arg)
|
|
if arg == INT_MAX:
|
|
return BoundIntegerExpression(self, [INT_MIN, INT_MAX - 1])
|
|
elif arg == INT_MIN:
|
|
return BoundIntegerExpression(self, [INT_MIN + 1, INT_MAX])
|
|
else:
|
|
return BoundIntegerExpression(
|
|
self, [INT_MIN,
|
|
CapInt64(arg - 1),
|
|
CapInt64(arg + 1), INT_MAX])
|
|
else:
|
|
return BoundIntegerExpression(self - arg, [INT_MIN, -1, 1, INT_MAX])
|
|
|
|
|
|
class _ProductCst(IntegerExpression):
|
|
"""Represents the product of a IntegerExpression by a constant."""
|
|
|
|
def __init__(self, expr, coef):
|
|
AssertIsInt64(coef)
|
|
if isinstance(expr, _ProductCst):
|
|
self.__expr = expr.Expression()
|
|
self.__coef = expr.Coefficient() * coef
|
|
else:
|
|
self.__expr = expr
|
|
self.__coef = coef
|
|
|
|
def __str__(self):
|
|
if self.__coef == -1:
|
|
return '-' + str(self.__expr)
|
|
else:
|
|
return '(' + str(self.__coef) + ' * ' + str(self.__expr) + ')'
|
|
|
|
def __repr__(self):
|
|
return 'ProductCst(' + repr(self.__expr) + ', ' + repr(self.__coef) + ')'
|
|
|
|
def Coefficient(self):
|
|
return self.__coef
|
|
|
|
def Expression(self):
|
|
return self.__expr
|
|
|
|
|
|
class _SumArray(IntegerExpression):
|
|
"""Represents the sum of a list of IntegerExpression and a constant."""
|
|
|
|
def __init__(self, array):
|
|
self.__array = []
|
|
self.__constant = 0
|
|
for x in array:
|
|
if isinstance(x, numbers.Integral):
|
|
AssertIsInt64(x)
|
|
self.__constant += x
|
|
elif isinstance(x, IntegerExpression):
|
|
self.__array.append(x)
|
|
else:
|
|
raise TypeError('Not an integer expression: ' + str(x))
|
|
|
|
def __str__(self):
|
|
if self.__constant == 0:
|
|
return '({})'.format(' + '.join(map(str, self.__array)))
|
|
else:
|
|
return '({} + {})'.format(' + '.join(map(str, self.__array)),
|
|
self.__constant)
|
|
|
|
def __repr__(self):
|
|
return 'SumArray({}, {})'.format(', '.join(map(repr, self.__array)),
|
|
self.__constant)
|
|
|
|
def Array(self):
|
|
return self.__array
|
|
|
|
def Constant(self):
|
|
return self.__constant
|
|
|
|
|
|
class IntVar(IntegerExpression):
|
|
"""Represents a IntegerExpression containing only a single variable."""
|
|
|
|
def __init__(self, model, bounds, name, is_present_index=None):
|
|
"""See CpModel.NewIntVar and .NewOptionalIntVar below."""
|
|
self.__model = model
|
|
self.__index = len(model.variables)
|
|
self.__var = model.variables.add()
|
|
self.__var.domain.extend(bounds)
|
|
self.__var.name = name
|
|
self.__negation = None
|
|
if is_present_index is not None:
|
|
self.__var.enforcement_literal.append(is_present_index)
|
|
|
|
def Index(self):
|
|
return self.__index
|
|
|
|
def __str__(self):
|
|
return self.__var.name
|
|
|
|
def __repr__(self):
|
|
return '%s(%s)' % (self.__var.name, DisplayBounds(self.__var.domain))
|
|
|
|
def Not(self):
|
|
for bound in self.__var.domain:
|
|
if bound < 0 or bound > 1:
|
|
raise TypeError('Cannot call Not on a non boolean variable: %s' % self)
|
|
if not self.__negation:
|
|
self.__negation = _NotBooleanVariable(self)
|
|
return self.__negation
|
|
|
|
|
|
class _NotBooleanVariable(IntegerExpression):
|
|
|
|
def __init__(self, boolvar):
|
|
self.__boolvar = boolvar
|
|
|
|
def Index(self):
|
|
return -self.__boolvar.Index() - 1
|
|
|
|
def Not(self):
|
|
return self.__boolvar
|
|
|
|
|
|
class _Product(IntegerExpression):
|
|
"""Represents the product of two IntegerExpressions."""
|
|
|
|
def __init__(self, left, right):
|
|
self.__left = left
|
|
self.__right = right
|
|
|
|
def __str__(self):
|
|
return '(' + str(self.__left) + ' * ' + str(self.__right) + ')'
|
|
|
|
def __repr__(self):
|
|
return 'Product(' + repr(self.__left) + ', ' + repr(self.__right) + ')'
|
|
|
|
def Left(self):
|
|
return self.__left
|
|
|
|
def Right(self):
|
|
return self.__right
|
|
|
|
|
|
class BoundIntegerExpression(object):
|
|
"""Represents a constraint: IntegerExpression in domain."""
|
|
|
|
def __init__(self, expr, bounds):
|
|
self.__expr = expr
|
|
self.__bounds = bounds
|
|
|
|
def __str__(self):
|
|
if len(self.__bounds) == 2:
|
|
lb = self.__bounds[0]
|
|
ub = self.__bounds[1]
|
|
if lb > INT_MIN and ub < INT_MAX:
|
|
if lb == ub:
|
|
return str(self.__expr) + ' == ' + str(lb)
|
|
else:
|
|
return str(lb) + ' <= ' + str(self.__expr) + ' <= ' + str(ub)
|
|
elif lb > INT_MIN:
|
|
return str(self.__expr) + ' >= ' + str(lb)
|
|
elif ub < INT_MAX:
|
|
return str(self.__expr) + ' <= ' + str(ub)
|
|
else:
|
|
return 'True (unbounded expr ' + str(self.__expr) + ')'
|
|
else:
|
|
return str(self.__expr) + ' in [' + DisplayBounds(self.__bounds) + ']'
|
|
|
|
def Expression(self):
|
|
return self.__expr
|
|
|
|
def Bounds(self):
|
|
return self.__bounds
|
|
|
|
|
|
class Constraint(object):
|
|
"""Base class for constraints."""
|
|
|
|
def __init__(self, constraints):
|
|
self.__index = len(constraints)
|
|
self.__constraint = constraints.add()
|
|
|
|
def OnlyEnforceIf(self, boolvar):
|
|
self.__constraint.enforcement_literal.append(boolvar.Index())
|
|
|
|
def Index(self):
|
|
return self.__index
|
|
|
|
def ConstraintProto(self):
|
|
return self.__constraint
|
|
|
|
|
|
class IntervalVar(object):
|
|
"""Represents a Interval variable."""
|
|
|
|
def __init__(self, model, start_index, size_index, end_index,
|
|
is_present_index, name):
|
|
self.__model = model
|
|
self.__index = len(model.constraints)
|
|
self.__ct = self.__model.constraints.add()
|
|
self.__ct.interval.start = start_index
|
|
self.__ct.interval.size = size_index
|
|
self.__ct.interval.end = end_index
|
|
if is_present_index is not None:
|
|
self.__ct.enforcement_literal.append(is_present_index)
|
|
if name:
|
|
self.__ct.name = name
|
|
|
|
def Index(self):
|
|
return self.__index
|
|
|
|
def __str__(self):
|
|
return self.__ct.name
|
|
|
|
def __repr__(self):
|
|
interval = self.__ct.interval
|
|
if self.__ct.enforcement_literal:
|
|
return '%s(start = %s, size = %s, end = %s, is_present = %s)' % (
|
|
self.__ct.name, ShortName(self.__model, interval.start),
|
|
ShortName(self.__model, interval.size),
|
|
ShortName(self.__model, interval.end),
|
|
ShortName(self.__model, self.__ct.enforcement_literal[0]))
|
|
else:
|
|
return '%s(start = %s, size = %s, end = %s)' % (
|
|
self.__ct.name, ShortName(self.__model, interval.start),
|
|
ShortName(self.__model, interval.size),
|
|
ShortName(self.__model, interval.end))
|
|
|
|
|
|
class CpModel(object):
|
|
"""Wrapper class around the cp_model proto."""
|
|
|
|
def __init__(self):
|
|
self.__model = cp_model_pb2.CpModelProto()
|
|
self.__constant_map = {}
|
|
self.__optional_constant_map = {}
|
|
|
|
# Integer variable.
|
|
|
|
def NewIntVar(self, lb, ub, name):
|
|
"""Creates an integer variable with domain [lb, ub]."""
|
|
return IntVar(self.__model, [lb, ub], name)
|
|
|
|
def NewEnumeratedIntVar(self, bounds, name):
|
|
"""Creates an integer variable with an enumerated domain.
|
|
|
|
Args:
|
|
bounds: A flattened list of disjoint intervals.
|
|
name: The name of the variable.
|
|
|
|
Returns:
|
|
a variable whose domain is union[bounds[2*i]..bounds[2*i + 1]].
|
|
|
|
To create a variable with domain [1, 2, 3, 5, 7, 8], pass in the
|
|
array [1, 3, 5, 5, 7, 8].
|
|
"""
|
|
return IntVar(self.__model, bounds, name)
|
|
|
|
def NewOptionalIntVar(self, lb, ub, is_present, name):
|
|
"""Creates an optional integer variable."""
|
|
is_present_index = self.GetOrMakeBooleanIndex(is_present)
|
|
return IntVar(self.__model, [lb, ub], name, is_present_index)
|
|
|
|
def NewOptionalEnumeratedIntVar(self, bounds, is_present, name):
|
|
"""Creates an optional enumerated integer variable."""
|
|
is_present_index = self.GetOrMakeBooleanIndex(is_present)
|
|
return IntVar(self.__model, bounds, name, is_present_index)
|
|
|
|
def NewBoolVar(self, name):
|
|
"""Creates a 0-1 variable with the given name."""
|
|
return IntVar(self.__model, [0, 1], name)
|
|
|
|
# Integer constraints.
|
|
|
|
def AddLinearConstraint(self, terms, lb, ub):
|
|
"""Adds the constraints lb <= sum(terms) <= ub, where term = (var, coef)."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
for t in terms:
|
|
if not isinstance(t[0], IntVar):
|
|
raise TypeError('Wrong argument' + str(t))
|
|
AssertIsInt64(t[1])
|
|
model_ct.linear.vars.append(t[0].Index())
|
|
model_ct.linear.coeffs.append(t[1])
|
|
model_ct.linear.domain.extend([lb, ub])
|
|
return ct
|
|
|
|
def AddSumConstraint(self, variables, lb, ub):
|
|
"""Adds the constraints lb <= sum(variables) <= ub."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
for v in variables:
|
|
model_ct.linear.vars.append(v.Index())
|
|
model_ct.linear.coeffs.append(1)
|
|
model_ct.linear.domain.extend([lb, ub])
|
|
return ct
|
|
|
|
def AddLinearConstraintWithBounds(self, terms, bounds):
|
|
"""Adds the constraints sum(terms) in bounds, where term = (var, coef)."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
for t in terms:
|
|
if not isinstance(t[0], IntVar):
|
|
raise TypeError('Wrong argument' + str(t))
|
|
AssertIsInt64(t[1])
|
|
model_ct.linear.vars.append(t[0].Index())
|
|
model_ct.linear.coeffs.append(t[1])
|
|
model_ct.linear.domain.extend(bounds)
|
|
return ct
|
|
|
|
def Add(self, ct):
|
|
"""Adds a BoundIntegerExpression to the model."""
|
|
if isinstance(ct, BoundIntegerExpression):
|
|
coeffs_map, constant = ct.Expression().GetVarValueMap()
|
|
bounds = [CapSub(x, constant) for x in ct.Bounds()]
|
|
return self.AddLinearConstraintWithBounds(iteritems(coeffs_map), bounds)
|
|
else:
|
|
raise TypeError('Not supported: CpModel.Add(' + str(ct) + ')')
|
|
|
|
def AddAllDifferent(self, variables):
|
|
"""Adds AllDifferent(variables)."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.all_diff.vars.extend([self.GetOrMakeIndex(x) for x in variables])
|
|
return ct
|
|
|
|
def AddElement(self, index, variables, target):
|
|
"""Adds the element constraint: variables[index] == target."""
|
|
|
|
if not variables:
|
|
raise ValueError('AddElement expects a non empty variables array')
|
|
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.element.index = self.GetOrMakeIndex(index)
|
|
model_ct.element.vars.extend([self.GetOrMakeIndex(x) for x in variables])
|
|
model_ct.element.target = self.GetOrMakeIndex(target)
|
|
return ct
|
|
|
|
def AddCircuit(self, nexts):
|
|
"""Adds Circuit(nexts)."""
|
|
if not nexts:
|
|
raise ValueError('AddCircuit expects a non empty nexts array')
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.circuit.nexts.extend([self.GetOrMakeIndex(x) for x in nexts])
|
|
return ct
|
|
|
|
def AddAllowedAssignments(self, variables, tuples):
|
|
"""Adds AllowedAssignments(variables, [tuples])."""
|
|
if not variables:
|
|
raise ValueError('AddAllowedAssignments expects a non empty variables '
|
|
'array')
|
|
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.table.vars.extend([self.GetOrMakeIndex(x) for x in variables])
|
|
arity = len(variables)
|
|
for t in tuples:
|
|
if len(t) != arity:
|
|
raise TypeError('Tuple ' + str(t) + ' has the wrong arity')
|
|
for v in t:
|
|
AssertIsInt64(v)
|
|
model_ct.table.values.extend(t)
|
|
|
|
def AddForbiddenAssignments(self, variables, tuples):
|
|
"""Adds AddForbiddenAssignments(variables, [tuples])."""
|
|
|
|
if not variables:
|
|
raise ValueError('AddForbiddenAssignments expects a non empty variables '
|
|
'array')
|
|
|
|
index = len(self.__model.constraints)
|
|
self.AddAllowedAssignments(variables, tuples)
|
|
self.__model.constraints[index].table.negated = True
|
|
|
|
def AddAutomata(self, transition_variables, starting_state, final_states,
|
|
transition_triples):
|
|
"""Adds an automata constraint.
|
|
|
|
Args:
|
|
transition_variables: A non empty list of variables whose values
|
|
correspond to the labels of the arcs traversed by the automata.
|
|
starting_state: The initial state of the automata.
|
|
final_states: a non empty list of admissible final states.
|
|
transition_triples: A list of transition for the automata, in the
|
|
following format (current_state, variable_value, next_state).
|
|
|
|
|
|
Raises:
|
|
ValueError: if transition_variables, final_states, or transition_triples
|
|
are empty.
|
|
"""
|
|
|
|
if not transition_variables:
|
|
raise ValueError('AddAutomata expects a non empty transition_variables '
|
|
'array')
|
|
if not final_states:
|
|
raise ValueError('AddAutomata expects some final states')
|
|
|
|
if not transition_triples:
|
|
raise ValueError('AddAutomata expects some transtion triples')
|
|
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.automata.vars.extend(
|
|
[self.GetOrMakeIndex(x) for x in transition_variables])
|
|
AssertIsInt64(starting_state)
|
|
model_ct.automata.starting_state = starting_state
|
|
for v in final_states:
|
|
AssertIsInt64(v)
|
|
model_ct.automata.final_states.append(v)
|
|
for t in transition_triples:
|
|
if len(t) != 3:
|
|
raise TypeError('Tuple ' + str(t) + ' has the wrong arity (!= 3)')
|
|
AssertIsInt64(t[0])
|
|
AssertIsInt64(t[1])
|
|
AssertIsInt64(t[2])
|
|
model_ct.automata.transition_head.append(t[0])
|
|
model_ct.automata.transition_label.append(t[0])
|
|
model_ct.automata.transition_tail.append(t[0])
|
|
|
|
def AddInverse(self, variables, inverse_variables):
|
|
"""Adds Inverse(variables, inverse_variables)."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.inverse.f_direct.extend(
|
|
[self.GetOrMakeIndex(x) for x in variables])
|
|
model_ct.inverse.f_inverse.extend(
|
|
[self.GetOrMakeIndex(x) for x in inverse_variables])
|
|
return ct
|
|
|
|
def AddReservoirConstraint(self, times, demands, min_level, max_level):
|
|
"""Adds Reservoir(times, demands, min_level, max_level)."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.reservoir.times.extend([self.GetOrMakeIndex(x) for x in times])
|
|
model_ct.reservoir.demands.extend(demands)
|
|
model_ct.reservoir.min_level = min_level
|
|
model_ct.reservoir.max_level = max_level
|
|
return ct
|
|
|
|
def AddMapDomain(self, var, bool_var_array, offset=0):
|
|
"""Creates var == i + offset <=> bool_var_array[i] == true for all i."""
|
|
|
|
for i, bool_var in enumerate(bool_var_array):
|
|
b_index = bool_var.Index()
|
|
var_index = var.Index()
|
|
model_ct = self.__model.constraints.add()
|
|
model_ct.linear.vars.append(var_index)
|
|
model_ct.linear.coeffs.append(1)
|
|
model_ct.linear.domain.extend([offset + i, offset + i])
|
|
model_ct.enforcement_literal.append(b_index)
|
|
|
|
model_ct = self.__model.constraints.add()
|
|
model_ct.linear.vars.append(var_index)
|
|
model_ct.linear.coeffs.append(1)
|
|
model_ct.enforcement_literal.append(-b_index - 1)
|
|
if offset + i - 1 >= INT_MIN:
|
|
model_ct.linear.domain.extend([INT_MIN, offset + i - 1])
|
|
if offset + i + 1 <= INT_MAX:
|
|
model_ct.linear.domain.extend([offset + i + 1, INT_MAX])
|
|
|
|
def AddImplication(self, a, b):
|
|
"""Adds a => b."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.bool_or.literals.append(self.GetOrMakeBooleanIndex(b))
|
|
model_ct.enforcement_literal.append(self.GetOrMakeBooleanIndex(a))
|
|
return ct
|
|
|
|
def AddBoolOr(self, literals):
|
|
"""Adds Or(literals) == true."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.bool_or.literals.extend(
|
|
[self.GetOrMakeBooleanIndex(x) for x in literals])
|
|
return ct
|
|
|
|
def AddBoolAnd(self, literals):
|
|
"""Adds And(literals) == true."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.bool_and.literals.extend(
|
|
[self.GetOrMakeBooleanIndex(x) for x in literals])
|
|
return ct
|
|
|
|
def AddBoolXOr(self, literals):
|
|
"""Adds XOr(literals) == true."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.bool_xor.literals.extend(
|
|
[self.GetOrMakeBooleanIndex(x) for x in literals])
|
|
return ct
|
|
|
|
def AddMinEquality(self, target, variables):
|
|
"""Adds target == Min(variables)."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.int_min.vars.extend([self.GetOrMakeIndex(x) for x in variables])
|
|
model_ct.int_min.target = target.Index()
|
|
return ct
|
|
|
|
def AddMaxEquality(self, target, args):
|
|
"""Adds target == Max(variables)."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.int_max.vars.extend([self.GetOrMakeIndex(x) for x in args])
|
|
model_ct.int_max.target = target.Index()
|
|
return ct
|
|
|
|
def AddDivisionEquality(self, target, num, denom):
|
|
"""Creates target == num // denom."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.int_div.vars.extend(
|
|
[self.GetOrMakeIndex(num),
|
|
self.GetOrMakeIndex(denom)])
|
|
model_ct.int_div.target = target.Index()
|
|
return ct
|
|
|
|
def AddModuloEquality(self, target, var, mod):
|
|
"""Creates target = var % mod."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.int_mod.vars.extend(
|
|
[self.GetOrMakeIndex(var),
|
|
self.GetOrMakeIndex(mod)])
|
|
model_ct.int_mod.target = target.Index()
|
|
return ct
|
|
|
|
def AddProdEquality(self, target, args):
|
|
"""Creates target == PROD(args)."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.int_prod.vars.extend([self.GetOrMakeIndex(x) for x in args])
|
|
model_ct.int_prod.target = target.Index()
|
|
return ct
|
|
|
|
# Scheduling support
|
|
|
|
def NewIntervalVar(self, start, size, end, name):
|
|
start_index = self.GetOrMakeIndex(start)
|
|
size_index = self.GetOrMakeIndex(size)
|
|
end_index = self.GetOrMakeIndex(end)
|
|
self.AssertIntVarIsNotOptional(start_index)
|
|
self.AssertIntVarIsNotOptional(size_index)
|
|
self.AssertIntVarIsNotOptional(end_index)
|
|
return IntervalVar(self.__model, start_index, size_index, end_index, None,
|
|
name)
|
|
|
|
def NewOptionalIntervalVar(self, start, size, end, is_present, name):
|
|
is_present_index = self.GetOrMakeBooleanIndex(is_present)
|
|
start_index = self.GetOrMakeOptionalIndex(start, is_present)
|
|
size_index = self.GetOrMakeIndex(size) # Currently, not optional.
|
|
end_index = self.GetOrMakeOptionalIndex(end, is_present)
|
|
self.CheckOptionalIntVarLiteral(start_index, is_present_index)
|
|
self.CheckOptionalIntVarLiteral(end_index, is_present_index)
|
|
return IntervalVar(self.__model, start_index, size_index, end_index,
|
|
is_present_index, name)
|
|
|
|
def AddNoOverlap(self, interval_vars):
|
|
"""Adds NoOverlap(interval_vars)."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.no_overlap.intervals.extend(
|
|
[self.GetIntervalIndex(x) for x in interval_vars])
|
|
return ct
|
|
|
|
def AddNoOverlap2D(self, x_intervals, y_intervals):
|
|
"""Adds NoOverlap2D(x_intervals, y_intervals)."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.no_overlap_2d.x_intervals.extend(
|
|
[self.GetIntervalIndex(x) for x in x_intervals])
|
|
model_ct.no_overlap_2d.y_intervals.extend(
|
|
[self.GetIntervalIndex(x) for x in y_intervals])
|
|
return ct
|
|
|
|
def AddCumulative(self, intervals, demands, capacity):
|
|
"""Adds Cumulative(intervals, demands, capacity)."""
|
|
ct = Constraint(self.__model.constraints)
|
|
model_ct = self.__model.constraints[ct.Index()]
|
|
model_ct.cumulative.intervals.extend(
|
|
[self.GetIntervalIndex(x) for x in intervals])
|
|
model_ct.cumulative.demands.extend(
|
|
[self.GetOrMakeIndex(x) for x in demands])
|
|
model_ct.cumulative.capacity = self.GetOrMakeIndex(capacity)
|
|
return ct
|
|
|
|
# Helpers.
|
|
|
|
def __str__(self):
|
|
return str(self.__model)
|
|
|
|
def ModelProto(self):
|
|
return self.__model
|
|
|
|
def Negated(self, index):
|
|
return -index - 1
|
|
|
|
def GetOrMakeIndex(self, arg):
|
|
if isinstance(arg, IntVar):
|
|
return arg.Index()
|
|
elif (isinstance(arg, _ProductCst) and
|
|
isinstance(arg.Expression(), IntVar) and arg.Coefficient() == -1):
|
|
return -arg.Expression().Index() - 1
|
|
elif isinstance(arg, numbers.Integral):
|
|
AssertIsInt64(arg)
|
|
return self.GetOrMakeIndexFromConstant(arg)
|
|
else:
|
|
raise TypeError('NotSupported: model.GetOrMakeIndex(' + str(arg) + ')')
|
|
|
|
def GetOrMakeOptionalIndex(self, arg, is_present):
|
|
if isinstance(arg, IntVar):
|
|
return arg.Index()
|
|
elif (isinstance(arg, _ProductCst) and
|
|
isinstance(arg.Expression(), IntVar) and arg.Coefficient() == -1):
|
|
return -arg.Expression().Index() - 1
|
|
elif isinstance(arg, numbers.Integral):
|
|
AssertIsInt64(arg)
|
|
return self.GetOrMakeOptionalIndexFromConstant(arg, is_present)
|
|
else:
|
|
raise TypeError('NotSupported: model.GetOrMakeOptionalIndex(%s, %s)' %
|
|
(arg, is_present))
|
|
|
|
def GetOrMakeBooleanIndex(self, arg):
|
|
if isinstance(arg, IntVar):
|
|
self.AssertIsBooleanVariable(arg)
|
|
return arg.Index()
|
|
elif isinstance(arg, _NotBooleanVariable):
|
|
self.AssertIsBooleanVariable(arg.Not())
|
|
return arg.Index()
|
|
elif isinstance(arg, numbers.Integral):
|
|
AssertIsBoolean(arg)
|
|
return self.GetOrMakeIndexFromConstant(arg)
|
|
else:
|
|
raise TypeError('NotSupported: model.GetOrMakeBooleanIndex(' + str(arg) +
|
|
')')
|
|
|
|
def GetIntervalIndex(self, arg):
|
|
if not isinstance(arg, IntervalVar):
|
|
raise TypeError('NotSupported: model.GetIntervalIndex(%s)' % arg)
|
|
return arg.Index()
|
|
|
|
def GetOrMakeIndexFromConstant(self, value):
|
|
if value in self.__constant_map:
|
|
return self.__constant_map[value]
|
|
index = len(self.__model.variables)
|
|
var = self.__model.variables.add()
|
|
var.domain.extend([value, value])
|
|
self.__constant_map[value] = index
|
|
return index
|
|
|
|
def GetOrMakeOptionalIndexFromConstant(self, value, is_present):
|
|
self.AssertIsBooleanVariable(is_present)
|
|
if (is_present, value) in self.__optional_constant_map:
|
|
return self.__optional_constant_map[(is_present, value)]
|
|
index = len(self.__model.variables)
|
|
var = self.__model.variables.add()
|
|
var.domain.extend([value, value])
|
|
var.enforcement_literal.append(self.GetOrMakeBooleanIndex(is_present))
|
|
self.__optional_constant_map[(is_present, value)] = index
|
|
return index
|
|
|
|
def VarIndexToVarProto(self, var_index):
|
|
if var_index > 0:
|
|
return self.__model.variables[var_index]
|
|
else:
|
|
return self.__model.variables[-var_index - 1]
|
|
|
|
def AssertIntVarIsNotOptional(self, var_index):
|
|
var = self.VarIndexToVarProto(var_index)
|
|
if var.enforcement_literal:
|
|
raise TypeError('Variable %s should not be marked as optional' %
|
|
ShortName(self.__model, var_index))
|
|
|
|
def CheckOptionalIntVarLiteral(self, var_index, lit_index):
|
|
var = self.VarIndexToVarProto(var_index)
|
|
if (len(var.enforcement_literal) != 1 or
|
|
var.enforcement_literal[0] != lit_index):
|
|
raise TypeError('Variable %s should be marked optional with literal %s' %
|
|
(ShortName(self.__model, var_index),
|
|
ShortName(self.__model, lit_index)))
|
|
|
|
def _SetObjective(self, obj, minimize):
|
|
"""Sets the objective of the model."""
|
|
if isinstance(obj, IntVar):
|
|
self.__model.ClearField('objective')
|
|
self.__model.objective.coeffs.append(1)
|
|
self.__model.objective.offset = 0
|
|
if minimize:
|
|
self.__model.objective.vars.append(obj.Index())
|
|
self.__model.objective.scaling_factor = 1
|
|
else:
|
|
self.__model.objective.vars.append(self.Negated(obj.Index()))
|
|
self.__model.objective.scaling_factor = -1
|
|
elif isinstance(obj, IntegerExpression):
|
|
coeffs_map, constant = obj.GetVarValueMap()
|
|
self.__model.ClearField('objective')
|
|
if minimize:
|
|
self.__model.objective.scaling_factor = 1
|
|
self.__model.objective.offset = constant
|
|
else:
|
|
self.__model.objective.scaling_factor = -1
|
|
self.__model.objective.offset = -constant
|
|
for v, c, in iteritems(coeffs_map):
|
|
self.__model.objective.coeffs.append(c)
|
|
if minimize:
|
|
self.__model.objective.vars.append(v.Index())
|
|
else:
|
|
self.__model.objective.vars.append(self.Negated(v.Index()))
|
|
else:
|
|
raise TypeError('TypeError: ' + str(obj) + ' is not a valid objective')
|
|
|
|
def Minimize(self, obj):
|
|
"""Sets the objective of the model to minimize(obj)."""
|
|
self._SetObjective(obj, minimize=True)
|
|
|
|
def Maximize(self, obj):
|
|
"""Sets the objective of the model to maximize(obj)."""
|
|
self._SetObjective(obj, minimize=False)
|
|
|
|
def HasObjective(self):
|
|
return self.__model.HasField('objective')
|
|
|
|
def AssertIsBooleanVariable(self, x):
|
|
if isinstance(x, IntVar):
|
|
var = self.__model.variables[x.Index()]
|
|
if len(var.domain) != 2 or var.domain[0] < 0 or var.domain[1] > 1:
|
|
raise TypeError('TypeError: ' + str(x) + ' is not a boolean variable')
|
|
elif not isinstance(x, _NotBooleanVariable):
|
|
raise TypeError('TypeError: ' + str(x) + ' is not a boolean variable')
|
|
|
|
|
|
def EvaluateExpression(expression, solution):
|
|
"""Evaluate an integer expression against a solution."""
|
|
value = 0
|
|
to_process = [(expression, 1)]
|
|
while to_process:
|
|
expr, coef = to_process.pop()
|
|
if isinstance(expr, _ProductCst):
|
|
to_process.append((expr.Expression(), coef * expr.Coefficient()))
|
|
elif isinstance(expr, _SumArray):
|
|
for e in expr.Array():
|
|
to_process.append((e, coef))
|
|
value += expr.Constant() * coef
|
|
elif isinstance(expr, IntVar):
|
|
value += coef * solution.solution[expr.Index()]
|
|
elif isinstance(expr, _NotBooleanVariable):
|
|
raise TypeError('Cannot interpret literals in a integer expression.')
|
|
return value
|
|
|
|
|
|
class CpSolverSolutionCallback(pywrapsat.PySolutionCallback):
|
|
"""Nicer solution callback that uses the CpSolver class."""
|
|
|
|
def __init__(self):
|
|
self.__current_solution = None
|
|
|
|
def Wrap(self, solution_proto):
|
|
self.__current_solution = solution_proto
|
|
self.NewSolution()
|
|
|
|
def Value(self, expression):
|
|
"""Returns the value of an integer expression."""
|
|
if not self.__current_solution:
|
|
raise RuntimeError('Solve() has not be called.')
|
|
return EvaluateExpression(expression, self.__current_solution)
|
|
|
|
def ObjectiveValue(self):
|
|
"""Returns the value of the objective."""
|
|
return self.__current_solution.objective_value
|
|
|
|
def NewSolution(self):
|
|
pass
|
|
|
|
|
|
class CpSolver(object):
|
|
"""Main solver class."""
|
|
|
|
def __init__(self):
|
|
self.__model = None
|
|
self.__solution = None
|
|
self.parameters = sat_parameters_pb2.SatParameters()
|
|
|
|
def Solve(self, model):
|
|
"""Solves the given model and returns the solve status."""
|
|
self.__solution = pywrapsat.SatHelper.SolveWithParameters(
|
|
model.ModelProto(), self.parameters)
|
|
return self.__solution.status
|
|
|
|
def SolveWithSolutionObserver(self, model, callback):
|
|
"""Solves a problem and pass each solution found to the callback."""
|
|
parameters = sat_parameters_pb2.SatParameters()
|
|
self.__solution = (
|
|
pywrapsat.SatHelper.SolveWithParametersAndSolutionObserver(
|
|
model.ModelProto(), parameters, callback))
|
|
return self.__solution.status
|
|
|
|
def SearchForAllSolutions(self, model, callback):
|
|
"""Search for all solutions of a satisfiability problem."""
|
|
if model.HasObjective():
|
|
raise TypeError('Search for all solutions is only defined on '
|
|
'satisfiability problems')
|
|
# Store old values.
|
|
enumerate_all = self.parameters.enumerate_all_solutions
|
|
presolve = self.parameters.cp_model_presolve
|
|
self.parameters.enumerate_all_solutions = True
|
|
self.parameters.cp_model_presolve = False
|
|
self.__solution = (
|
|
pywrapsat.SatHelper.SolveWithParametersAndSolutionObserver(
|
|
model.ModelProto(), self.parameters, callback))
|
|
# Restore parameters.
|
|
self.parameters.enumerate_all_solutions = enumerate_all
|
|
self.parameters.cp_model_presolve = presolve
|
|
return self.__solution.status
|
|
|
|
def Value(self, expression):
|
|
"""Returns the value of an integer expression."""
|
|
if not self.__solution:
|
|
raise RuntimeError('Solve() has not be called.')
|
|
return EvaluateExpression(expression, self.__solution)
|
|
|
|
def ObjectiveValue(self):
|
|
"""Returns the objective value found after solve."""
|
|
return self.__solution.objective_value
|
|
|
|
def StatusName(self, status):
|
|
return cp_model_pb2.CpSolverStatus.Name(status)
|
|
|
|
def NumBooleans(self):
|
|
return self.__solution.num_booleans
|
|
|
|
def NumConflicts(self):
|
|
return self.__solution.num_conflicts
|
|
|
|
def NumBranches(self):
|
|
return self.__solution.num_branches
|
|
|
|
def WallTime(self):
|
|
return self.__solution.wall_time
|