  Let $\mathbf{K}$ be a field of characteristic zero with
$\overline{\mathbf{K}}$ its algebraic closure. Given a sequence of polynomials
$\mathbf{g} = (g_1, \ldots, g_s) \in \mathbf{K}[x_1, \ldots , x_n]^s$ and a
polynomial matrix $\mathbf{F} = [f_{i,j}] \in \mathbf{K}[x_1, \ldots, x_n]^{p
\times q}$, with $p \leq q$, we are interested in determining the isolated
points of $V_p(\mathbf{F},\mathbf{g})$, the algebraic set of points in
$\overline{\mathbf{K}}$ at which all polynomials in $\mathbf{g}$ and all
$p$-minors of $\mathbf{F}$ vanish, under the assumption $n = q - p + s + 1$.
Such polynomial systems arise in a variety of applications including for
example polynomial optimization and computational geometry. We design a
randomized sparse homotopy algorithm for computing the isolated points in
$V_p(\mathbf{F},\mathbf{g})$ which takes advantage of the determinantal
structure of the system defining $V_p(\mathbf{F}, \mathbf{g})$. Its complexity
is polynomial in the maximum number of isolated solutions to such systems
sharing the same sparsity pattern and in some combinatorial quantities attached
to the structure of such systems. It is the first algorithm which takes
advantage both on the determinantal structure and sparsity of input
polynomials. We also derive complexity bounds for the particular but important
case where $\mathbf{g}$ and the columns of $\mathbf{F}$ satisfy weighted degree
constraints. Such systems arise naturally in the computation of critical points
of maps restricted to algebraic sets when both are invariant by the action of
the symmetric group.
