Add prime_factorise_range(n)