permlib
0.2.6
Library for permutation computations
|
00001 // --------------------------------------------------------------------------- 00002 // 00003 // This file is part of PermLib. 00004 // 00005 // Copyright (c) 2009-2011 Thomas Rehn <thomas@carmen76.de> 00006 // All rights reserved. 00007 // 00008 // Redistribution and use in source and binary forms, with or without 00009 // modification, are permitted provided that the following conditions 00010 // are met: 00011 // 1. Redistributions of source code must retain the above copyright 00012 // notice, this list of conditions and the following disclaimer. 00013 // 2. Redistributions in binary form must reproduce the above copyright 00014 // notice, this list of conditions and the following disclaimer in the 00015 // documentation and/or other materials provided with the distribution. 00016 // 3. The name of the author may not be used to endorse or promote products 00017 // derived from this software without specific prior written permission. 00018 // 00019 // THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR 00020 // IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES 00021 // OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. 00022 // IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT, 00023 // INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT 00024 // NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, 00025 // DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY 00026 // THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT 00027 // (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF 00028 // THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 00029 // 00030 // --------------------------------------------------------------------------- 00031 00032 00033 #ifndef SETSTABILIZEREFINEMENT_H_ 00034 #define SETSTABILIZEREFINEMENT_H_ 00035 00036 #include <permlib/predicate/pointwise_stabilizer_predicate.h> 00037 00038 #include <permlib/search/partition/partition.h> 00039 #include <permlib/search/partition/refinement.h> 00040 00041 #include <algorithm> 00042 #include <boost/dynamic_bitset.hpp> 00043 #include <boost/foreach.hpp> 00044 00045 namespace permlib { 00046 namespace partition { 00047 00049 template<class PERM> 00050 class SetStabilizeRefinement : public Refinement<PERM> { 00051 public: 00053 template<class InputIterator> 00054 SetStabilizeRefinement(unsigned long n, InputIterator begin, InputIterator end); 00055 00056 virtual unsigned int apply(Partition& pi) const; 00057 00058 virtual bool init(Partition& pi); 00059 private: 00060 std::vector<unsigned long> toStab; 00061 }; 00062 00063 template<class PERM> 00064 template<class InputIterator> 00065 SetStabilizeRefinement<PERM>::SetStabilizeRefinement(unsigned long n, InputIterator begin, InputIterator end) 00066 : Refinement<PERM>(n, Default), toStab(begin, end) 00067 { 00068 std::sort(toStab.begin(), toStab.end()); 00069 PERMLIB_DEBUG(print_iterable(toStab.begin(), toStab.end(), 0, "to stab");) 00070 } 00071 00072 template<class PERM> 00073 unsigned int SetStabilizeRefinement<PERM>::apply(Partition& pi) const { 00074 BOOST_ASSERT( this->initialized() ); 00075 unsigned int ret = 0; 00076 BOOST_FOREACH(unsigned int cell, Refinement<PERM>::m_cellPairs) { 00077 PERMLIB_DEBUG(std::cout << "apply set stab " << cell << std::endl;) 00078 if (pi.intersect(toStab.begin(), toStab.end(), cell)) 00079 ++ret; 00080 } 00081 return ret; 00082 } 00083 00084 template<class PERM> 00085 bool SetStabilizeRefinement<PERM>::init(Partition& pi) { 00086 for (unsigned int c = 0; c < pi.cells(); ++c) { 00087 if (pi.intersect(toStab.begin(), toStab.end(), c)) 00088 Refinement<PERM>::m_cellPairs.push_back(c); 00089 } 00090 if (!Refinement<PERM>::m_cellPairs.empty()) { 00091 typename Refinement<PERM>::RefinementPtr ref(new SetStabilizeRefinement<PERM>(*this)); 00092 Refinement<PERM>::m_backtrackRefinements.push_back(ref); 00093 return true; 00094 } 00095 return false; 00096 } 00097 00098 } 00099 } 00100 00101 #endif // -- SETSTABILIZEREFINEMENT_H_