gecoder-with-gecode 0.7.1

data/ext/gecode-1.3.1/gecode/set/view/complement.icc

@@ -0,0 +1,525 @@
+/*
+ *  Main authors:
+ *     Guido Tack <tack@gecode.org>
+ *
+ *  Contributing authors:
+ *     Christian Schulte <schulte@gecode.org>
+ *
+ *  Copyright:
+ *     Guido Tack, 2004
+ *     Christian Schulte, 2004
+ *
+ *  Last modified:
+ *     $Date: 2006-08-24 11:25:05 +0200 (Thu, 24 Aug 2006) $ by $Author: schulte $
+ *     $Revision: 3559 $
+ *
+ *  This file is part of Gecode, the generic constraint
+ *  development environment:
+ *     http://www.gecode.org
+ *
+ *  See the file "LICENSE" for information on usage and
+ *  redistribution of this file, and for a
+ *     DISCLAIMER OF ALL WARRANTIES.
+ *
+ */
+
+namespace Gecode { namespace Set {
+
+  template <class View>
+  forceinline
+  ComplementView<View>::ComplementView(void) {}
+
+  template <class View>
+  forceinline
+  ComplementView<View>::ComplementView(View& s0)
+    : DerivedViewBase<View>(s0) {}
+
+  template <class View>
+  forceinline ModEvent
+  ComplementView<View>::me_negateset(ModEvent me) {
+    switch(me) {
+    case ME_SET_LUB : return ME_SET_GLB;
+    case ME_SET_GLB : return ME_SET_LUB;
+    case ME_SET_CLUB : return ME_SET_CGLB;
+    case ME_SET_CGLB : return ME_SET_CLUB;
+    default: return me;
+    }
+  }
+
+  template <class View>
+  forceinline PropCond
+  ComplementView<View>::pc_negateset(PropCond pc) {
+    switch(pc) {
+    case PC_SET_CLUB  : return PC_SET_CGLB;
+    case PC_SET_CGLB  : return PC_SET_CLUB;
+    default: return pc;
+    }
+  }
+
+  template <class View>
+  forceinline bool
+  ComplementView<View>::assigned(void) const { return view.assigned(); }
+
+  template <class View>
+  forceinline unsigned int
+  ComplementView<View>::glbSize(void) const {
+    return Limits::Set::card_max - view.lubSize();
+  }
+
+  template <class View>
+  forceinline unsigned int
+  ComplementView<View>::lubSize(void) const {
+    return Limits::Set::card_max - view.glbSize();
+  }
+
+  template <class View>
+  forceinline unsigned int
+  ComplementView<View>::unknownSize(void) const {
+    return lubSize() - glbSize();
+  }
+
+  template <class View>
+  forceinline bool
+  ComplementView<View>::contains(int n) const { return view.notContains(n); }
+
+  template <class View>
+  forceinline bool
+  ComplementView<View>::notContains(int n) const { return view.contains(n); }
+
+  template <class View>
+  forceinline unsigned int
+  ComplementView<View>::cardMin() const {
+    return Limits::Set::card_max - view.cardMax();
+  }
+
+  template <class View>
+  forceinline unsigned int
+  ComplementView<View>::cardMax() const {
+    return Limits::Set::card_max - view.cardMin();
+  }
+
+  template <class View>
+  forceinline int
+  ComplementView<View>::lubMin() const {
+    GlbRanges<View> lb(view);
+    RangesCompl<GlbRanges<View> > lbc(lb);
+    if (lbc()) {
+      return lbc.min();
+    } else {
+      return BndSet::MIN_OF_EMPTY;
+    }
+  }
+
+  template <class View>
+  forceinline int
+  ComplementView<View>::lubMax() const {
+    GlbRanges<View> lb(view);
+    RangesCompl<GlbRanges<View> > lbc(lb);
+    if (lbc()) {
+      while(lbc()) ++lbc;
+      return lbc.max();
+    } else {
+      return BndSet::MAX_OF_EMPTY;
+    }
+  }
+
+  template <class View>
+  forceinline int
+  ComplementView<View>::glbMin() const {
+    LubRanges<View> ub(view);
+    RangesCompl<LubRanges<View> > ubc(ub);
+    if (ubc()) {
+      return ubc.min();
+    } else {
+      return BndSet::MIN_OF_EMPTY;
+    }
+  }
+
+  template <class View>
+  forceinline int
+  ComplementView<View>::glbMax() const {
+    LubRanges<View> ub(view);
+    RangesCompl<LubRanges<View> > ubc(ub);
+    if (ubc()) {
+      while(ubc()) ++ubc;
+      return ubc.max();
+    } else {
+      return BndSet::MAX_OF_EMPTY;
+    }
+  }
+
+  template <class View>
+  forceinline ModEvent
+  ComplementView<View>::cardMin(Space* home, unsigned int c) {
+    if (c < Limits::Set::card_max)
+      return me_negateset(view.cardMax(home, Limits::Set::card_max - c));
+    return ME_SET_NONE;
+  }
+
+  template <class View>
+  forceinline ModEvent
+  ComplementView<View>::cardMax(Space* home, unsigned int c) {
+    if (c < Limits::Set::card_max)
+      return me_negateset(view.cardMin(home, Limits::Set::card_max - c));
+    return ME_SET_NONE;
+  }
+
+  template <class View>
+  forceinline ModEvent
+  ComplementView<View>::include(Space* home, int c) {
+    return me_negateset((view.exclude(home, c)));
+  }
+
+  template <class View>
+  forceinline ModEvent
+  ComplementView<View>::exclude(Space* home, int c) {
+    return me_negateset((view.include(home, c)));
+  }
+
+  template <class View>
+  forceinline ModEvent
+  ComplementView<View>::intersect(Space* home, int c) {
+    Iter::Ranges::Singleton si(c,c);
+    return me_negateset((view.includeI(home, si)));
+  }
+
+  template <class View>
+  forceinline ModEvent
+  ComplementView<View>::intersect(Space* home, int i, int j) {
+    Iter::Ranges::Singleton si(i,j);
+    return me_negateset((view.includeI(home, si)));
+  }
+
+  template <class View>
+  forceinline ModEvent
+  ComplementView<View>::include(Space* home, int j, int k) {
+    return me_negateset(view.exclude(home,j,k));
+  }
+
+  template <class View>
+  forceinline ModEvent
+  ComplementView<View>::exclude(Space* home, int j, int k) {
+    return me_negateset(view.include(home,j,k));
+  }
+
+  template <class View>
+  template <class I> ModEvent
+  ComplementView<View>::excludeI(Space* home,I& iter) {
+    return me_negateset(view.includeI(home,iter));
+  }
+
+  template <class View>
+  template <class I> ModEvent
+  ComplementView<View>::includeI(Space* home,I& iter) {
+    return me_negateset(view.excludeI(home,iter));
+  }
+
+  template <class View>
+  template <class I> ModEvent
+  ComplementView<View>::intersectI(Space* home,I& iter) {
+    RangesCompl<I> c(iter);
+    return me_negateset(view.includeI(home,c));
+  }
+
+  template <class View>
+  forceinline void
+  ComplementView<View>::subscribe(Space* home, Propagator* p, PropCond pc,
+				  bool process) {
+    view.subscribe(home,p, pc_negateset(pc),process);
+  }
+
+  template <class View>
+  forceinline void
+  ComplementView<View>::cancel(Space* home, Propagator* p, PropCond pc) {
+    view.cancel(home,p, pc_negateset(pc));
+  }
+
+  template <class View>
+  forceinline ModEvent
+  ComplementView<View>::pme(const Propagator* p) {
+    return me_negateset(View::pme(p));
+  }
+
+  template <class View>
+  forceinline PropModEvent
+  ComplementView<View>::pme(ModEvent me) {
+    return me_negateset(View::pme(me));
+  }
+
+  template <class View>
+  forceinline void
+  ComplementView<View>::update(Space* home, bool share, 
+			       ComplementView& y) {
+    view.update(home,share,y.view);
+  }
+
+  /*
+   * Specialization for double negation
+   *
+   */
+
+  template <class View>
+  forceinline
+  ComplementView<ComplementView<View> >::ComplementView(void) {}
+
+  template <class View>
+  forceinline
+  ComplementView<ComplementView<View> >::
+  ComplementView(ComplementView<View>& s0)
+    : View(s0.base()) {}
+
+
+  /**
+   * \brief %Range iterator for least upper bound of complement set views
+   * \ingroup TaskActorSetView
+   */
+  template <class View>
+  class LubRanges<ComplementView<View> > {
+  private:
+    GlbRanges<View> lb;
+    RangesCompl<GlbRanges<View> > lbc;
+  public:
+    /// \name Constructors and initialization
+    //@{
+    /// Default constructor
+    LubRanges(void) {}
+    /// Initialize with ranges for view \a x
+    LubRanges(const ComplementView<View>& x);
+    /// Initialize with ranges for view \a x
+    void init(const ComplementView<View>& x);
+
+    /// \name Iteration control
+    //@{
+    /// Test whether iterator is still at a range or done
+    bool operator()(void) const;
+    /// Move iterator to next range (if possible)
+    void operator++(void);
+    //@}
+
+    /// \name Range access
+    //@{
+    /// Return smallest value of range
+    int min(void) const;
+    /// Return largest value of range
+    int max(void) const;
+    /// Return width of ranges (distance between minimum and maximum)
+    unsigned int width(void) const;
+    //@}
+  };
+
+  template <class View>
+  forceinline
+  LubRanges<ComplementView<View> >::LubRanges(const ComplementView<View>& s)
+    : lb(s.base()), lbc(lb) {}
+
+  template <class View>
+  forceinline void
+  LubRanges<ComplementView<View> >::init(const ComplementView<View>& s) {
+    lb.init(s.base());
+    lbc.init(lb);
+  }
+
+  template <class View>
+  forceinline bool
+  LubRanges<ComplementView<View> >::operator()(void) const { return lbc(); }
+
+  template <class View>
+  forceinline void
+  LubRanges<ComplementView<View> >::operator++(void) { return ++lbc; }
+
+  template <class View>
+  forceinline int
+  LubRanges<ComplementView<View> >::min(void) const { return lbc.min(); }
+
+  template <class View>
+  forceinline int
+  LubRanges<ComplementView<View> >::max(void) const { return lbc.max(); }
+
+  template <class View>
+  forceinline unsigned int
+  LubRanges<ComplementView<View> >::width(void) const { return lbc.width(); }
+
+  /**
+   * \brief Range iterator for the least upper bound of double-complement-views
+   *
+   * This class provides (by specialization) a range iterator 
+   * for the least upper bounds of complements of complement set views.
+   *
+   * \ingroup TaskActorSet
+   */  
+  template <class View>
+  class LubRanges<ComplementView<ComplementView<View> > > :
+  public LubRanges<View> {
+  public:
+    /// \name Constructors and initialization
+    //@{
+    /// Default constructor
+    LubRanges(void) {}
+    /// Initialize with ranges for view \a x
+    LubRanges(const ComplementView<ComplementView<View> >& x);
+    /// Initialize with ranges for view \a x
+    void init(const ComplementView<ComplementView<View> >& x);
+    //@}
+  };
+
+  template <class View>
+  forceinline
+  LubRanges<ComplementView<ComplementView<View> > >::
+  LubRanges(const ComplementView<ComplementView<View> >& x) :
+  LubRanges<View>(x) {}
+
+  template <class View>
+  forceinline void
+  LubRanges<ComplementView<ComplementView<View> > >::
+  init(const ComplementView<ComplementView<View> >& x) {
+    LubRanges<View>::init(x);
+  }
+
+  /**
+   * \brief %Range iterator for greatest lower bound of complement set views
+   * \ingroup TaskActorSetView
+   */
+  template <class View>
+  class GlbRanges<ComplementView<View> > {
+  private:
+    LubRanges<View> ub;
+    RangesCompl<LubRanges<View> > ubc;
+  public:
+    /// \name Constructors and initialization
+    //@{
+    /// Default constructor
+    GlbRanges(void) {}
+    /// Initialize with ranges for view \a x
+    GlbRanges(const ComplementView<View> & x);
+    /// Initialize with ranges for view \a x
+    void init(const ComplementView<View> & x);
+
+    /// \name Iteration control
+    //@{
+    /// Test whether iterator is still at a range or done
+    bool operator()(void) const;
+    /// Move iterator to next range (if possible)
+    void operator++(void);
+    //@}
+
+    /// \name Range access
+    //@{
+    /// Return smallest value of range
+    int min(void) const;
+    /// Return largest value of range
+    int max(void) const;
+    /// Return width of ranges (distance between minimum and maximum)
+    unsigned int width(void) const;
+    //@}
+  };
+
+  template <class View>
+  forceinline
+  GlbRanges<ComplementView<View> >::GlbRanges(const ComplementView<View> & s)
+    : ub(s.base()), ubc(ub) {}
+
+  template <class View>
+  forceinline void
+  GlbRanges<ComplementView<View> >::init(const ComplementView<View> & s) {
+    ub.init(s.base());
+    ubc.init(ub);
+  }
+
+  template <class View>
+  forceinline bool
+  GlbRanges<ComplementView<View> >::operator()(void) const { return ubc(); }
+
+  template <class View>
+  forceinline void
+  GlbRanges<ComplementView<View> >::operator++(void) { return ++ubc; }
+
+  template <class View>
+  forceinline int
+  GlbRanges<ComplementView<View> >::min(void) const { return ubc.min(); }
+
+  template <class View>
+  forceinline int
+  GlbRanges<ComplementView<View> >::max(void) const { return ubc.max(); }
+
+  template <class View>
+  forceinline unsigned int
+  GlbRanges<ComplementView<View> >::width(void) const { return ubc.width(); }
+  
+  /**
+   * \brief Range iterator for the greatest lower bound of double-complement-views
+   *
+   * This class provides (by specialization) a range iterator 
+   * for the greatest lower bounds of complements of complement set views.
+   *
+   * \ingroup TaskActorSet
+   */  
+  template <class View>
+  class GlbRanges<ComplementView<ComplementView<View> > > :
+  public GlbRanges<View> {
+  public:
+    /// \name Constructors and initialization
+    //@{
+    /// Default constructor
+    GlbRanges(void) {}
+    /// Initialize with ranges for view \a x
+    GlbRanges(const ComplementView<ComplementView<View> >& x);
+    /// Initialize with ranges for view \a x
+    void init(const ComplementView<ComplementView<View> >& x);
+    //@}
+  };
+
+  template <class View>
+  forceinline
+  GlbRanges<ComplementView<ComplementView<View> > >::
+  GlbRanges(const ComplementView<ComplementView<View> >& x) :
+  GlbRanges<View>(x) {}
+
+  template <class View>
+  forceinline void
+  GlbRanges<ComplementView<ComplementView<View> > >::
+  init(const ComplementView<ComplementView<View> >& x) {
+    GlbRanges<View>::init(x);
+  }
+
+}
+
+
+  /*
+   * Testing
+   *
+   */
+  template <class View>
+  forceinline bool
+  same(const Set::ComplementView<View>& x, 
+       const Set::ComplementView<View>& y) {
+    return same(x.base(),y.base());
+  }
+  template <class View>
+  forceinline bool
+  before(const Set::ComplementView<View>& x, 
+	 const Set::ComplementView<View>& y) {
+    return before(x.base(),y.base());
+  }
+  template <class View>
+  forceinline bool
+  same(const Set::ComplementView<Set::ComplementView<View> >& x, 
+       const Set::ComplementView<Set::ComplementView<View> >& y) {
+    return same(x,y);
+  }
+  template <class View>
+  forceinline bool
+  before(const Set::ComplementView<Set::ComplementView<View> >& x, 
+	 const Set::ComplementView<Set::ComplementView<View> >& y) {
+    return before(x,y);
+  }
+
+}
+
+template <class View>
+forceinline
+std::ostream&
+operator<<(std::ostream& os, const Gecode::Set::ComplementView<View>& s) {
+  return os << "(" << s.base() << ")^C";
+}
+
+// STATISTICS: set-var