如何用C++20基于Rcpp::NumericVector生成重复元素的向量？

针对你的需求，这里提供几种比嵌套循环更简洁或高效的实现方案，均适配C++20标准：

方案1：直接使用Rcpp内置的rep函数（最推荐） #

Rcpp封装了R语言的rep函数，能直接实现向量重复，底层是R的高效实现，代码极简且性能优异：

// 若x是Rcpp::NumericVector
Rcpp::NumericVector y = Rcpp::rep(x, n);

// 若x是std::vector<double>，先转换为Rcpp::NumericVector再重复
Rcpp::NumericVector x_rcpp(x.begin(), x.end());
Rcpp::NumericVector y = Rcpp::rep(x_rcpp, n);

这个方案完全不需要手动处理循环和内存分配，可读性和性能都是最优的。

方案2：用C++标准库std::copy简化循环 #

如果倾向于纯C++标准库实现，可以用std::copy批量复制整个向量，替代内层循环，代码更简洁，且编译器会对std::copy做优化，性能和嵌套循环相当甚至更好：

// 处理Rcpp::NumericVector的情况
const int x_size = x.size();
Rcpp::NumericVector y(x_size * n);
auto dest_it = y.begin();
const auto x_begin = x.begin();
const auto x_end = x.end();

for (int i = 0; i < n; ++i) {
    dest_it = std::copy(x_begin, x_end, dest_it);
}

// 处理std::vector<double>的情况
const int x_size = x.size();
Rcpp::NumericVector y(x_size * n);
auto dest_it = y.begin();
const auto x_begin = x.cbegin();
const auto x_end = x.cend();

for (int i = 0; i < n; ++i) {
    dest_it = std::copy(x_begin, x_end, dest_it);
}

方案3：底层内存复制（极致性能） #

对于double这类POD类型，可以用std::memcpy直接批量复制内存块，在某些场景下能获得极致性能：

// 处理Rcpp::NumericVector的情况
const int x_size = x.size();
const size_t bytes = x_size * sizeof(double);
Rcpp::NumericVector y(x_size * n);

double* dest_ptr = y.begin();
const double* x_ptr = x.begin();

for (int i = 0; i < n; ++i) {
    std::memcpy(dest_ptr, x_ptr, bytes);
    dest_ptr += x_size;
}

关于动态插入的方案（不推荐） #

你提到的动态调整大小并插入自身副本的方式，虽然可行但性能较差——每次insert可能触发内存重新分配。如果一定要用这种方式，代码如下：

// 处理Rcpp::NumericVector
Rcpp::NumericVector y = x;
for (int i = 1; i < n; ++i) {
    y.insert(y.end(), x.begin(), x.end());
}

但该方案内存分配次数多，整体性能不如预分配内存的方案，仅适合小向量场景。

内容的提问来源于stack exchange，提问作者user