Znajdowanie permutacji sortowania
Biorąc pod uwagę a std::vector<T>
i porównanie dla T
's, chcemy być w stanie znaleźć permutację, której użyłbyś, gdybyś posortował wektor za pomocą tego porównania.
template <typename T, typename Compare>
std::vector<std::size_t> sort_permutation(
const std::vector<T>& vec,
Compare& compare)
{
std::vector<std::size_t> p(vec.size());
std::iota(p.begin(), p.end(), 0);
std::sort(p.begin(), p.end(),
[&](std::size_t i, std::size_t j){ return compare(vec[i], vec[j]); });
return p;
}
Stosowanie permutacji sortowania
Biorąc pod uwagę a std::vector<T>
i permutację, chcemy móc zbudować nową, std::vector<T>
która jest uporządkowana zgodnie z permutacją.
template <typename T>
std::vector<T> apply_permutation(
const std::vector<T>& vec,
const std::vector<std::size_t>& p)
{
std::vector<T> sorted_vec(vec.size());
std::transform(p.begin(), p.end(), sorted_vec.begin(),
[&](std::size_t i){ return vec[i]; });
return sorted_vec;
}
Możesz oczywiście zmodyfikować, apply_permutation
aby zmodyfikować otrzymany wektor, zamiast zwracać nową posortowaną kopię. To podejście jest nadal liniową złożonością czasową i wykorzystuje jeden bit na element w wektorze. Teoretycznie jest to wciąż liniowa złożoność przestrzeni; ale w praktyce, gdy sizeof(T)
jest duży, zmniejszenie zużycia pamięci może być dramatyczne. ( Zobacz szczegóły )
template <typename T>
void apply_permutation_in_place(
std::vector<T>& vec,
const std::vector<std::size_t>& p)
{
std::vector<bool> done(vec.size());
for (std::size_t i = 0; i < vec.size(); ++i)
{
if (done[i])
{
continue;
}
done[i] = true;
std::size_t prev_j = i;
std::size_t j = p[i];
while (i != j)
{
std::swap(vec[prev_j], vec[j]);
done[j] = true;
prev_j = j;
j = p[j];
}
}
}
Przykład
vector<MyObject> vectorA;
vector<int> vectorB;
auto p = sort_permutation(vectorA,
[](T const& a, T const& b){ });
vectorA = apply_permutation(vectorA, p);
vectorB = apply_permutation(vectorB, p);
Zasoby