org.chocosolver.graphsolver.cstrs.symmbreaking.PropIncrementalAdjacencyUndirectedMatrix Maven / Gradle / Ivy

Go to download
Show more of this group Show more artifacts with this name
Show all versions of choco-graph Show documentation
A graph module for Choco Solver.
There is a newer version: 4.2.2
package org.chocosolver.graphsolver.cstrs.symmbreaking;

import org.chocosolver.graphsolver.variables.GraphEventType;
import org.chocosolver.graphsolver.variables.IUndirectedGraphVar;
import org.chocosolver.graphsolver.variables.delta.IGraphDeltaMonitor;
import org.chocosolver.solver.constraints.Propagator;
import org.chocosolver.solver.constraints.PropagatorPriority;
import org.chocosolver.solver.exception.ContradictionException;
import org.chocosolver.solver.variables.BoolVar;
import org.chocosolver.solver.variables.Variable;
import org.chocosolver.solver.variables.events.IntEventType;
import org.chocosolver.util.ESat;
import org.chocosolver.util.objects.setDataStructures.ISet;
import org.chocosolver.util.procedure.PairProcedure;
import org.chocosolver.util.tools.ArrayUtils;

import java.util.HashSet;
import java.util.Set;

/**
 * @author Моклев Вячеслав
 */
public class PropIncrementalAdjacencyUndirectedMatrix extends Propagator {

	private IUndirectedGraphVar graph;
	private IGraphDeltaMonitor gdm;
	private PairProcedure enforce;
	private PairProcedure remove;
	private int n;
	private BoolVar[] t;

	public PropIncrementalAdjacencyUndirectedMatrix(IUndirectedGraphVar graphVar, BoolVar[] t) {
		super(ArrayUtils.append(new Variable[]{graphVar}, t), PropagatorPriority.LINEAR, true);
		graph = graphVar;
		gdm = graph.monitorDelta(this);
		enforce = (PairProcedure) (from, to) -> {
			t[from + to * n].instantiateTo(1, this);
			t[to + from * n].instantiateTo(1, this);
		};
		remove = (PairProcedure) (from, to) -> {
			t[from + to * n].instantiateTo(0, this);
			t[to + from * n].instantiateTo(0, this);
		};
		n = graphVar.getNbMaxNodes();
		this.t = t;
	}

	@Override
	public void propagate(int evtmask) throws ContradictionException {
		// first, nonincremental propagation
		for (int i = 0; i < n; i++) {
			t[i + i * n].instantiateTo(0, this);
		}
		propagateGraphChanged();
		propagateTChanged();
		// initializing incremental data-structures
		gdm.unfreeze();
	}

	private void propagateGraphChanged() throws ContradictionException {
		for (int i = 0; i < n; i++) {
			for (int j = 0; j < n; j++) {
				if (t[i + j * n].isInstantiatedTo(1)) {
					graph.enforceArc(i, j, this);
				}
				if (t[i + j * n].isInstantiatedTo(0)) {
					graph.removeArc(i, j, this);
				}
			}
		}
	}

	private void propagateTChanged() throws ContradictionException {
		for (int u = 0; u < n; u++) {
			for (int v: graph.getMandNeighOf(u)) {
				t[u + v * n].instantiateTo(1, this);
			}
		}
		for (int u = 0; u < n; u++) {
			Set set = new HashSet<>();
			for (int v: graph.getPotNeighOf(u)) {
				set.add(v);
			}
			for (int v = 0; v < n; v++) {
				if (!set.contains(v)) {
					t[u + v * n].instantiateTo(0, this);
				}
			}
		}
	}

	@Override
	public void propagate(int idxVarInProp, int mask) throws ContradictionException {
		// incremental propagation
		if (idxVarInProp == 0) {
			gdm.freeze();
			gdm.forEachArc(enforce, GraphEventType.ADD_ARC);
			gdm.forEachArc(remove, GraphEventType.REMOVE_ARC);
			gdm.unfreeze();
		} else {
			propagateTChanged();
		}
	}

	@Override
	public int getPropagationConditions(int vIdx) {
		if (vIdx == 0) {
			return GraphEventType.ADD_ARC.getMask() | GraphEventType.REMOVE_ARC.getMask();
		} else {
			return IntEventType.boundAndInst();
		}
	}

	@Override
	public ESat isEntailed() {
		for (int i = 0; i < n; i++) {
			ISet children = graph.getMandNeighOf(i);
			for (int j = 0; j < n; j++) {
				if ((t[i + j * n].isInstantiatedTo(0) || t[j + i * n].isInstantiatedTo(0)) && children.contains(j)) {
					return ESat.FALSE;
				}
			}
		}
		for (int i = 0; i < n; i++) {
			ISet children = graph.getPotNeighOf(i);
			for (int j = 0; j < n; j++) {
				if ((t[i + j * n].isInstantiatedTo(1) || t[j + i * n].isInstantiatedTo(1)) && !children.contains(j)) {
					return ESat.FALSE;
				}
			}
		}
		if (graph.isInstantiated()) {
			return ESat.TRUE;
		}
		return ESat.UNDEFINED;
	}
}