First draft of a function that lazily returns all factorizations