public Sudoku() :
base(1, 9)
{
int size = 9;
m_Matrix = new IntVarMatrix(m_Solver, size, size, new IntInterval(1, size));
for (int idx = 0; idx < size; ++idx)
{
IntVarListAllDifferent constraint = m_Matrix.Row(idx).AllDifferent();
constraint.Level = PropagateLevel.High;
m_Solver.Add(constraint);
}
for (int idx = 0; idx < size; ++idx)
{
IntVarListAllDifferent constraint = m_Matrix.Col(idx).AllDifferent();
constraint.Level = PropagateLevel.High;
m_Solver.Add(constraint);
}
for (int subRow = 0; subRow < 3; ++subRow)
{
for (int subCol = 0; subCol < 3; ++subCol)
{
IntVarListAllDifferent constraint = m_Matrix.Matrix(subRow * 3, subCol * 3, 3, 3).VarList.AllDifferent();
constraint.Level = PropagateLevel.High;
m_Solver.Add(constraint);
}
}
}
public MagicSquare(int count) :
base(0, 10000)
{
m_Matrix = new IntVarMatrix(m_Solver, count, count, new IntInterval(1, count * count));
m_MagicConstant = (count * (count * count + 1)) / 2;
for (int idx = 0; idx < count; ++idx)
{
m_Solver.Add(m_Matrix.Row(idx).Sum(m_MagicConstant));
m_Solver.Add(m_Matrix.Col(idx).Sum(m_MagicConstant));
}
bool mm = false;
if (mm)
{
if (count % 2 == 0)
{
int magicConstant2 = m_MagicConstant / 2;
for (int idx = 0; idx < count; ++idx)
{
m_Solver.Add(m_Matrix.Row(idx, 0, count / 2).Sum(magicConstant2));
m_Solver.Add(m_Matrix.Row(idx, count / 2, count).Sum(magicConstant2));
m_Solver.Add(m_Matrix.Col(idx, 0, count / 2).Sum(magicConstant2));
m_Solver.Add(m_Matrix.Col(idx, count / 2, count).Sum(magicConstant2));
}
}
}
m_Solver.Add(m_Matrix.DiagLeftTopToBottomRight().Sum(m_MagicConstant));
m_Solver.Add(m_Matrix.DiagRightTopToBottomLeft().Sum(m_MagicConstant));
IntVarListAllDifferent ad = m_Matrix.VarList.AllDifferent();
m_Solver.Add(ad);
// remove symmetry
m_Solver.Add(m_Matrix[0, 0] < m_Matrix[0, count - 1]);
m_Solver.Add(m_Matrix[0, 0] < m_Matrix[count - 1, count - 1]);
m_Solver.Add(m_Matrix[0, 0] < m_Matrix[count - 1, 0]);
m_Solver.Add(m_Matrix[0, count - 1] < m_Matrix[count - 1, 0]);
}
