Solving systems of polynomial congruences modulo a large prime

We consider the following polynomial congruences problem: given a prime p, and a set of polynomials f₁,...,f_m∈ F_p[x₁,...,x_n] of total degree at most d, solve the system f₁=...=f_m=0 for solution(s) in F_pⁿ. We give a randomized algorithm for the decision version of this problem. When the system has F_p-rational solutions our algorithm finds one of them as well as an approximation of the total number of such solutions. For a fixed number of variables, the algorithm runs in random polynomial time with parallel complexity poly-logarithmic in d, m and p, using a polynomial number of processors. As an essential step of the algorithm, we also formulate an algebraic homotopy method for extracting components of all dimensions of an algebraic set. The method is efficiently parallelizable