williamsp plus 1 factorization method