动态调整大小并在C++中填充矢量矢量

Question

我有一个矢量，它具有未知数量的行和3列。向量的构造如下：进行统计测试，如果它通过了阈值，则向量应该存储关于它的信息。 什么我做的是：

vector< vector <int> > validated_edge_list; 
validated_edge_list.resize(1); 
validated_edge_list.at(1).resize(3); 

for(int i = 0; i < e ; i++) 
{ 
    p = gsl_cdf_hypergeometric_P(edge_list[i][2], 
           k_vec[edge_list[i][1]], 
           M-k_vec[edge_list[i][1]], 
           N_vec[edge_list[i][0]]); // n2_matrix[i][j] = M-k_matrix[i][j] 

    if (p <= bonferroni_lvl) 
    { 
     validated_edge_list[c][0] = edge_list[i][0]; 
     validated_edge_list[c][1] = edge_list[i][1]; 
     validated_edge_list[c][2] = edge_list[i][2]; 
     c = c + 1; 
     validated_edge_list.resize(c+1); 
     validated_edge_list.at(c+1).resize(3); 
    } 
} 

正如你看到的我是一个新的原始每次手动添加。它给了我下面的错误：

terminate called after throwing an instance of 'std::out_of_range' 
what(): vector::_M_range_check: __n (which is 1) >= this->size() (which is 1) 
Aborted (core dumped) 

我认为我做错了什么，我也认为我应该使用的push_back选项，但我不知道怎么办。

我该如何解决这个问题？ （我是新的C++。）

Answer 1

就像你说的，你应该使用push_back。请勿使用resize。 push_back是为你的向量添加一个元素，并负责处理所有事情。然而，resize只会增加或减少向量的容量。

可能的解决办法：（我没有测试过，但是这应该给你的总体思路）

vector< vector <int> > validated_edge_list; 

for(int i = 0; i < e ; i++) 
{ 
    p = gsl_cdf_hypergeometric_P (edge_list[i][2],k_vec[edge_list[i][1]],M-k_vec[edge_list[i][1]],N_vec[edge_list[i][0]]); // n2_matrix[i][j] = M-k_matrix[i][j] 
    if (p <= bonferroni_lvl) 
    { 
     vector<int> single_edge_list = vector<int>(3); // Create a vector a 3 int's 
     single_edge_list[0] = edge_list[i][0] ; // Fill the vector. 
     single_edge_list[1] = edge_list[i][1] ; // Fill the vector. 
     single_edge_list[2] = edge_list[i][2] ; // Fill the vector. 

     validated_edge_list.push_back(single_edge_list); // Add it to validated_edge_list. 
     c++; // You don't really need this anymore 
    } 

} 

需要注意的是，因为里面validated_edge_list所有矢量有3个，你不的长度不需要使用矢量矢量，你可以使用structure（或类），你可以调用EdgeList。然而，这不是必要的。编辑：你可以找到更有效率和更好的方式做同样的事情（如YSC下面做的），但如果你是一个初学者在一个小程序上工作，而你并不介意降低效率，那么这应该很容易编程和足够好。

Answer 2

我强烈建议在内部线性化数据;向量矢量的开销很大（当外部向量重新定位它的内存时），读取和写入使用它们的代码是很痛苦的。

这是一个很小的想法：

template<class T, std::size_t width> 
struct Matrix 
{ 
    Matrix(std::size_t height) { _data.reserve(width*height); } 
    const T& operator()(std::size_t i, std::size_t j) const { return _data[i*width+j]; } 
      T& operator()(std::size_t i, std::size_t j)  { return _data[i*width+j]; } 

private: 
    std::vector<T> _data; 
}; 

Usage example on coliru