Question

我使用multimap来存储编码原子反应的对象。 multimap看起来像这样：

std::multimap<ReactionElement, ReactionElement> reaction_map;

键是反应物，值是产品。然后我遇到了一种情况，我找到了两个反应物原子，我可以在地图中查找这些反应物可能形成的可能产物。

ReactionElements类的基本内容如下所示：

class ReactionElement {

  friend bool operator==(const ReactionElement& lhs, const ReactionElement& rhs);
  friend bool operator<(const ReactionElement& lhs, const ReactionElement& rhs);

public:
  // Some methods here ...

private:
  // The Atom class tracks the element
  Atom atom_a;
  Atom atom_b;

  // SiteSpecies and NeighborSpecies classes track the reaction geometry
  SiteSpecies site_species_a;
  NeighborSpecies neighbor_species_b;

  // int members track reaction energetics
  int e_initial, e_transition, e_final;

  double reac_distance;

}; // ReactionElement

我从文件中读取了一堆ReactionElement-ReactionElement对，并将每个对插入到多图中。

问题在于：当我去检索它们时，只有一些反应是可以访问的。对于某些反应，我可以调用reaction_map.find（反应物）并且什么也得不到。但是，我可以遍历多图，我看到我插入的所有对。我用它作为基本诊断：

for(multimap<ReactionElement, ReactionElement>::iterator it = reaction_map.begin();
    it != reaction_map.end(); ++it)
  std::cout << reaction_map.count(it->first) << '\n';

这打印了一些，两个和零。这怎么可能？

有什么想法吗？

（编辑）比较运算符详述如下。我相当确定每个成员类都是有序的。我会证实。

// This directly compares each member (using a tolerance of .01 for reac_distance).
bool operator==(const ReactionElement& lhs, const ReactionElement& rhs) {
  return (lhs.atom_a == rhs.atom_a and lhs.atom_b == rhs.atom_b and
      lhs.site_species_a == rhs.site_species_a and
      lhs.neighbor_species_b == rhs.neighbor_species_b and
      lhs.e_initial == rhs.e_initial and lhs.e_transition == rhs.e_transition and
      lhs.e_final == rhs.e_final and
      fabs(lhs.hop_distance-rhs.hop_distance) <= 0.01);
}

// < orders on members in the order they appear in the class definition.
bool operator<(const ReactionElement& lhs, const ReactionElement& rhs) {
  if (lhs.atom_a < rhs.atom_a)
    return true;
  else if (lhs.atom_a == rhs.atom_a and lhs.atom_b < rhs.atom_b) 
    return true;
  else if (lhs.atom_a == rhs.atom_a and lhs.atom_b == rhs.atom_b and 
           lhs.site_species_a < rhs.site_species_a)
    return true;
  // etc for the remaining members.

  return false;
}

Answer 1

迭代与count / find之间的区别在于count和find必须比较元素（迭代不会）。默认情况下，std::multimap使用std::less进行比较。这相当于<。 std::multimap要求此关系为strict weak ordering。这意味着特别是。以下规则：

!(x < x)
如果(x < y)，则!(y < x)

确保您的operator<实施符合这些规则。如果没有，请更改它，或为Comparator提供自定义std::multimap。

C ++多地图查找失败

1 个答案: