我的小std::独特的例子不起作用

My small std::unique example not working

本文关键字：不起作用 std 我的      更新时间：2023-10-16

我正在尝试使用cxx-11的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;    
}

但是当我运行这个代码时（链接到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
---

唯一数组仍然包含一个重复的数组（这是错误的）！

我做错了什么？

让我们尝试更精确的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;});

怎么样：

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

实时演示链接

相关文章：