library(gtools)
n=8
k=5
a=permutations(n,n)
b=permutations(n,k)
seged=function(aa,bb){
  prod((which(aa==bb[c(1:(length(bb)-1))])<which(aa==bb[c(1:(length(bb)-1))+1])))
}
seged2=function(aaa,b){
  seged2=sapply(c(1:dim(b)[1]), function(i){seged(aaa,b[i,])})
}
print(date())
c=sapply(c(1:dim(a)[1]), function(i){seged2(a[i,],b)})
print(date())

print(date())
d=c%*%t(c)
print(date())

print(date())
A=eigen(d,symmetric = TRUE)
print(date())

sert=A$values
sertkl=kmeans(sert,50)
print(date())
library(RSpectra)
sert2=eigs(as(d,"dgCMatrix"),20,which = "LR")

############másik változat########
n=8
k=5
a=permutations(n,n)
b=permutations(n,k)

#az "x" vektorból megnézzük, hogy az "y" melyik helyen van
s1=function(x,y){
  s1=which(x==y)
}

#az "x" vektorból megnézzük, hogy az "yvektor" elemei mely helyeken vannak
s2=function(x,yvektor){
  s2=sapply(c(1:length(yvektor)),function(i){s1(x,yvektor[i])})
}

#Megnézzük, hogy a "z" vektorban az elemek növekvő sorrendben vannak-e
s3=function(z){
  s3=prod(sapply(c(1:(length(z)-1)),function(i){(z[i]<z[i+1])}))
}

#Megnézzük, hogy az "x" vektorban az "yvektor" elemei ugyanabban a sorrendben vannak-e
s4=function(x,yvektor){
  s4=s3(s2(x,yvektor))
}

#Megnézzük, hogy az "x" vektorban a variációk elemei elemei ugyanabban a sorrendben vannak-e
s5=function(x,b){
  s5=sapply(c(1:dim(b)[1]),function(i){s4(x,b[i,])})
}

#Megnézzük, hogy a permutációkban a variációk elemei elemei ugyanabban a sorrendben vannak-e
s6=function(a,b){
  s6=sapply(c(1:dim(a)[1]),function(i){s5(a[i,],b)})
}


print(date())
cmod=s6(a,b)
print(date())

print(date())
dmod=cmod%*%t(cmod)
print(date())

print(date())
Amod=eigen(dmod,symmetric = TRUE)
print(date())

save(cmod,file = paste("matrix_k",k,"_n",n,sep=""))
save(Amod,file = paste("spectra_k",k,"_n",n,sep=""))

sertmod=Amod$values

print(date())
kidob=sapply(c(1:10),function(i){s5(a[i,],b)})
print(date())