For the special case where all $c_{i,j}$'s are equal to 1 and $m=2n-1$, take a look at Ramanujam's paper: http://ramanujan.sirinudi.org/Volumes/published/ram03.pdf. Needless to say, it is an ingenious method. The steps are the following:
(1) The key idea is to recognize that the coefficients (w.r.t $\theta$, upto the $2n-1$ order) of $\sum_{k=1}^n\frac{x_k}{1-\theta y_k}$ in the series expansion would be the LHS of the nonlinear equations.
(2) Now, $\sum_{k=1}^n\frac{x_k}{1-\theta y_k} = d_1 + d_2\theta + \cdots + d_{2n}\theta^{2n-1} + \cdots = \frac{A_1 + A_2\theta + \cdots + A_{n-1}\theta^{n-1}}{B_1 + B_2\theta + \cdots + B_{n-1}\theta^{n-1}}$.
(3) Multiplying by $B_1 + B_2\theta + \cdots + B_{n-1}\theta^{n-1}$ on either side and comparing coefficients, gives linear equations in $A_i$'s and $B_i$'s.
(4) If a solution exists, one can then determine partial fractions (in $\theta$), and the coefficients would be the answers.
For your case, you might want to think on the same lines with Pade approximants.
All this is only a suggestion for the most general case. Hope it helps.