cp是一个2*n的矩阵，treecut是一个大矩阵groupfinding(m,n)是一个输入参数为正整数的函数，返回值是一个整数（用在后面的程序里定位向量内的元素）：
groupfinding <- function(m,n){
    return(2205-max(which(treecut[m,]==treecut[n,])))
}


现在需要将cp的每一行的两个元素作为groupfinding的两个输入参数进行运行，我知道这可以用for循环来实现：
for(n in 1:(length(cp)/2)){
    for(n in 1:(length(cp)/2)){
        m<-cp[n,1];
        n<-cp[n,2];
        Haddingup<-Haddingup+HC$height[groupfinding2(m,n)];
    }


其中Haddingup的初始值是0
但是for循环太慢了，有没有什么方法可以用apply来代替for循环？