Complex polynomials won't work on their own, but a generalization can be made which allows one to produce an explicit formula (in stereographic coordinates) for any combination of indices which sum to $2$.
We can cover the Riemann sphere with two sets of opposite stereographic coordinates in the standard way, so that the transition function between them is $\tau(z)=1/z$, thinking of the tangent spaces $T_z\mathbb{C}$ as copies of $\mathbb{C}$ in the standard way, the differential of the transition function $d_z\tau:T_z\mathbb{C}\to T_{1/z}\mathbb{C}$ is just multiplication by $-1/z^2$.
We may treat a function $F:\mathbb{C}\to\mathbb{C}$ as a vector field in stereographic coordinates, and to see if it extends smoothly to the Riemann sphere, we can write its representative in the opposite coordinates, centered at infinity:
$$
\widetilde{F}(z)=-z^2F(1/z)
$$
We see that a complex polynomial of degree $>2$ will have a nonremovable singularity at infinity, but for degree $\le 2$ things behave as expected: each root has index equal to its algebreic multiplicity, and each polynomial of degree $d\in\{0,1\}$ has an additional zero of index $2-d$ at infinity. To get the other cases, we have to allow for negative indices; one way to do this is by complex conjugation.
Choose $a_1,\cdots,a_{n+2},b_1,\cdots,b_{n}\in\mathbb{C}\setminus\{0\}$, and define
$$
F(z)=\frac{(z-a_1)\cdots(z-a_{n+2})(\bar{z}-b_1)\cdots(\bar{z}-b_n)}{1+|z|^{2n}}
$$
This function is smooth on $\mathbb{C}$, and applying the transition function and doing a bit of algebra we obtain
$$
\widetilde{F}(z)=-a_1\cdots a_{n+2}b_1\cdots b_n\frac{(z-1/a_1)\cdots(z-1/a_{n+2})(\bar{z}-1/b_1)\cdots(\bar{z}-1/b_n)}{1+|z|^{2n}}
$$
And so $F$ extends smoothly to the Riemann sphere, and does not vanish at infinity. Additionally, the indices of the zeros of $F$ can be easily computed: for any zero $z_0\in\mathbb{C}$, the index of $z_0$ is the number of $a_i$ equal to $z_0$ minus the number of $b_i$ equal to $z_0$. This can be seen by first showing that the index of the vector field $V(z)=z^p\bar{z^q}$ is $p-q$ using polar decomposition, and looking at $F(z_0+z)$ in a small neighborhood around the origin.
Using this functional form, to obtain zeros of indices $I_1,I_2,\cdots,I_m\in\mathbb{Z}$ at distinct coordinates $z_1,\cdots,z_m\in\mathbb{C}$ we can choose positive integers $r_1,\cdots,r_m,s_1,\cdots s_n\in\mathbb{N}$ such that $r_i-s_i=I_i$ and then choose $a_1,\cdots,a_{n+2}$ to contain $r_i$ copies of each $z_i$ and $b_1,\cdots,b_n$ to contain $s_i$ copies of each $z_i$.