// Licensed to the Apache Software Foundation (ASF) under one // or more contributor license agreements. See the NOTICE file // distributed with this work for additional information // regarding copyright ownership. The ASF licenses this file // to you under the Apache License, Version 2.0 (the // "License"); you may not use this file except in compliance // with the License. You may obtain a copy of the License at // // http://www.apache.org/licenses/LICENSE-2.0 // // Unless required by applicable law or agreed to in writing, // software distributed under the License is distributed on an // "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY // KIND, either express or implied. See the License for the // specific language governing permissions and limitations // under the License. #pragma once #include #include #include #include #include #include namespace arrow { namespace internal { template > std::vector ArgSort(const std::vector& values, Cmp&& cmp = {}) { std::vector indices(values.size()); std::iota(indices.begin(), indices.end(), 0); std::sort(indices.begin(), indices.end(), [&](int64_t i, int64_t j) -> bool { return cmp(values[i], values[j]); }); return indices; } template size_t Permute(const std::vector& indices, std::vector* values) { if (indices.size() <= 1) { return indices.size(); } // mask indicating which of values are in the correct location std::vector sorted(indices.size(), false); size_t cycle_count = 0; for (auto cycle_start = sorted.begin(); cycle_start != sorted.end(); cycle_start = std::find(cycle_start, sorted.end(), false)) { ++cycle_count; // position in which an element belongs WRT sort auto sort_into = static_cast(cycle_start - sorted.begin()); if (indices[sort_into] == sort_into) { // trivial cycle sorted[sort_into] = true; continue; } // resolve this cycle const auto end = sort_into; for (int64_t take_from = indices[sort_into]; take_from != end; take_from = indices[sort_into]) { std::swap(values->at(sort_into), values->at(take_from)); sorted[sort_into] = true; sort_into = take_from; } sorted[sort_into] = true; } return cycle_count; } } // namespace internal } // namespace arrow