Add factorize()