For comparision with Haskell, we can define
R for any given list of factors by folding the merge operator through the list of lists of multiples:
r = foldr merge [] . map (\n -> map (*n) [1..])
We use a duplicate-eating merge:
merge [] ys = ys
merge xs [] = xs
merge (x:xs) (y:ys)
| x < y = x : merge xs (y:ys)
| y < x = y : merge (x:xs) ys
| otherwise = x : merge xs ys
Sample output:
> take 60 (r [2,3,5])
[2,3,4,5,6,8,9,10,12,14,15,16,18,20,21,22,24,25,26,27,28,30,32,33,34,3
+5,36,38,39,40,42,44,45,46,48,50,51,52,54,55,56,57,58,60,62,63,64,65,6
+6,68,69,70,72,74,75,76,78,80,81,82]