midSquareRand <- function(seed, len) {
   randvector <- NULL
   for(i in 1:len) {
      value <- seed * seed 
      randvector <- c(randvector, value)
      seed <- (value %/% 1000) %% 1000000
   }   
   return(randvector)
}

midSquareRand(675248, 100)

675248 * 675248
455959861504
   959861