斐波纳契系列使用Rcpp的意外结果

时间:2013-09-18 16:25:23

标签: r rcpp

我刚刚开始使用Rcpp所以很抱歉,如果我错过了一个简单的步骤或类似的东西......我试过?sourceCpp

library(Rcpp)
sourceCpp(code='
  #include <Rcpp.h>

  // [[Rcpp::export]]
  int fibonacci(const int x) {
    if (x == 0) return(0);
    if (x == 1) return(1);
    return (fibonacci(x - 1)) + fibonacci(x - 2);
  }'
)

最多fibonacci(46)一切都很好,但后来我得到了:

> fibonacci(47)
[1] -1323752223
> fibonacci(48)
[1] 512559680
> fibonacci(49)
[1] -811192543
> fibonacci(50)
[1] -298632863

根据这个page,上面应该是:

47 : 2971215073
48 : 4807526976
49 : 7778742049
50 : 12586269025

你得到相同的结果吗?

1 个答案:

答案 0 :(得分:5)

对于有符号整数,你超过了maximum limit(技术上,我认为这是long int)。请改用double ...

library(Rcpp)
sourceCpp(code='
  #include <Rcpp.h>

  // [[Rcpp::export]]
  double fibonacci(const double x) {
    if (x == 0) return(0);
    if (x == 1) return(1);
    return (fibonacci(x - 1)) + fibonacci(x - 2);
  }'
)

fibonacci(47)
#[1] 2971215073