gistfile1.hs

import Control.Monad
import System.Random
import Data.List

size :: Int
size = 3

steps :: Int
steps = 10000

main :: IO ()
main = do
  results <- replicateM steps (simulate size)
  print $ avg results

simulate :: Int -> IO Int
simulate k = do
  gen <- newStdGen
  let ns = randomRs (1,k) gen
  return $ length $ takeWhile ((<k) . length . nub) $ inits ns

avg :: [Int] -> Float
avg xs = fromIntegral (sum xs) / genericLength xs