Unordered_map（C ++）中的“错误：分配只读位置”

Question

我有一个笨拙的哈希表（特别是unordered_map），其中包含int键和vector< vector< int >>数据。 我定期需要更新此int的二维向量中的元素。 没有内在的原因我不应该对吧？ 我改用了一种较新的g ++编译器，抱怨在下面指定的行上分配了只读位置。

typedef std::tr1::unordered_map< int, vector< vector< int > > > pimap;

vector< Strain * > liveStrains;
pimap phenotypeIs;
int NUM_DEMES = 3;

...
vector< Strain * >::const_iterator lsItr;
for ( lsItr = liveStrains.begin(); lsItr != liveStrains.end(); ++lsItr ) {
 int thisP = (*lsItr)->getPhenotype();
 pimap::iterator piItr = phenotypeIs.begin();
 piItr = phenotypeIs.find( thisP );
 if ( piItr != phenotypeIs.end() ) {
   for ( int d = 0; d < NUM_DEMES; d++ ) {
      ( piItr -> second )[ thisStep ].at( d ) = (*lsItr)->getI( d );  // error here
   }
 }
}

我是C ++的新手，所以没有什么太明显的。 感谢您的任何帮助。

---

遵循蒂姆的建议

我已将以下代码的相关部分替换为以下内容：

  pimap::iterator piItr = phenotypeIs.find( thisP );
  if ( piItr != phenotypeIs.end() ) {
    for ( int d = 0; d < NUM_DEMES; d++ ) {
      vector< vector< int > > & thisVec2 = piItr->second;
      vector<int> & thisVec = thisVec2.at( thisStep );
      int & ii = thisVec.at( d );
      ii = (*lsItr)->getI( d );
      // ( piItr -> second )[ thisStep ].at( d ) = (*lsItr)->getI( d ); // error was here
    }

该代码可以编译而不会出现错误，并且可以正常运行。 像蒂姆一样，我仍然不太明白为什么该修复程序有效。 该错误先前在gcc版本4.1.2 20080704（Red Hat 4.1.2-44）中出现，但在gcc版本4.0.1（Apple Inc.内部版本5465）中没有出现。 如果我没有紧迫的期限，我将尝试更仔细地分析错误！

Answer 1

你确定真的有thisStep + 1中的每一级的矢量元素和NUM_DEMES元素在每一个二级矢量？

如果我正确阅读，实际上并没有分配给地图迭代器，因此我怀疑错误在于矢量访问。

将最后一个语句分解为多个语句可能会有所帮助，以便对每个语句仅执行一项操作以缩小问题所在。 例如，

Strain* strain = *lsItr;
vector<vector<int> >& vv = piItr->second;
vector<int>& v = vv[thisStep];
int& i = v.at(d);     // <-- My bet is that the error occurs here or the prev. line
i = strain->getI( d );

顺便说一句， piItr = phenotypeIs.begin(); 在这里没有影响，可能很简单：

pimap::iterator piItr = phenotypeIs.find( thisP );

Answer 2

( piItr -> second )[ thisStep ].at( d )

at（）将迭代器返回到内部向量中，而不访问该值。 你想要的是

 *(( piItr -> second )[ thisStep ].at( d ))