在我的第一个问题(详细的实验研究)中:Appropriate container for the fast insertion and lookup of n-dimensional real vectors (initial benchmarking provided)我使用未排序的集来管理我的初始N维浮点数组(可能设计不佳的Hash函数),因此获得了非常奇怪的行为:

#include <iostream>
#include <chrono>
#include <random>
#include <array>
#include <unordered_set>

const int N = 3;  // Dimensionality of the arrays

std::array<double, N> getRandomArray() {
  // Engines and distributions retain state, thus defined as static
  static std::default_random_engine e;                    // engine
  static std::uniform_real_distribution<double> d(0, 1);  // distribution
  std::array<double, N> ret;
  for (size_t i = 0; i < N; ++i) {
    ret[i] = d(e);
  }
  return ret;
}

// Return Squared Euclidean Distance
template <typename InputIt1, typename InputIt2>
double EuclideanDistance2(InputIt1 beg1, InputIt1 end1, InputIt2 beg2) {
  double val = 0.0;
  while (beg1 != end1) {
    double dist = (*beg1++) - (*beg2++);
    val += dist*dist;
  }
  return val;
}

struct ArrayHash {  // Hash Function
  std::size_t operator() (const std::array<double, N>& arr) const {
    std::size_t ret = 0;
    for (const double elem : arr) {
      ret += std::hash<double>()(elem);
    }
    return ret;
  }
};

struct ArrayEqual {  // Equivalence Criterion
  bool operator() (const std::array<double, N>& arr1,
                          const std::array<double, N>& arr2) const {
    return EuclideanDistance2(arr1.begin(), arr1.end(), arr2.begin()) < tol*tol;
  }
 private:
  static constexpr double tol = 1e-6;  // Comparison tolerance
};


int main() {
  // create a unordered set of double arrays (usda)
  std::unordered_set<std::array<double, N>, ArrayHash, ArrayEqual> usda;
  // record start time
  auto start = std::chrono::steady_clock::now();
  // Generate and insert one hundred thousands new double arrays
  for (size_t i = 0; i < 100000; ++i) {
    // Get a new random double array (da)
    std::array<double, N> da = getRandomArray();
    usda.insert(da);
  }
  // record finish time
  auto end = std::chrono::steady_clock::now();
  std::chrono::duration<double> diff = end - start;
  std::cout << "Time to generate and insert unique elements into UNORD. SET: "
            << diff.count() << " s\n";
  std::cout << "unord. set size() = " << usda.size() << std::endl;
  return 0;
}

两件最奇怪的事情是:
  • 即使容差(1e-1)宽松也没有任何优化标志的运行实验,几乎所有随机 vector (实现为N维数组)都被识别为唯一。我没有使用vectorssets观察到这一点。
  • 启用-O3优化标志时,唯一元素的数量与未进行优化的数量明显不同,这可以肯定地说我的方法有问题。

  • 编辑:解决了第二个问题,同时考虑到@WhozCraig备注。

    因此,我的问题是:这是奇怪的行为,因为我的哈希函数设计错误吗?如果是这样,您能否建议如何为我的案例制作更好的哈希函数?

    最佳答案

    您的程序表现出未定义的行为(未初始化的std::size_t ret变量除外)。

    首先,ArrayEqual不产生等价关系。它不是可传递的-存在三个数组abc,使得ab足够“接近”,而bc足够“接近”,但是ac不够“足够”。

    其次,对于ArrayHash声明相等的两个数组,ArrayEqual可能不会返回相同的哈希值。

    这两个都是std::unordered_set模板参数的先决条件。

    07-28 09:09