如何用C++在matlab中实现unique的行为？

Question

在Matlab中，[C，ia，ic]=unique(A)函数返回一个排序数组C，其中A中的重复项被删除，ia和ic数组包含索引，使得C=A(ia)和A=C(ic)。例如:

A = [1, 2, 3, -1, -2, 0, 0, 0, 4]; 
[C, ia, ic] = unique(A); 
C = [-2, -1, 0, 1, 2, 3, 4]; 
% The indices are incremented by 1 to accomodate the C++ convection.
ia = [4, 3, 7, 0, 1, 2, 8]; % size(ia) = size(C) 
ic = [3, 4, 5, 1, 0, 2, 2, 2, 6]; % size(ic) = size(A) 

下面是我的实现:

  std::vector<int> A = {1, 2, 3, -1, -2, 0, 0, 0, 4};
  std::set<int> C(A.begin(), A.end());

  std::vector<int> ic;
  ic.reserve(A.size());

  std::transform(A.begin(), A.end(), std::back_inserter(ic),
                 [&](int x) { return (std::distance(C.begin(), C.find(x))); });

我现在可以正确地得到阵列C和ic，但是我不知道如何得到阵列IA。

有人能帮我一下吗？

Answer 1

另一种方法在语法上与MATLAB类似:

#include <algorithm>
#include <vector>

std::vector<int> A = {1, 2, 3, -1, -2, 0, 0, 0, 4};
std::sort(A.begin(), A.end());
auto last = std::unique(A.begin(), A.end());
A.erase(last, A.end()); 

编译器资源管理器链接

Answer 2

答案与您计算ic的方式非常相似。只需在transform调用中交换a和c:

std::vector<int> ia;
ia.reserve(C.size());

std::transform(C.begin(), C.end(), 
               std::back_inserter(ia),
               [&](int x) { return std::distance(A.begin(), 
                                     std::find(A.begin(), A.end(), x)); 
});

您还需要使用std::find，因为std::vector不提供.find()成员函数。