n=11 p=c(5,2,11,9,8,1,3,10,4,7,6) kappa=function(p,N){ kappa=rep(0,N) for (u in 1:N){ v=p[u] kappa[v]=1 w=1 if(v