我刚刚开始使用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);
}'
)
Up to fibonacci(46)
一切都很好,但后来我得到:
> fibonacci(47)
[1] -1323752223
> fibonacci(48)
[1] 512559680
> fibonacci(49)
[1] -811192543
> fibonacci(50)
[1] -298632863
根据这个page http://www.maths.surrey.ac.uk/hosted-sites/R.Knott/Fibonacci/fibtable.html上面应该是:
47 : 2971215073
48 : 4807526976
49 : 7778742049
50 : 12586269025
你得到同样的结果吗?
您超出了最大限制 http://www.cplusplus.com/reference/climits/对于有符号整数(从技术上讲,这将是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
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)