00001 /* 00002 * <partial_mapped_xover.cpp> 00003 * Copyright (C) DOLPHIN Project-Team, INRIA Futurs, 2006-2007 00004 * (C) OPAC Team, LIFL, 2002-2007 00005 * 00006 * Sebastien Cahon, Alexandru-Adrian Tantar 00007 * 00008 * This software is governed by the CeCILL license under French law and 00009 * abiding by the rules of distribution of free software. You can use, 00010 * modify and/ or redistribute the software under the terms of the CeCILL 00011 * license as circulated by CEA, CNRS and INRIA at the following URL 00012 * "http://www.cecill.info". 00013 * 00014 * As a counterpart to the access to the source code and rights to copy, 00015 * modify and redistribute granted by the license, users are provided only 00016 * with a limited warranty and the software's author, the holder of the 00017 * economic rights, and the successive licensors have only limited liability. 00018 * 00019 * In this respect, the user's attention is drawn to the risks associated 00020 * with loading, using, modifying and/or developing or reproducing the 00021 * software by the user in light of its specific status of free software, 00022 * that may mean that it is complicated to manipulate, and that also 00023 * therefore means that it is reserved for developers and experienced 00024 * professionals having in-depth computer knowledge. Users are therefore 00025 * encouraged to load and test the software's suitability as regards their 00026 * requirements in conditions enabling the security of their systems and/or 00027 * data to be ensured and, more generally, to use and operate it in the 00028 * same conditions as regards security. 00029 * The fact that you are presently reading this means that you have had 00030 * knowledge of the CeCILL license and that you accept its terms. 00031 * 00032 * ParadisEO WebSite : http://paradiseo.gforge.inria.fr 00033 * Contact: paradiseo-help@lists.gforge.inria.fr 00034 * 00035 */ 00036 00037 #include <assert.h> 00038 00039 #include <utils/eoRNG.h> 00040 00041 #include "partial_mapped_xover.h" 00042 #include "mix.h" 00043 00044 void PartialMappedXover :: repair (Route & __route, unsigned __cut1, unsigned __cut2) { 00045 00046 unsigned v [__route.size ()] ; // Number of times a cities are visited ... 00047 00048 for (unsigned i = 0 ; i < __route.size () ; i ++) 00049 v [i] = 0 ; 00050 00051 for (unsigned i = 0 ; i < __route.size () ; i ++) 00052 v [__route [i]] ++ ; 00053 00054 std :: vector <unsigned> vert ; 00055 00056 for (unsigned i = 0 ; i < __route.size () ; i ++) 00057 if (! v [i]) 00058 vert.push_back (i) ; 00059 00060 mix (vert) ; 00061 00062 for (unsigned i = 0 ; i < __route.size () ; i ++) 00063 if (i < __cut1 || i >= __cut2) 00064 if (v [__route [i]] > 1) { 00065 __route [i] = vert.back () ; 00066 vert.pop_back () ; 00067 } 00068 } 00069 00070 bool PartialMappedXover :: operator () (Route & __route1, Route & __route2) { 00071 00072 unsigned cut1 = rng.random (__route1.size ()), cut2 = rng.random (__route2.size ()) ; 00073 00074 if (cut2 < cut1) 00075 std :: swap (cut1, cut2) ; 00076 00077 // Between the cuts 00078 for (unsigned i = cut1 ; i < cut2 ; i ++) 00079 std :: swap (__route1 [i], __route2 [i]) ; 00080 00081 // Outside the cuts 00082 repair (__route1, cut1, cut2) ; 00083 repair (__route2, cut1, cut2) ; 00084 00085 __route1.invalidate () ; 00086 __route2.invalidate () ; 00087 00088 return true ; 00089 }