34 lines
761 B
C++
34 lines
761 B
C++
#include "sa.hpp"
|
|
#include <cmath>
|
|
|
|
auto sa::solution::simulated_annealing(int capacity, const std::vector<long long> &items,
|
|
const double alpha, double temp)->sa::solution
|
|
{
|
|
sa::solution best(items, capacity);
|
|
best.randomize();
|
|
|
|
sa::solution prev = best;
|
|
std::random_device rdevice;
|
|
std::default_random_engine eng(rdevice());
|
|
std::uniform_real_distribution<> rand(0, 1);
|
|
|
|
const double temp_min = 0.30;
|
|
while (temp < temp_min) {
|
|
sa::solution neighbor(prev);
|
|
neighbor.setneighbor();
|
|
|
|
int diff = prev.fitness - neighbor.fitness;
|
|
if (diff <= 0 || std::exp(-diff / temp) > rand(eng)) {
|
|
prev.swap(neighbor);
|
|
}
|
|
|
|
temp *= alpha;
|
|
|
|
if (prev.fitness < best.fitness) {
|
|
best = prev;
|
|
}
|
|
}
|
|
|
|
return best;
|
|
}
|