prods = [a*b | a <- numbers, b <- numbers] where numbers = [1..100]
stats = map (\d -> (d, length $ filter (\p -> p `mod` d == 0) prods)) [1..10]
main = print stats