1
有問題試圖平均係數時間(pr)的值來分解乘以2的素數。 目前有一個迭代工作,但不確定如何讓它與n次迭代一起工作。如何在R中創建平均系統時間?
require("gmp")
prime<-function(a) {
pr <- vector(mode = "numeric", length = a)
sz <- vector(mode = "numeric", length = a)
for(i in 1:a){
y<-urand.bigz(nb = 2, size = i)
p1<-nextprime(as.bigz(y[1]))
p2<-nextprime(as.bigz(y[2]))
p<-p1*p2
pr[i] <-system.time(factorize(p))[1]
sz[i] <-sizeinbase(p)
}
return(cbind(sz,pr))
}
times <- prime(48)
plot(times,xlab="Size of Psuedoprime ",ylab="Time taken to factorize", main="Size of Vs Time taken to factorize of a given Pseudoprime")
任何幫助,將不勝感激
非常感謝,一切似乎都奏效。 – Stuck