All Downloads are FREE. Search and download functionalities are using the official Maven repository.

org.cicirello.search.operators.permutations.BlockInterchangeMutation Maven / Gradle / Ivy

Go to download

Chips-n-Salsa is a Java library of customizable, hybridizable, iterative, parallel, stochastic, and self-adaptive local search algorithms. The library includes implementations of several stochastic local search algorithms, including simulated annealing, hill climbers, as well as constructive search algorithms such as stochastic sampling. Chips-n-Salsa now also includes genetic algorithms as well as evolutionary algorithms more generally. The library very extensively supports simulated annealing. It includes several classes for representing solutions to a variety of optimization problems. For example, the library includes a BitVector class that implements vectors of bits, as well as classes for representing solutions to problems where we are searching for an optimal vector of integers or reals. For each of the built-in representations, the library provides the most common mutation operators for generating random neighbors of candidate solutions, as well as common crossover operators for use with evolutionary algorithms. Additionally, the library provides extensive support for permutation optimization problems, including implementations of many different mutation operators for permutations, and utilizing the efficiently implemented Permutation class of the JavaPermutationTools (JPT) library. Chips-n-Salsa is customizable, making extensive use of Java's generic types, enabling using the library to optimize other types of representations beyond what is provided in the library. It is hybridizable, providing support for integrating multiple forms of local search (e.g., using a hill climber on a solution generated by simulated annealing), creating hybrid mutation operators (e.g., local search using multiple mutation operators), as well as support for running more than one type of search for the same problem concurrently using multiple threads as a form of algorithm portfolio. Chips-n-Salsa is iterative, with support for multistart metaheuristics, including implementations of several restart schedules for varying the run lengths across the restarts. It also supports parallel execution of multiple instances of the same, or different, stochastic local search algorithms for an instance of a problem to accelerate the search process. The library supports self-adaptive search in a variety of ways, such as including implementations of adaptive annealing schedules for simulated annealing, such as the Modified Lam schedule, implementations of the simpler annealing schedules but which self-tune the initial temperature and other parameters, and restart schedules that adapt to run length.

There is a newer version: 7.0.1
Show newest version
/*
 * Chips-n-Salsa: A library of parallel self-adaptive local search algorithms.
 * Copyright (C) 2002-2022 Vincent A. Cicirello
 *
 * This file is part of Chips-n-Salsa (https://chips-n-salsa.cicirello.org/).
 *
 * Chips-n-Salsa is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * Chips-n-Salsa is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see .
 */

package org.cicirello.search.operators.permutations;

import org.cicirello.math.rand.RandomSampler;
import org.cicirello.permutations.Permutation;
import org.cicirello.search.operators.IterableMutationOperator;
import org.cicirello.search.operators.MutationIterator;
import org.cicirello.search.operators.UndoableMutationOperator;

/**
 * This class implements a block interchange mutation on permutations, where one mutation consists
 * in swapping two randomly chosen non-overlapping "blocks" (i.e., subsequences). The block
 * interchange is chosen uniformly at random from among all possible block interchanges. The two
 * random blocks are not required to be of the same length, but are required to be non-overlapping
 * (i.e., cannot share elements).
 *
 * 

As an example, consider the permutation: p1 = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9]. If a block * interchange swaps the blocks [1, 2] and [5, 6, 7, 8], the result is: [0, 5, 6, 7, 8, 3, 4, 1, 2, * 9]. This mutation operator is related to the {@link BlockMoveMutation}, which swaps a pair of * randomly selected adjacent blocks. * *

The runtime (worst case and average case) of both the {@link #mutate(Permutation) mutate} and * {@link #undo(Permutation) undo} methods is O(n), where n is the length of the permutation. There * are a variety of ways of demonstrating the worst case behavior, one of which is if the exchanged * blocks are at opposite ends of the permutation and are of differing lengths, which would result * in movement of every permutation element. On average, approximately 0.6 n elements are moved by * this mutation operator. * * @author Vincent A. Cicirello, https://www.cicirello.org/ */ public class BlockInterchangeMutation implements UndoableMutationOperator, IterableMutationOperator { // needed to implement undo private final int[] indexes; /** Constructs a BlockInterchangeMutation mutation operator. */ public BlockInterchangeMutation() { indexes = new int[4]; } @Override public final void mutate(Permutation c) { if (c.length() >= 2) { generateIndexes(c.length(), indexes); c.swapBlocks(indexes[0], indexes[1], indexes[2], indexes[3]); } } @Override public final void undo(Permutation c) { if (c.length() >= 2) { c.swapBlocks( indexes[0], indexes[0] + indexes[3] - indexes[2], indexes[3] - indexes[1] + indexes[0], indexes[3]); } } @Override public BlockInterchangeMutation split() { return new BlockInterchangeMutation(); } /** * {@inheritDoc} * *

The worst case runtime of the {@link MutationIterator#hasNext} and the {@link * MutationIterator#setSavepoint} methods of the {@link MutationIterator} created by this method * is O(1). The worst case runtime of the {@link MutationIterator#nextMutant} and {@link * MutationIterator#rollback} methods is O(n), where n is the length of the Permutation. */ @Override public MutationIterator iterator(Permutation p) { return new BlockInterchangeIterator(p); } /* * This package access method allows the window limited version * implemented as a subclass to change how indexes are generated * without modifying the mutate method. */ void generateIndexes(int n, int[] indexes) { // The RandomSampler.sampleInsertion method puts result in sorted order // to begin with. RandomSampler.sampleInsertion(n + 2, 4, indexes); // All index values generated by above are unique. // However, block size 1 should be allowed, which would require // duplicated indexes. We handle this by passing n+2 (above). An index // equal to n is treated as a duplicate of indexes[0] (i.e., a block of size 1). // An index equal to n+1 is likewise treated as a duplicate to define the right // block to be of size 1. if (indexes[3] == n) { indexes[3] = indexes[2]; indexes[2] = indexes[1]; indexes[1] = indexes[0]; } else if (indexes[2] == n) { indexes[3] = indexes[2] = indexes[1]; indexes[1] = indexes[0]; } else if (indexes[3] == n + 1) { indexes[3] = indexes[2]; } } }





© 2015 - 2025 Weber Informatics LLC | Privacy Policy