gecoder-with-gecode 0.9.0 → 0.9.1

data/ext/gecode-2.2.0/examples/magic-square.cc

@@ -0,0 +1,138 @@
+/* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
+/*
+ *  Main authors:
+ *     Christian Schulte <schulte@gecode.org>
+ *
+ *  Copyright:
+ *     Christian Schulte, 2001
+ *
+ *  Last modified:
+ *     $Date: 2007-11-30 13:58:34 +0100 (Fri, 30 Nov 2007) $ by $Author: tack $
+ *     $Revision: 5524 $
+ *
+ *  This file is part of Gecode, the generic constraint
+ *  development environment:
+ *     http://www.gecode.org
+ *
+ *  Permission is hereby granted, free of charge, to any person obtaining
+ *  a copy of this software and associated documentation files (the
+ *  "Software"), to deal in the Software without restriction, including
+ *  without limitation the rights to use, copy, modify, merge, publish,
+ *  distribute, sublicense, and/or sell copies of the Software, and to
+ *  permit persons to whom the Software is furnished to do so, subject to
+ *  the following conditions:
+ *
+ *  The above copyright notice and this permission notice shall be
+ *  included in all copies or substantial portions of the Software.
+ *
+ *  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+ *  EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+ *  MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+ *  NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+ *  LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+ *  OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+ *  WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+ *
+ */
+
+#include "examples/support.hh"
+#include "gecode/minimodel.hh"
+
+/**
+ * \brief %Example: Magic squares
+ *
+ * Compute magic squares of arbitrary size
+ *
+ * See problem 19 at http://www.csplib.org/.
+ *
+ * \ingroup ExProblem
+ *
+ */
+class MagicSquare : public Example {
+private:
+  /// Size of magic square
+  const int n;
+  /// Fields of square
+  IntVarArray x;
+
+public:
+  /// Post constraints
+  MagicSquare(const SizeOptions& opt)
+    : n(opt.size()), x(this,n*n,1,n*n) {
+    // Number of fields on square
+    const int nn = n*n;
+
+    // Sum of all a row, column, or diagonal
+    const int s  = nn*(nn+1) / (2*n);
+
+    // Matrix-wrapper for the square
+    MiniModel::Matrix<IntVarArray> m(x, n, n);
+
+    for (int i = n; i--; ) {
+        linear(this, m.row(i), IRT_EQ, s, opt.icl());
+        linear(this, m.col(i), IRT_EQ, s, opt.icl());
+    }
+    // Both diagonals must have sum s
+    {
+      IntVarArgs d1y(n);
+      IntVarArgs d2y(n);
+      for (int i = n; i--; ) {
+        d1y[i] = m(i,i);
+        d2y[i] = m(n-i-1,i);
+      }
+      linear(this, d1y, IRT_EQ, s, opt.icl());
+      linear(this, d2y, IRT_EQ, s, opt.icl());
+    }
+
+    // All fields must be distinct
+    distinct(this, x, opt.icl());
+
+    // Break some (few) symmetries
+    rel(this, m(0,0), IRT_GR, m(0,n-1));
+    rel(this, m(0,0), IRT_GR, m(n-1,0));
+
+    branch(this, x, INT_VAR_SIZE_MIN, INT_VAL_SPLIT_MIN);
+  }
+
+  /// Constructor for cloning \a s
+  MagicSquare(bool share, MagicSquare& s) : Example(share,s), n(s.n) {
+    x.update(this, share, s.x);
+  }
+
+  /// Copy during cloning
+  virtual Space*
+  copy(bool share) {
+    return new MagicSquare(share,*this);
+  }
+  /// Print solution
+  virtual void
+  print(std::ostream& os) {
+    // Matrix-wrapper for the square
+    MiniModel::Matrix<IntVarArray> m(x, n, n);
+    for (int i = 0; i<n; i++) {
+      os << "\t";
+      for (int j = 0; j<n; j++) {
+        os.width(2);
+        os << m(i,j) << "  ";
+      }
+      os << std::endl;
+    }
+  }
+
+};
+
+/** \brief Main-function
+ *  \relates MagicSquare
+ */
+int
+main(int argc, char* argv[]) {
+  SizeOptions opt("MagicSquare");
+  opt.iterations(1);
+  opt.size(7);
+  opt.parse(argc,argv);
+  Example::run<MagicSquare,DFS,SizeOptions>(opt);
+  return 0;
+}
+
+// STATISTICS: example-any
+

data/ext/gecode-2.2.0/examples/minesweeper.cc

@@ -0,0 +1,304 @@
+/* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
+/*
+ *  Main authors:
+ *     Guido Tack <tack@gecode.org>
+ *
+ *  Copyright:
+ *     Guido Tack, 2006
+ *
+ *  Last modified:
+ *     $Date: 2007-11-30 13:58:34 +0100 (Fri, 30 Nov 2007) $ by $Author: tack $
+ *     $Revision: 5524 $
+ *
+ *  This file is part of Gecode, the generic constraint
+ *  development environment:
+ *     http://www.gecode.org
+ *
+ *  Permission is hereby granted, free of charge, to any person obtaining
+ *  a copy of this software and associated documentation files (the
+ *  "Software"), to deal in the Software without restriction, including
+ *  without limitation the rights to use, copy, modify, merge, publish,
+ *  distribute, sublicense, and/or sell copies of the Software, and to
+ *  permit persons to whom the Software is furnished to do so, subject to
+ *  the following conditions:
+ *
+ *  The above copyright notice and this permission notice shall be
+ *  included in all copies or substantial portions of the Software.
+ *
+ *  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+ *  EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+ *  MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+ *  NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+ *  LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+ *  OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+ *  WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+ *
+ */
+
+#include <cctype>
+
+#include "examples/support.hh"
+#include "gecode/minimodel.hh"
+
+namespace {
+  extern const char *specs[];
+  extern const unsigned int n_examples;
+  int spec_size(const char *s);
+  int mineField(const char *s, int n, int i, int j);
+}
+
+/**
+ * \brief %Example: Minesweeper
+ *
+ * This is the classical MineSweeper game.
+ *
+ * The instances are taken from
+ *   http://www.janko.at/Raetsel/Minesweeper/index.htm
+ *
+ * \ingroup ExProblem
+ *
+ */
+class MineSweeper : public Example {
+private:
+  const char *spec;
+  int size;
+  BoolVarArray b;
+
+  /// Access position (\a h,\a w) in the matrix.
+  BoolVar pos(int h, int w) {
+    return b[h*size + w];
+  }
+
+  /// Return the fields in \a m around position (\a x,\a y)
+  BoolVarArgs fieldsAround(MiniModel::Matrix<BoolVarArray>& m,
+                           int x, int y) {
+    int bvsize=0;
+    for (int ix = std::max(0, x-1); ix<=x+1 && ix<size; ix++)
+      for (int iy = std::max(0, y-1); iy<=y+1 && iy<size; iy++)
+        bvsize++;
+    bvsize--; // remove (x,y) itself
+    BoolVarArgs b(bvsize);
+    int count=0;
+    for (int ix = std::max(0, x-1); ix<=x+1 && ix<size; ix++)
+      for (int iy = std::max(0, y-1); iy<=y+1 && iy<size; iy++)
+        if (ix != x || iy != y)
+          b[count++] = m(ix,iy);
+
+    return b;
+  }
+
+public:
+  /// Actual model
+  MineSweeper(const SizeOptions& opt)
+    : spec(specs[opt.size()]), 
+      size(spec_size(spec)),
+      b(this,size*size,0,1) {
+    MiniModel::Matrix<BoolVarArray> m(b, size, size);
+
+    // Initialize matrix and post constraints
+    for (int h=0; h<size; h++)
+      for (int w=0; w<size; w++) {
+        int v = mineField(spec, size, h, w);
+        if (v != -1) {
+          rel(this, m(w, h), IRT_EQ, 0);
+          linear(this, fieldsAround(m, w, h), IRT_EQ, v);
+        }
+      }
+
+    // Install branching
+    branch(this, b, INT_VAR_NONE, INT_VAL_MAX);
+  }
+
+  /// Print solution
+  virtual void
+  print(std::ostream& os) {
+    for (int h = 0; h < size; ++h) {
+      os << '\t';
+      for (int w = 0; w < size; ++w) {
+        int v = mineField(spec, size, h, w);
+        if ( v != -1)
+          os << v << " ";
+        else if (pos(h,w).val() == 1)
+          os << "* ";
+        else
+          os << ". ";
+      }
+      os << std::endl;
+    }
+    os << std::endl;
+  }
+
+  /// Constructor for cloning \a s
+  MineSweeper(bool share, MineSweeper& s) :
+    Example(share,s), spec(s.spec), size(s.size) {
+    b.update(this, share, s.b);
+  }
+  /// Copy space during cloning
+  virtual Space*
+  copy(bool share) {
+    return new MineSweeper(share,*this);
+  }
+
+};
+
+
+/** \brief Main-function
+ *  \relates MineSweeper
+ */
+int
+main(int argc, char* argv[]) {
+  SizeOptions opt("MineSweeper");
+  opt.size(0);
+  opt.parse(argc,argv);
+  if (opt.size() >= n_examples) {
+    std::cerr << "Error: size must be between 0 and "
+              << n_examples-1 << std::endl;
+    return 1;
+  }
+  Example::run<MineSweeper,DFS,SizeOptions>(opt);
+  return 0;
+}
+
+
+namespace  {
+
+  /** \name Minesweeper specifications
+    *
+    * A specification is a square matrix of characters. Alphanumeric characters represent
+    * the number of mines adjacent to that field. Dots represent fields with an unknown number
+    * of mines adjacent to it (or an actual mine).
+    *
+    * \relates MineSweeper
+    */
+  //@{
+  
+  /// The specifications
+  const char* specs[] = {
+      // 0
+      "..2.3."
+      "2....."
+      "..24.3"
+      "1.34.."
+      ".....3"
+      ".3.3..",
+      // 1
+      ".2.211.."
+      "..4.2..2"
+      "2..2..3."
+      "2.22.3.3"
+      "..1...4."
+      "1...2..3"
+      ".2.22.3."
+      "1.1..1.1",
+      // 2
+      "1..2.2.2.."
+      ".32...4..1"
+      "...13...4."
+      "3.1...3..."
+      ".21.1..3.2"
+      ".3.2..2.1."
+      "2..32..2.."
+      ".3...32..3"
+      "..3.33...."
+      ".2.2...22.",
+      // 3
+      "2...3.1."
+      ".5.4...1"
+      "..5..4.."
+      "2...4.5."
+      ".2.4...2"
+      "..5..4.."
+      "2...5.4."
+      ".3.3...2",
+      // 4
+      "0.0.1..11."
+      "1.2.2.22.."
+      "......2..2"
+      ".23.11...."
+      "0......2.1"
+      "...22.1..."
+      ".....3.32."
+      ".5.2...3.1"
+      ".3.1..3..."
+      ".2...12..0",
+      // 5
+      ".21.2.2..."
+      ".4..3...53"
+      "...4.44..3"
+      "4.4..5.6.."
+      "..45....54"
+      "34....55.."
+      "..4.4..5.5"
+      "2..33.6..."
+      "36...3..4."
+      "...4.2.21.",
+      // 6
+      ".32..1.."
+      "....1..3"
+      "3..2...4"
+      ".5...5.."
+      "..6...5."
+      "3...5..4"
+      "2..5...."
+      "..2..34.",
+      // 7
+      ".1.....3."
+      "...343..."
+      "244...443"
+      "...4.4..."
+      ".4.4.3.6."
+      "...4.3..."
+      "123...133"
+      "...322..."
+      ".2.....3.",
+      // 8
+      "......."
+      ".23435."
+      ".1...3."
+      "...5..."
+      ".1...3."
+      ".12234."
+      ".......",
+      // 9
+      "2...2...2"
+      ".4.4.3.4."
+      "..4...1.."
+      ".4.3.3.4."
+      "2.......2"
+      ".5.4.5.4."
+      "..3...3.."
+      ".4.3.5.6."
+      "2...1...2"
+    };
+
+  /// Number of specifications
+  const unsigned int n_examples = sizeof(specs)/sizeof(char*);
+
+  /// Compute the size of a specification
+  int spec_size(const char *s) {
+    int l = std::strlen(s);
+    int res = static_cast<int>(std::sqrt(static_cast<float>(l)));
+    return res;
+  }
+
+  /// Return value at position (\a i,\a j) in the example \a s of size \a n
+  int mineField(const char *s, int n, int i, int j) {
+    assert(spec_size(s) == n);
+    assert(i >= 0 && i < n);
+    assert(j >= 0 && j < n);
+    char c = s[i*n + j];
+    if (!std::isalnum(c))
+      return -1;
+    if (std::isdigit(c))
+      return c - '0';
+    if (std::islower(c))
+      c = static_cast<char>(std::toupper(c));
+    // std::alpha(c) == true
+    int res = (c - 'A') + 10;
+    if (res > n) return 0;
+    else return res;
+  }
+  //@}
+}
+
+// STATISTICS: example-any

data/ext/gecode-2.2.0/examples/money.cc

@@ -0,0 +1,105 @@
+/* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
+/*
+ *  Main authors:
+ *     Christian Schulte <schulte@gecode.org>
+ *
+ *  Copyright:
+ *     Christian Schulte, 2001
+ *
+ *  Last modified:
+ *     $Date: 2007-11-30 13:58:34 +0100 (Fri, 30 Nov 2007) $ by $Author: tack $
+ *     $Revision: 5524 $
+ *
+ *  This file is part of Gecode, the generic constraint
+ *  development environment:
+ *     http://www.gecode.org
+ *
+ *  Permission is hereby granted, free of charge, to any person obtaining
+ *  a copy of this software and associated documentation files (the
+ *  "Software"), to deal in the Software without restriction, including
+ *  without limitation the rights to use, copy, modify, merge, publish,
+ *  distribute, sublicense, and/or sell copies of the Software, and to
+ *  permit persons to whom the Software is furnished to do so, subject to
+ *  the following conditions:
+ *
+ *  The above copyright notice and this permission notice shall be
+ *  included in all copies or substantial portions of the Software.
+ *
+ *  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+ *  EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+ *  MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+ *  NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+ *  LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+ *  OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+ *  WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+ *
+ */
+
+#include "examples/support.hh"
+#include "gecode/minimodel.hh"
+
+/**
+ * \brief %Example: SEND+MORE=MONEY puzzle
+ *
+ * Well-known cryptoarithmetic puzzle.
+ * Henry Dudeney, Strand Magazine, July 1924.
+ *
+ * \ingroup ExProblem
+ *
+ */
+class Money : public Example {
+protected:
+  /// Number of letters
+  static const int nl = 8;
+  /// Array of letters
+  IntVarArray le;
+public:
+  /// Actual model
+  Money(const Options& opt) : le(this,nl,0,9) {
+    IntVar
+      s(le[0]), e(le[1]), n(le[2]), d(le[3]),
+      m(le[4]), o(le[5]), r(le[6]), y(le[7]);
+
+    rel(this, s, IRT_NQ, 0);
+    rel(this, m, IRT_NQ, 0);
+
+    post(this,            1000*s+100*e+10*n+d
+                +         1000*m+100*o+10*r+e
+               == 10000*m+1000*o+100*n+10*e+y,
+         opt.icl());
+
+    distinct(this, le, opt.icl());
+    branch(this, le, INT_VAR_SIZE_MIN, INT_VAL_MIN);
+  }
+  /// Print solution
+  virtual void
+  print(std::ostream& os) {
+    os << "\t" << le << std::endl;
+  }
+
+  /// Constructor for cloning \a s
+  Money(bool share, Money& s) : Example(share,s) {
+    le.update(this, share, s.le);
+  }
+  /// Copy during cloning
+  virtual Space*
+  copy(bool share) {
+    return new Money(share,*this);
+  }
+};
+
+/** \brief Main-function
+ *  \relates Money
+ */
+int
+main(int argc, char* argv[]) {
+  Options opt("SEND+MORE=MONEY");
+  opt.solutions(0);
+  opt.iterations(20000);
+  opt.parse(argc,argv);
+  Example::run<Money,DFS,Options>(opt);
+  return 0;
+}
+
+// STATISTICS: example-any
+