A subdivision approach to the solution of polynomial constraints over finite domains using the modified Bernstein form