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]);
}
