fix: primefactors(): result and algorithm, bound enabled again
add: warning about non-polynomial coeffs in subst in manual