faster ``_choose_factor``

For instance in ``minpoly(sum(sin(i*pi/13) for i in range(13)))``  ``_choose_factor`` is more than 2x faster