Rcpp生成随机数异常：调用无关函数后输出固定值

你遇到的问题核心在于R的随机数状态管理机制：

• R的全局随机数生成器状态由R runtime管控，当通过Rcpp的Function对象调用R导出的函数（比如示例中的f()）时，R会在进入函数前保存当前随机数状态，函数执行完毕后自动恢复该状态。
• 这就导致每次调用g()后，随机数状态都会回到调用g()之前的节点，后续调用R::runif()时，每次都会生成相同的随机数，甚至和第一个循环的起始值一致——因为第一个循环结束后的状态，在每次g()调用后都被强制恢复了。

方法1：改用Rcpp的随机数接口 #

Rcpp提供了自身的随机数函数（比如Rcpp::runif()），这些函数内部会正确处理与R随机数状态的同步，避免R函数调用导致的状态回滚问题。修改后的代码示例：

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
void f() {}

// [[Rcpp::export]]
void test(Function g) { 
  for(int j=0;j<5;++j) {
    Rcout << Rcpp::runif(1, 0, 1)[0] << " ";
  }
  Rcout << "\n";
  for(int j=0;j<5;++j) {
    g();
    Rcout << Rcpp::runif(1, 0, 1)[0] << " ";
  }
}

方法2：绕过R函数调用的状态恢复 #

如果必须调用R函数且要保留R::runif()，可以直接在C++中实现g()的逻辑（比如示例中f()是空函数，完全可以跳过R函数调用），避免触发R的状态保存与恢复机制。

方法3：手动管理随机数状态（不推荐，仅作备选） #

你可以在调用g()前保存当前随机数状态，调用后手动更新状态，强制推进随机数生成器的进度：

#include <Rcpp.h>
#include <R_ext/Random.h>
using namespace Rcpp;

// [[Rcpp::export]]
void f() {}

// [[Rcpp::export]]
void test(Function g) { 
  for(int j=0;j<5;++j) {
    Rcout << R::runif(0,1) << " ";
  }
  Rcout << "\n";
  
  // 保存当前随机数状态
  GetRNGstate();
  unsigned int seed1 = RNG_seed[0];
  int seed2 = RNG_seed[1];
  PutRNGstate();
  
  for(int j=0;j<5;++j) {
    g();
    // 手动恢复并推进状态
    GetRNGstate();
    RNG_seed[0] = seed1;
    RNG_seed[1] = seed2;
    PutRNGstate();
    // 生成一个随机数推进状态
    R::runif(0,1);
    // 更新保存的状态
    GetRNGstate();
    seed1 = RNG_seed[0];
    seed2 = RNG_seed[1];
    PutRNGstate();
    // 输出目标随机数
    Rcout << R::runif(0,1) << " ";
  }
}

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