Factorization of multivariate polynomials over finite fields