Rcpp 问题。坐等大牛

cyc

求一组向量set的前topN最小值以及对应的下标。

利用sourceCpp将下面的函数导入到R里面，多次运行sorttopNC函数，R会自动关闭。没有任何的error出现。[s:11]。大神求解额~ （注：maxC函数测试了，应该没有任何问题）

x=sample(1:100,20)

sorttopNC(x,5)

#************程序如下：

<br />
# include &lt Rcpp.h &gt<br />
using namespace Rcpp;</p>
<p>// [[Rcpp::export]]<br />
int maxC(NumericMatrix dset){<br />
int n=dset.nrow();<br />
int kindex=0;<br />
double maxv=dset(0,0);<br />
for(int i=1;i<n;i++){<br />
if(dset(i,0)>maxv){<br />
kindex=i;<br />
maxv=dset(i,0);<br />
}<br />
}<br />
return(kindex);<br />
}</p>
<p>// [[Rcpp::export]]<br />
NumericMatrix sorttopNC(NumericVector set, int topN)<br />
{<br />
   int setlen=set.size();<br />
  NumericMatrix result(topN,2);<br />
  for( int i=0;i<topN;i++){<br />
   result(i,0)=set[i];<br />
   result(i,1)=i; }<br />
int maxkey=maxC(result);<br />
double maxvalue=result(maxkey,0);<br />
int k=topN;</p>
<p>while(k<setlen){<br />
if(set[k]<maxvalue){<br />
result(maxkey,0)=set[k];<br />
result(maxkey,1)=k;<br />
maxkey=maxC(result);<br />
maxvalue=result(maxkey,0);<br />
}<br />
k++;<br />
          }<br />
for(int i=0;i<setlen;i++){<br />
result(i,1)+=1;<br />
}<br />
return(result);<br />
}<br />

</p>

superdesolator

回复第1楼的 cyc：问下int n=dset.nrow();和int setlen=set.size();

这2个返回函数是Rccp中自带的功能么？

superdesolator

回复第1楼的 cyc：略显无语,我的错误在哪?

#include<Rcpp.h><br />
using namespace Rcpp;<br />
// [[Rcpp::export]]<br />
NumericVector sort(NumericVector a)<br />
{<br />
int setlen=a.size();<br />
for(int i=0;i<setlen-1;i++)<br />
{<br />
for (int j=i;j<setlen-1-i;j++)<br />
{<br />
 int temp=0;<br />
 if(a[j]>a[j+1])<br />
{<br />
 temp=a[j];<br />
a[j]=a[j+1];<br />
a[j+1]=temp;<br />
}<br />
}}<br />
return(a);<br />
}

怎么感觉不行啊.看样子还得每次都比较所有的啊。。。。
</p>

cyc

回复第2楼的 superdesolator：int n=dset.nrow();和int setlen=set.size();分别是取矩阵的行数和向量的长度。是Rcpp里面的函数，应该Rcpp sugar的一种。

cyc

回复第3楼的 superdesolator：你这种算法是冒泡排序，时间复杂度是O(n^2),效率太低。远不如快速排序O(n*logn)等。而且你的代码还没有记录对应的下标。

我的问题是只需要找到向量的前topN个最小值和对应的下标。并不需要关心其他值和他们之间的顺序。所以并不需要对所有的数值进行排序。

我的这种只需对set进行一次循环。result记录topN个最小值，maxC函数查找result(第一列是值，第二列是下标)第一列的最大值对应的下标。对于set的第i个元素来说，如果比maxvalue大，则进行一次替换操作。时间复杂度最坏的情况是topN*n,最好的情况是n

cyc

回复第3楼的 superdesolator：

正确的冒泡排序是这个

cppFunction('

NumericVector sortC(NumericVector a)

{

int setlen=a.size();

for(int i=0;i<setlen;i++)

{

for (int j=1;j<setlen-i;j++)

{

double temp=0;

if(a[j]<a[j-1])

{

temp=a[j];

a[j]=a[j-1];

a[j-1]=temp;

}

}}

return(a);

}

')

常用的排序如下，R里面的sort就是快速排序~

http://www.caogenit.com/caogenxueyuan/yingyongfangxiang/rengongzhineng/1724.html

nan.xiao

珍爱绳命，请使用 stl 。。

superdesolator

回复第6楼的 cyc：你的代码好像没有完全显示...嗯,我懂你的那段代码什么意思了,我对算法快慢只有个大致的了解,却不知道原来冒泡不如你这种块。你的我运行了3次就R关闭了，我不知什么原因，但感觉到和内存处理释放变量机制有关吧，我纯不懂啊，因为我设计的时候，错把输出成了NumericVector，但是函数返回的是NumericMatrix,就出现了R自动关闭。。不知道你的是不是运行次数多了，内存混乱了。。。。当然，再次重申，我纯不懂。。。[s:11]

bjt

冒泡算法：http://www.numbertheory.nl/2013/05/14/much-more-efficient-bubble-sort-in-r-using-the-rcpp-and-inline-packages/

bjt

回复第7楼的肖楠：这个是什么意思？时序分解？

nan.xiao

回复第10楼的刘思喆：

C++ 标准模板库

http://gallery.rcpp.org/articles/sorting/

cyc

顶一下，还没有解决额[s:12]

lcf

我就改了一个地方就没有错了~

#include <Rcpp.h> 这个地方除了中间能有空格以外，其他地方都不能有~~

如果还不行，估计Rcpp没安装好~

cyc

回复第13楼的 lcf：不是这个的问题吧。问题是运行一次正常。运行多次R console自动关闭~

cyc

ding