Question

我试图使用cxx-11＆＃39; std::unique()来查找 array中的唯一元素：

#include <iostream>
#include <algorithm>
#include <vector>
#include <typeinfo>


int main(){
    const int n=11;
    double x[n],a3[n],a1[n];

    x[0]=-0.717778;
    x[1]=-0.496843;
    x[2]=-0.429063;
    x[3]=-0.3596;
    x[4]=-0.205607;
    x[5]=0.0730536;
    x[6]=0.138018;
    x[7]=0.585526;
    x[8]=2.40104;
    x[9]=3.75268;
    x[10]=4.55704;

    a3[0]=0.790832;
    a3[1]=0.569896;
    a3[2]=0.502116;
    a3[3]=0.432653;
    a3[4]=0.343625;
    a3[5]=0.512472;
    a3[6]=0.56708;
    a3[7]=1.01459;
    a3[8]=2.32799;
    a3[9]=3.67962;
    a3[10]=4.48398;

    std::cout.precision(10);
    std::copy(a3,a3+n,a1);
    for(int i=0;i<n;i++)            a1[i]+=x[i];            
    std::sort(a1,a1+n);
    for(int i=0;i<n;i++)            std::cout << a1[i] << std::endl;
    std::cout << "---" << std::endl;
    int n_1=std::unique(a1,a1+n)-a1;
    std::cout << "length of unique subvector " << n_1 << std::endl;
    std::cout << "---" << std::endl;
    for(int i=0;i<n_1;i++)         std::cout << a1[i] << std::endl;
    std::cout << "---" << std::endl;    
}

但是当我运行此code (link to coliru)时它返回：

original array
0.073053
0.073053
0.073053
0.073054
0.138018
0.5855256
0.705098
1.600116
4.72903
7.4323
9.04102
---
length of unique subarray 10
---
unique array
0.073053
0.073053
0.073054
0.138018
0.5855256
0.705098
1.600116
4.72903
7.4323
9.04102
---

唯一数组仍包含副本（因此错误）！

我做错了什么？

Answer 1

让我们更精确地尝试std::cout.precision(20)：

0.073052999999999979064
0.073053000000000034575
0.073053999999999952308
0.13801800000000000179
0.58552559999999997942
0.70509800000000000253
1.6001160000000000938
4.7290299999999998448
7.4322999999999996845
9.0410199999999996123

由于大多数小数部分无法用二进制浮点格式精确表示，因此稍微不同的舍入误差会导致稍微不同的结果。

通常，您不能指望不同浮点计算的结果完全相等，即使应用于数学实数的相应计算也是如此。

您可以测试“几乎相等”，仔细选择适合您的数字域的容差。 unique允许您指定自己的谓词，而不是简单的相等测试：

std::unique(a1,a1+n,[](double x, double y){return std::abs(x-y) < tolerance;});

Answer 2

怎么样：

int n_1 = std::unique(a1,a1+n,
            [](float a, float b)
            {
                return std::fabs(a-b) < 10e-9;
            }
          ) - a1;

Live demo link

我的小std :: unique例子不起作用

2 个答案: