10.9.4.4 Arithmetic-Logical Constraints

smt(:ConstraintBody)   since release 4.2,deprecated

The arithmetic, membership, and propositional constraints described earlier are transformed at compile time to conjunctions of library constraints. Although linear in the size of the source code, the expansion of a propositional formula over reifiable constraints to library goals can have time and memory overheads, and propagates disjunctions very weakly. Temporary variables holding intermediate values may have to be introduced, and the grain size of the constraint solver invocations can be rather small. As an alternative to the default propagation of such constraint formulas, this constraint is a front-end to the case/[3,4] propagator, which treats such a formula globally.

Although often convenient, this constraint is deprecated, because it cannot guarantee better performance than a decomposition, nor any particular level of consistency.

ConstraintBody should be of one of the following forms, or a propositional combination of such forms. See Syntax of Indexicals for the exact definition:

• var in ConstantRange
• element(var,CList,var)
• table([VList],CTable)
• LinExpr RelOp LinExpr
• var { X stands for X#=1 }

count(+Val,+List,+RelOp,?Count)   since release 4.0.5

where Val is an integer, List is a list of domain variables, Count a domain variable, and RelOp is a relational symbol as in Arithmetic Constraints. True if N is the number of occurrences of Val in List and N RelOp Count.

Corresponds to count_* and exactly in MiniZinc.

If you want to count the occurrences of multiple values, the following constraint is a better alternative.

global_cardinality(+Xs,+Vals)

global_cardinality(+Xs,+Vals,+Options)

where Xs = [X1,…,Xd] is a list of domain variables, and Vals = [K1-V1,…,Kn-Vn] is a list of pairs where each key Ki is a unique integer and Vi is a domain variable. True if every element of Xs is equal to some key and for each pair Ki-Vi, exactly Vi elements of Xs are equal to Ki.

If either Xs or Vals is ground, and in many other special cases, then global_cardinality/[2,3] maintains domain consistency, but generally, bounds consistency cannot be guaranteed. A domain-consistency algorithm [Regin 96] is used, roughly linear in the total size of the domains.

Corresponds to global_cardinality* and distribute in MiniZinc.

Options is a list of zero or more of the following:

consistency(Cons)

Which filtering algorithm to use. One of the following:

domain

The constraint will use the algorithm mentioned above. The default.

bounds

The constraint will use the algorithm mentioned above.

value

The constraint will use a simple algorithm, which prevents too few or too many of the Xs from taking values among the Vals.

cost(Cost,Matrix)

Cost is constrained to equal a cost computed as a function of Xs and Matrix, which should be a d*n matrix of integers, represented as a list of d lists, each of length n. Assume that each Xi equals K(pi). The cost is Matrix[1,p1]+…+Matrix[d,pd].

With this option, a domain-consistency algorithm [Regin 99] is used, the complexity of which is roughly O(d(m + n log n)) where m is the total size of the domains.

all_equal(+Variables)   reifiable   since release 4.7

where Variables is a list of domain variables. They are constrained to take the same value.

Corresponds to all_equal in MiniZinc.

all_equal_reif(+Variables, ?Reif))   since release 4.7

is the reified version of all_equal/1, i.e., Reif is 1 if all_equal/1 with the same argument holds; otherwise, Reif is 0.

all_different(+Variables)

all_different(+Variables, +Options)

all_distinct(+Variables)

all_distinct(+Variables, +Options)

where Variables is a list of domain variables. Each variable is constrained to take a value that is unique among the variables. Declaratively, this is equivalent to an inequality constraint for each pair of variables.

Corresponds to all_different in MiniZinc.

Options is a list of zero or more of the following:

L #= R   since release 4.3

where R should be an integer, and L an expressions of one of the following forms, where X1, …, Xj occur among Variables:

X1 + ... + Xj

X1*X1 + ... + Xj*Xj

X1 * ... * Xj

The given equation is a side constraint for the constraint to hold. A special bounds-consistency algorithm is used, which considers the main constraint and the side constraints globally. This option is only valid if the domains of X1, …, Xj consist of integers strictly greater than zero.

consistency(Cons)

Which algorithm to use, one of the following:

domain

The default for all_distinct/[1,2] and assignment/[2,3]. A domain-consistency algorithm [Regin 94] is used, roughly linear in the total size of the domains.

bounds

A bounds-consistency algorithm [Lopez-Ortiz 03] is used, which runs in O(n log n) time for n variables.

value

The default for all_different/[1,2]. An algorithm achieving exactly the same pruning as a set of pairwise inequality constraints is used, roughly linear in the number of variables.

circuit(Boolean)

If true, then circuit(Xs) must hold for the constraint to be true.

subcircuit(Boolean)   since release 4.6

If true, then subcircuit(Xs) must hold for the constraint to be true.

cost(Cost,Matrix)   since release 4.9.0

Cost is constrained to equal a cost computed as a function of Xs and Matrix, namely Matrix[1,X1]+…+Matrix[n,Xn]. Matrix should be an n*m matrix of integers, represented as a list of lists.

all_different_except_0(+Variables)   since release 4.6

all_distinct_except_0(+Variables)   since release 4.6

where Variables is a list of domain variables with finite bounds. The variables are constrained to be all different, except those elements that are assigned the value 0.

In terms of consistency, all_different_except_0/1 corresponds to all_different/1, and all_distinct_except_0/1 to all_distinct/1.

Corresponds to alldifferent_except_0 in MiniZinc.

nvalue(?N, +Variables)

where Variables is a list of domain variables with finite bounds, and N is a domain variable. True if N is the number of distinct values taken by Variables. Approximates bounds consistency in N and domain consistency in Variables. Can be thought of as a relaxed version of all_distinct/2.

Corresponds to nvalue in MiniZinc.

The following is a constraint over two lists of length n of variables. Each variable is constrained to take a value in 1..n that is unique for its list. Furthermore, the lists are dual in a sense described below.

assignment(+Xs, +Ys)

assignment(+Xs, +Ys, +Options)

where Xs = [X1,…,Xn] and Ys = [Y1,…,Yn] are lists of domain variables. True if all Xi, Yi are in 1..n and Xi=j iff Yj=i.

Corresponds to inverse in MiniZinc.

Options is a list of zero or more of the following, where Boolean must be true or false (false is the default):

consistency(Cons)

Same meaning as for all_different/2.

circuit(Boolean)

If true, then circuit(Xs,Ys) must hold for the constraint to be true.

subcircuit(Boolean)   since release 4.6

If true, then subcircuit(Xs,Ys) must hold for the constraint to be true.

cost(Cost,Matrix)

Cost is constrained to equal a cost computed as a function of Xs and Matrix, namely Matrix[1,X1]+…+Matrix[n,Xn]. Matrix should be an n*n matrix of integers, represented as a list of lists.

The following is a constraint over a list of length n of variables, which all are constrained to 1..n.

symmetric_all_different(+Xs)   since release 4.6

symmetric_all_distinct(+Xs)   since release 4.6

True if the list is all different, and that all Xi in 1..n, and Xi=j iff Xj=i.

Corresponds to symmetric_all_different in MiniZinc. symmetric_all_distinct/1 maintains stronger consistency.

The following constraint captures the relation between a list of values, a list of the values in ascending order, and their positions in the original list:

sorting(+Xs,+Ps,+Ys)

where Xs = [X1,…,Xn], Ps = [P1,…,Pn], and Ys = [Y1,…,Yn] are lists of domain variables. The constraint holds if the following are true:

• Ys is in ascending order.
• Ps is a permutation of 1..n.
• for all i in 1..n : Xi = Y(Pi)

In practice, the underlying algorithm [Mehlhorn 00] is likely to achieve bounds consistency, and is guaranteed to do so if Ps is ground or completely free.

Corresponds to sort in MiniZinc.

The following constraint is a generalization of sorting/3, namely:

• It sorts not domain variables, but tuples of them.
• The tuples are split into a key prefix and a value suffix. They are sorted wrt. the key part only.
• The sort is stable: if two tuples have identical keys, then their relative order is preserved in the output.

keysorting(+Xs,+Ys)   since release 4.3.1

keysorting(+Xs,+Ys,+Options)

where Xs and Ys are lists of tuples of domain variables. Both lists should be of the same length n, and all tuples should have the same length m. The constraint holds if the following are true:

• Ps is a permutation of 1..n.
• for all i in 1..n, j in 1..m : Ys[i,j] = Xs[Ps[i],j].
• [[Ys[i,1],...,Ys[i,k],Ps[i]] | i in 1..n] is in lex ascending order, where k equals Keys as defined below in the options.

The filtering algorithm is based on [Mehlhorn 00] and endeavors to achieve bounds consistency, but does not guarantee it.

Corresponds to Prolog’s keysort/2. In particular, the sort is stable by definition. Corresponds to arg_sort in MiniZinc.

Options is a list of zero or more of the following:

keys(Keys)

where Keys should be a positive integer, denoting the length of the key prefix. The default is 1.

permutation(Ps)

where Ps should be a list of length n of domain variables. Its meaning is described above.

The following constraints express the fact that several vectors of domain variables are in ascending lexicographic order:

lex_chain(+Vectors)

lex_chain(+Vectors,+Options)

where Vectors is a list of vectors (lists) of domain variables with finite bounds. The constraint holds if Vectors are in ascending lexicographic order.

Corresponds to *lex2, lex_greater*, lex_less* in MiniZinc.

Options is a list of zero or more of the following:

op(Op)

If Op is the atom #=< (the default), then the constraints holds if Vectors are in non-descending lexicographic order. If Op is the atom #<, then the constraints holds if Vectors are in strictly ascending lexicographic order.

increasing

This option imposes the additional constraint that each vector in Vectors be sorted in strictly ascending order.

among(Least,Most,Values)

If given, then Least and Most should be integers such that 0 <= Least <= Most and Values should be a list of distinct integers. This option imposes the additional constraint on each vector in Vectors that at least Least and at most Most elements belong to Values.

global(Boolean)   since release 4.2.1

if true, then a more expensive algorithm [Carlsson & Beldiceanu 02], which guaranteed domain consistency unless the increasing/0 or among/3 options are given, will be used.

In the following constraints, a literal is either a term X or a term #\ X, where X is a 0/1 variable. They maintain domain consistency:

bool_and(+Lits, +Lit)   since release 4.3

where Lits is a list of literals [L0,...,Lj] and Lit is a literal. True if Lit equals the Boolean conjunction of Lits, and usually more efficient than the equivalent L0#/\...#/\Lj #<=> Lit.

bool_or(+Lits, +Lit)   since release 4.3

where Lits is a list of literals [L0,...,Lj] and Lit is a literal. True if Lit equals the Boolean disjunction of Lits, and usually more efficient than the equivalent L0#\/...#\/Lj #<=> Lit.

bool_xor(+Lits, +Lit)   since release 4.3

where Lits is a list of literals [L0,...,Lj] and Lit is a literal. True if Lit equals the Boolean exclusive or of Lits, and usually more efficient than the equivalent L0#\...#\Lj #<=> Lit.

bool_channel(+Lits, ?Y, +Relop, +Offset)   since release 4.3

where Lits is a list of literals [L0,...,Lj], Y is a domain variable, RelOp is an arithmetic comparison as in Syntax of Arithmetic Expressions, and Offset is an integer. Expresses the constraint Li #<=> (Y RelOp i+Offset) for i in 0..j. Usually more efficient than a bunch of reified comparisons between a given variable and a sequence of integers.