library(gtools)
library(parallel)


############m?sik v?ltozat########
for (n in 5:8){
  print(date())
  print(n)
k=4
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)})
}

no_cores=detectCores()-1
cl=makeCluster(no_cores)
cl=makeCluster(no_cores)
clusterExport(cl, "a")
clusterExport(cl, "b")
clusterExport(cl, "s1")
clusterExport(cl, "s2")
clusterExport(cl, "s3")
clusterExport(cl, "s4")
clusterExport(cl, "s5")
clusterExport(cl, "s6")

print(date())
cmod=parSapply(cl,c(1:dim(a)[1]),function(i){s5(a[i,],b)})
print(date())
stopCluster(cl)

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(dmod,file = paste("matrixnegyzetes_k",k,"_n",n,sep=""))
save(Amod,file = paste("spectra_k",k,"_n",n,sep=""))
sertmod=Amod$values
sertE=round(sertmod)
hiba=max(abs(sertE-sertmod))
sertek=table(sertE)
sertek=sertek[c(length(sertek):1)]
write.table(sertek, paste("sertekEgesz_k_",k,"_n_",n,".res",sep=""),sep = "   ",col.names = TRUE,row.names = FALSE)

}
sertmod=Amod$values
sertE=round(sertmod)
hiba=max(abs(sertE-sertmod))