C++ 在C+中按代理排序(或:按另一个容器的内容排序)+;
我有一组数据,分成两个数组(我们称它们为C++ 在C+中按代理排序(或:按另一个容器的内容排序)+;,c++,sorting,stl,C++,Sorting,Stl,我有一组数据,分成两个数组(我们称它们为data和keys)。也就是说,对于索引为i的任何给定项,我可以使用data[i]访问该项的数据,并使用keys[i]访问该项的键。我无法更改此结构(例如,将密钥和数据交织到单个数组中),因为我需要将数据数组传递给需要特定数据布局的库函数 如何根据键数组的内容对两个数组进行排序(最好使用标准库函数)?您可以使用映射: int main() { vector<int> keys; vector<string> data;
data
和keys
)。也就是说,对于索引为i
的任何给定项,我可以使用data[i]
访问该项的数据,并使用keys[i]
访问该项的键。我无法更改此结构(例如,将密钥和数据交织到单个数组中),因为我需要将数据数组传递给需要特定数据布局的库函数
如何根据键
数组的内容对两个数组进行排序(最好使用标准库函数)?您可以使用映射:
int main() {
vector<int> keys;
vector<string> data;
keys.push_back(5); data.push_back("joe");
keys.push_back(2); data.push_back("yaochun");
keys.push_back(1); data.push_back("holio");
// load the keys and data to the map (they will automatically be inserted in sorted order by key)
map<int, string> sortedVals;
for(int i = 0; i < (int)keys.size(); ++i) {
sortedVals[keys[i]] = data[i];
}
// copy the map values back to vectors
int ndx=0;
for(map<int, string>::iterator it = sortedVals.begin(); it != sortedVals.end(); ++it) {
keys[ndx] = it->first;
data[ndx] = it->second;
++ndx;
}
// verify
for(int i = 0; i < (int)keys.size(); ++i) {
cout<<keys[i]<<" "<<data[i]<<endl;
}
return 0;
}
创建包含两个数组索引的对象向量。定义运算符您可以使用函子进行排序,例如:
template <class T>
struct IndexFunctor {
IndexFunctor(const std::vector<T>& v_) : v(v_) {}
bool operator ()(int a, int b) const {
return v[a] < v[b];
}
const std::vector<T>& v;
};
template <class K, class D>
void SortByKeys(std::vector<K>& keys, std::vector<D>& data) {
// Calculate the valid order inside a permutation array p.
const int n = static_cast<int>(keys.size());
std::vector<int> p(n);
for (int i = 0; i < n; ++i) p[i] = i;
std::sort(p.begin(), p.end(), IndexFunctor(keys));
// Reorder the keys and data in temporary arrays, it cannot be done in place.
std::vector<K> aux_keys(n);
std::vector<D> aux_data(n);
for (int i = 0; i < n; ++i) {
aux_keys[i] = keys[p[i]];
aux_data[i] = data[p[i]];
}
// Assign the ordered vectors by swapping, it should be faster.
keys.swap(aux_keys);
data.swap(aux_data);
}
模板
结构索引器{
IndexFunctor(const std::vector&v_uu):v(v_u){
布尔运算符()(int a,int b)常量{
返回v[a]
这个问题真的让我深思。我提出了一个解决方案,它利用一些C++0x特性来获得一个非常类似STL的并行排序算法。为了“就地”执行排序,我必须编写一个back\u remove\u迭代器
,作为back\u insert\u迭代器
的对应项,以允许算法读取和写入同一容器。你可以跳过这些部分,直接进入有趣的内容
我没有对它进行任何核心测试,但它在时间和空间上似乎都相当有效,主要是因为使用了std::move()
,以防止不必要的复制
#include <algorithm>
#include <iostream>
#include <string>
#include <vector>
//
// An input iterator that removes elements from the back of a container.
// Provided only because the standard library neglects one.
//
template<class Container>
class back_remove_iterator :
public std::iterator<std::input_iterator_tag, void, void, void, void> {
public:
back_remove_iterator() : container(0) {}
explicit back_remove_iterator(Container& c) : container(&c) {}
back_remove_iterator& operator=
(typename Container::const_reference value) { return *this; }
typename Container::value_type operator*() {
typename Container::value_type value(container->back());
container->pop_back();
return value;
} // operator*()
back_remove_iterator& operator++() { return *this; }
back_remove_iterator operator++(int) { return *this; }
Container* container;
}; // class back_remove_iterator
//
// Equivalence operator for back_remove_iterator. An iterator compares equal
// to the end iterator either if it is default-constructed or if its
// container is empty.
//
template<class Container>
bool operator==(const back_remove_iterator<Container>& a,
const back_remove_iterator<Container>& b) {
return !a.container ? !b.container || b.container->empty() :
!b.container ? !a.container || a.container->empty() :
a.container == b.container;
} // operator==()
//
// Inequivalence operator for back_remove_iterator.
//
template<class Container>
bool operator!=(const back_remove_iterator<Container>& a,
const back_remove_iterator<Container>& b) {
return !(a == b);
} // operator!=()
//
// A handy way to default-construct a back_remove_iterator.
//
template<class Container>
back_remove_iterator<Container> back_remover() {
return back_remove_iterator<Container>();
} // back_remover()
//
// A handy way to construct a back_remove_iterator.
//
template<class Container>
back_remove_iterator<Container> back_remover(Container& c) {
return back_remove_iterator<Container>(c);
} // back_remover()
//
// A comparison functor that sorts std::pairs by their first element.
//
template<class A, class B>
struct sort_pair_by_first {
bool operator()(const std::pair<A, B>& a, const std::pair<A, B>& b) {
return a.first < b.first;
} // operator()()
}; // struct sort_pair_by_first
//
// Performs a parallel sort of the ranges [keys_first, keys_last) and
// [values_first, values_last), preserving the ordering relation between
// values and keys. Sends key and value output to keys_out and values_out.
//
// This works by building a vector of std::pairs, sorting them by the key
// element, then returning the sorted pairs as two separate sequences. Note
// the use of std::move() for a vast performance improvement.
//
template<class A, class B, class I, class J, class K, class L>
void parallel_sort(I keys_first, I keys_last, J values_first, J values_last,
K keys_out, L values_out) {
typedef std::vector< std::pair<A, B> > Pairs;
Pairs sorted;
while (keys_first != keys_last)
sorted.push_back({std::move(*keys_first++), std::move(*values_first++)});
std::sort(sorted.begin(), sorted.end(), sort_pair_by_first<A, B>());
for (auto i = sorted.begin(); i != sorted.end(); ++i)
*keys_out++ = std::move(i->first),
*values_out++ = std::move(i->second);
} // parallel_sort()
int main(int argc, char** argv) {
//
// There is an ordering relation between keys and values,
// but the sets still need to be sorted. Sounds like a job for...
//
std::vector<int> keys{0, 3, 1, 2};
std::vector<std::string> values{"zero", "three", "one", "two"};
//
// parallel_sort! Unfortunately, the key and value types do need to
// be specified explicitly. This could be helped with a utility
// function that accepts back_remove_iterators.
//
parallel_sort<int, std::string>
(back_remover(keys), back_remover<std::vector<int>>(),
back_remover(values), back_remover<std::vector<std::string>>(),
std::back_inserter(keys), std::back_inserter(values));
//
// Just to prove that the mapping is preserved.
//
for (unsigned int i = 0; i < keys.size(); ++i)
std::cout << keys[i] << ": " << values[i] << '\n';
return 0;
} // main()
#包括
#包括
#包括
#包括
//
//从容器后面删除元素的输入迭代器。
//仅因标准库忽略了一个而提供。
//
模板
类返回\u移除\u迭代器:
公共标准:迭代器{
公众:
back_remove_迭代器():容器(0){}
显式back_remove_迭代器(Container&c):Container&c{}
后退\删除\迭代器和运算符=
(typename容器::const_引用值){return*this;}
typename容器::值\类型运算符*(){
typename容器::值\类型值(容器->返回());
容器->弹出返回();
返回值;
}//运算符*()
back_remove_迭代器和运算符++(){return*this;}
back_remove_迭代器运算符++(int){return*this;}
集装箱*集装箱;
}; // 类back\u remove\u迭代器
//
//back_remove_迭代器的等价运算符。迭代器比较相等
//如果它是默认构造的,或者
//容器是空的。
//
模板
布尔运算符==(常量返回\u删除\u迭代器&a,
const back_remove_迭代器&b){
return!a.container?!b.container | | b.container->empty():
!b.container?!a.container | | a.container->empty():
a、 容器==b.容器;
}//运算符==()
//
//back_remove_迭代器的非等价运算符。
//
模板
接线员=(const back_remove_迭代器&a,
const back_remove_迭代器&b){
返回!(a==b);
}//接线员=()
//
//默认构造back_remove_迭代器的简便方法。
//
模板
back\u remove\u迭代器back\u remover(){
返回_remove_迭代器();
}//后卸料器()
//
//构造back_remove_迭代器的简便方法。
//
模板
后退移除迭代器后退移除器(容器和c){
返回\删除\迭代器(c);
}//后卸料器()
//
//按std::pairs的第一个元素排序的比较函子。
//
模板
结构排序\u对\u按\u优先{
布尔运算符()(常数std::pair&a,常数std::pair&b){
返回a.firstPairs;
成对排序;
while(keys\u first!=keys\u last)
已排序。向后推({std::move(*keys_first++),std::move(*values_first++)});
排序(sorted.begin()、sorted.end()、sort_pair_by_first());
for(自动i=排序的.begin();i!=排序的.end();++i)
*按键_out++=std::move(i->first),
*值_out++=std::move(i->second);
}//并行排序()
int main(int argc,字符**argv){
//
//键和值之间存在顺序关系,
//但是布景还需要整理。听起来像是。。。
//
向量键{0,3,1,2};
向量值{“零”、“三”、“一”、“二”};
//
//并行排序!不幸的是,键和值类型确实需要
//被显式指定。这可以通过实用程序得到帮助
//接受back\u remove\u迭代器的函数。
//
并行排序
(后拆卸器(键),后拆卸器(),
背面去除器(值),背面
template <class T>
struct IndexFunctor {
IndexFunctor(const std::vector<T>& v_) : v(v_) {}
bool operator ()(int a, int b) const {
return v[a] < v[b];
}
const std::vector<T>& v;
};
template <class K, class D>
void SortByKeys(std::vector<K>& keys, std::vector<D>& data) {
// Calculate the valid order inside a permutation array p.
const int n = static_cast<int>(keys.size());
std::vector<int> p(n);
for (int i = 0; i < n; ++i) p[i] = i;
std::sort(p.begin(), p.end(), IndexFunctor(keys));
// Reorder the keys and data in temporary arrays, it cannot be done in place.
std::vector<K> aux_keys(n);
std::vector<D> aux_data(n);
for (int i = 0; i < n; ++i) {
aux_keys[i] = keys[p[i]];
aux_data[i] = data[p[i]];
}
// Assign the ordered vectors by swapping, it should be faster.
keys.swap(aux_keys);
data.swap(aux_data);
}
#include <algorithm>
#include <iostream>
#include <string>
#include <vector>
//
// An input iterator that removes elements from the back of a container.
// Provided only because the standard library neglects one.
//
template<class Container>
class back_remove_iterator :
public std::iterator<std::input_iterator_tag, void, void, void, void> {
public:
back_remove_iterator() : container(0) {}
explicit back_remove_iterator(Container& c) : container(&c) {}
back_remove_iterator& operator=
(typename Container::const_reference value) { return *this; }
typename Container::value_type operator*() {
typename Container::value_type value(container->back());
container->pop_back();
return value;
} // operator*()
back_remove_iterator& operator++() { return *this; }
back_remove_iterator operator++(int) { return *this; }
Container* container;
}; // class back_remove_iterator
//
// Equivalence operator for back_remove_iterator. An iterator compares equal
// to the end iterator either if it is default-constructed or if its
// container is empty.
//
template<class Container>
bool operator==(const back_remove_iterator<Container>& a,
const back_remove_iterator<Container>& b) {
return !a.container ? !b.container || b.container->empty() :
!b.container ? !a.container || a.container->empty() :
a.container == b.container;
} // operator==()
//
// Inequivalence operator for back_remove_iterator.
//
template<class Container>
bool operator!=(const back_remove_iterator<Container>& a,
const back_remove_iterator<Container>& b) {
return !(a == b);
} // operator!=()
//
// A handy way to default-construct a back_remove_iterator.
//
template<class Container>
back_remove_iterator<Container> back_remover() {
return back_remove_iterator<Container>();
} // back_remover()
//
// A handy way to construct a back_remove_iterator.
//
template<class Container>
back_remove_iterator<Container> back_remover(Container& c) {
return back_remove_iterator<Container>(c);
} // back_remover()
//
// A comparison functor that sorts std::pairs by their first element.
//
template<class A, class B>
struct sort_pair_by_first {
bool operator()(const std::pair<A, B>& a, const std::pair<A, B>& b) {
return a.first < b.first;
} // operator()()
}; // struct sort_pair_by_first
//
// Performs a parallel sort of the ranges [keys_first, keys_last) and
// [values_first, values_last), preserving the ordering relation between
// values and keys. Sends key and value output to keys_out and values_out.
//
// This works by building a vector of std::pairs, sorting them by the key
// element, then returning the sorted pairs as two separate sequences. Note
// the use of std::move() for a vast performance improvement.
//
template<class A, class B, class I, class J, class K, class L>
void parallel_sort(I keys_first, I keys_last, J values_first, J values_last,
K keys_out, L values_out) {
typedef std::vector< std::pair<A, B> > Pairs;
Pairs sorted;
while (keys_first != keys_last)
sorted.push_back({std::move(*keys_first++), std::move(*values_first++)});
std::sort(sorted.begin(), sorted.end(), sort_pair_by_first<A, B>());
for (auto i = sorted.begin(); i != sorted.end(); ++i)
*keys_out++ = std::move(i->first),
*values_out++ = std::move(i->second);
} // parallel_sort()
int main(int argc, char** argv) {
//
// There is an ordering relation between keys and values,
// but the sets still need to be sorted. Sounds like a job for...
//
std::vector<int> keys{0, 3, 1, 2};
std::vector<std::string> values{"zero", "three", "one", "two"};
//
// parallel_sort! Unfortunately, the key and value types do need to
// be specified explicitly. This could be helped with a utility
// function that accepts back_remove_iterators.
//
parallel_sort<int, std::string>
(back_remover(keys), back_remover<std::vector<int>>(),
back_remover(values), back_remover<std::vector<std::string>>(),
std::back_inserter(keys), std::back_inserter(values));
//
// Just to prove that the mapping is preserved.
//
for (unsigned int i = 0; i < keys.size(); ++i)
std::cout << keys[i] << ": " << values[i] << '\n';
return 0;
} // main()
#ifndef HDR_PAIRED_ITERATOR
#define HDR_PAIRED_ITERATOR
#include <iterator>
/// pair_view mostly looks like a std::pair,
/// and can decay to a std::pair, but is really a pair of references
template <typename ItA, typename ItB>
struct pair_view {
typedef typename ItA::value_type first_type;
typedef typename ItB::value_type second_type;
typedef std::pair<first_type, second_type> pair_type;
pair_view() {}
pair_view(const ItA &a, const ItB &b):
first(*a), second(*b) {}
pair_view &operator=(const pair_view &x)
{ first = x.first; second = x.second; return *this; }
pair_view &operator=(const pair_type &x)
{ first = x.first; second = x.second; return *this; }
typename ItA::reference first;
typename ItB::reference second;
operator pair_type() const
{ return std::make_pair(first, second); }
friend bool operator==(const pair_view &a, const pair_view &b)
{ return (a.first == b.first) && (a.second == b.second); }
friend bool operator<(const pair_view &a, const pair_view &b)
{ return (a.first < b.first) || ((a.first == b.first) && (a.second < b.second)); }
friend bool operator!=(const pair_view &a, const pair_view &b)
{ return !(a == b); }
friend bool operator>(const pair_view &a, const pair_view &b)
{ return (b < a); }
friend bool operator<=(const pair_view &a, const pair_view &b)
{ return !(b < a); }
friend bool operator>=(const pair_view &a, const pair_view &b)
{ return !(a < b); }
friend bool operator==(const pair_view &a, const pair_type &b)
{ return (a.first == b.first) && (a.second == b.second); }
friend bool operator<(const pair_view &a, const pair_type &b)
{ return (a.first < b.first) || ((a.first == b.first) && (a.second < b.second)); }
friend bool operator!=(const pair_view &a, const pair_type &b)
{ return !(a == b); }
friend bool operator>(const pair_view &a, const pair_type &b)
{ return (b < a); }
friend bool operator<=(const pair_view &a, const pair_type &b)
{ return !(b < a); }
friend bool operator>=(const pair_view &a, const pair_type &b)
{ return !(a < b); }
friend bool operator==(const pair_type &a, const pair_type &b)
{ return (a.first == b.first) && (a.second == b.second); }
friend bool operator<(const pair_type &a, const pair_type &b)
{ return (a.first < b.first) || ((a.first == b.first) && (a.second < b.second)); }
friend bool operator!=(const pair_type &a, const pair_type &b)
{ return !(a == b); }
friend bool operator>(const pair_type &a, const pair_type &b)
{ return (b < a); }
friend bool operator<=(const pair_type &a, const pair_type &b)
{ return !(b < a); }
friend bool operator>=(const pair_type &a, const pair_type &b)
{ return !(a < b); }
};
template <typename ItA, typename ItB>
struct paired_iterator {
// --- standard iterator traits
typedef typename pair_view<ItA, ItB>::pair_type value_type;
typedef pair_view<ItA, ItB> reference;
typedef paired_iterator<ItA,ItB> pointer;
typedef typename std::iterator_traits<ItA>::difference_type difference_type;
typedef std::random_access_iterator_tag iterator_category;
// --- methods not required by the Random Access Iterator concept
paired_iterator(const ItA &a, const ItB &b):
a(a), b(b) {}
// --- iterator requirements
// default construction
paired_iterator() {}
// copy construction and assignment
paired_iterator(const paired_iterator &x):
a(x.a), b(x.b) {}
paired_iterator &operator=(const paired_iterator &x)
{ a = x.a; b = x.b; return *this; }
// pre- and post-increment
paired_iterator &operator++()
{ ++a; ++b; return *this; }
paired_iterator operator++(int)
{ paired_iterator tmp(*this); ++(*this); return tmp; }
// pre- and post-decrement
paired_iterator &operator--()
{ --a; --b; return *this; }
paired_iterator operator--(int)
{ paired_iterator tmp(*this); --(*this); return tmp; }
// arithmetic
paired_iterator &operator+=(const difference_type &n)
{ a += n; b += n; return *this; }
friend paired_iterator operator+(const paired_iterator &x, const difference_type &n)
{ return paired_iterator(x.a+n, x.b+n); }
friend paired_iterator operator+(const difference_type &n, const paired_iterator &x)
{ return paired_iterator(x.a+n, x.b+n); }
paired_iterator &operator-=(const difference_type &n)
{ a -= n; b -= n; return *this; }
friend paired_iterator operator-(const paired_iterator &x, const difference_type &n)
{ return paired_iterator(x.a-n, x.b-n); }
friend difference_type operator-(const paired_iterator &x, const paired_iterator &y)
{ return (x.a - y.a); }
// (in-)equality and ordering
friend bool operator==(const paired_iterator &x, const paired_iterator &y)
{ return (x.a == y.a) && (x.b == y.b); }
friend bool operator<(const paired_iterator &x, const paired_iterator &y)
{ return (x.a < y.a); }
// derived (in-)equality and ordering operators
friend bool operator!=(const paired_iterator &x, const paired_iterator &y)
{ return !(x == y); }
friend bool operator>(const paired_iterator &x, const paired_iterator &y)
{ return (y < x); }
friend bool operator<=(const paired_iterator &x, const paired_iterator &y)
{ return !(y < x); }
friend bool operator>=(const paired_iterator &x, const paired_iterator &y)
{ return !(x < y); }
// dereferencing and random access
reference operator*() const
{ return reference(a,b); }
reference operator[](const difference_type &n) const
{ return reference(a+n, b+n); }
private:
ItA a;
ItB b;
};
template <typename ItA, typename ItB>
paired_iterator<ItA, ItB> make_paired_iterator(const ItA &a, const ItB &b)
{ return paired_iterator<ItA, ItB>(a, b); }
#endif
#include <vector>
#include <algorithm>
#include <iostream>
template <typename ItA, typename ItB>
void print_kvs(const ItA &k0, const ItB &v0, const ItA &kn, const ItB &vn) {
ItA k(k0);
ItB v(v0);
while (k != kn || v != vn) {
if (k != kn && v != vn)
std::cout << "[" << *k << "] = " << *v << "\n";
else if (k != kn)
std::cout << "[" << *k << "]\n";
else if (v != vn)
std::cout << "[?] = " << *v << "\n";
if (k != kn) ++k;
if (v != vn) ++v;
}
std::cout << std::endl;
}
int main() {
std::vector<int> keys;
std::vector<std::string> data;
keys.push_back(0); data.push_back("zero");
keys.push_back(1); data.push_back("one");
keys.push_back(2); data.push_back("two");
keys.push_back(3); data.push_back("three");
keys.push_back(4); data.push_back("four");
keys.push_back(5); data.push_back("five");
keys.push_back(6); data.push_back("six");
keys.push_back(7); data.push_back("seven");
keys.push_back(8); data.push_back("eight");
keys.push_back(9); data.push_back("nine");
print_kvs(keys.begin(), data.begin(), keys.end(), data.end());
std::cout << "Shuffling\n";
std::random_shuffle(
make_paired_iterator(keys.begin(), data.begin()),
make_paired_iterator(keys.end(), data.end())
);
print_kvs(keys.begin(), data.begin(), keys.end(), data.end());
std::cout << "Sorting\n";
std::sort(
make_paired_iterator(keys.begin(), data.begin()),
make_paired_iterator(keys.end(), data.end())
);
print_kvs(keys.begin(), data.begin(), keys.end(), data.end());
std::cout << "Sort descending\n";
std::sort(
make_paired_iterator(keys.begin(), data.begin()),
make_paired_iterator(keys.end(), data.end()),
std::greater< std::pair<int,std::string> >()
);
print_kvs(keys.begin(), data.begin(), keys.end(), data.end());
return 0;
}
#include <iostream>
#include <iomanip>
#include <type_traits>
#include <utility>
#include <iterator>
#include <algorithm>
#include <numeric>
#include <deque>
#include <forward_list>
#include <vector>
#include <cstdlib>
#include <cassert>
template< typename pattern_iterator, typename target_iterator >
void
pattern_sort(pattern_iterator pbeg, pattern_iterator pend, target_iterator tbeg, target_iterator tend)
{
//assert(std::distance(pbeg, pend) == std::distance(tbeg, tend));
using pattern_traits = std::iterator_traits< pattern_iterator >;
using target_traits = std::iterator_traits< target_iterator >;
static_assert(std::is_base_of< std::forward_iterator_tag, typename pattern_traits::iterator_category >{});
static_assert(std::is_base_of< std::forward_iterator_tag, typename target_traits::iterator_category >{});
struct iterator_adaptor
{
iterator_adaptor(typename pattern_traits::reference pattern,
typename target_traits::reference target)
: p(&pattern)
, t(&target)
{ ; }
iterator_adaptor(iterator_adaptor &&)
: p(nullptr)
, t(nullptr)
{ ; }
void
operator = (iterator_adaptor && rhs) &
{
if (!!rhs.p) {
assert(!!rhs.t);
std::swap(p, rhs.p);
std::iter_swap(t, rhs.t);
}
}
bool
operator < (iterator_adaptor const & rhs) const
{
return (*p < *rhs.p);
}
private :
typename pattern_traits::pointer p;
typename target_traits::pointer t;
};
std::deque< iterator_adaptor > proxy_; // std::vector can be used instead
//proxy_.reserve(static_cast< std::size_t >(std::distance(pbeg, pend))); // it's (maybe) worth it if proxy_ is std::vector and if walking through whole [tbeg, tend) range is not too expensive operation (in case if target_iterator is worse then RandomAccessIterator)
auto t = tbeg;
auto p = pbeg;
while (p != pend) {
assert(t != tend);
proxy_.emplace_back(*p, *t);
++p;
++t;
}
std::sort(std::begin(proxy_), std::end(proxy_));
}
int
main()
{
std::forward_list< int > keys{5, 4, 3, 2, 1};
std::vector< std::size_t > indices(static_cast< std::size_t >(std::distance(std::cbegin(keys), std::cend(keys))));
std::iota(std::begin(indices), std::end(indices), std::size_t{0}); // indices now: 0, 1, 2, 3, 4
std::copy(std::cbegin(keys), std::cend(keys), std::ostream_iterator< int >(std::cout, " ")); std::cout << std::endl;
std::copy(std::cbegin(indices), std::cend(indices), std::ostream_iterator< std::size_t >(std::cout, " ")); std::cout << std::endl;
pattern_sort(std::cbegin(keys), std::cend(keys), std::begin(indices), std::end(indices)); std::cout << std::endl;
std::copy(std::cbegin(keys), std::cend(keys), std::ostream_iterator< int >(std::cout, " ")); std::cout << std::endl;
std::copy(std::cbegin(indices), std::cend(indices), std::ostream_iterator< std::size_t >(std::cout, " ")); std::cout << std::endl;
// now one can use indices to access keys and data to as ordered containers
return EXIT_SUCCESS;
}