In this paper, we investigate the global existence of classical solutions to the following parabolic–parabolic–parabolic two-species chemotaxis-competition system with singular sensitivity and nonlinear productions \(\begin{aligned} {\left\{ \begin{array}{ll} u_t=\Delta u-\chi _1\nabla \cdot (\frac{u}{w}\nabla w)+u(a_1-b_1u^{\alpha _1}-c_1v^{\beta _2}), & t>0,~x\in \Omega ,\\ v_t=\Delta v-\chi _2\nabla \cdot (\frac{v}{w}\nabla w)+v(a_2-b_2v^{\alpha _2}-c_2u^{\beta _1}), & t>0,~x\in \Omega ,\\ w_t=\Delta w-w+u^{\gamma _1}+v^{\gamma _2}, & t>0,~x\in \Omega , \quad \quad \quad \quad \quad \quad \quad \quad {\mathrm{(0.1)}}\\ \frac{\partial u}{\partial \nu }=\frac{\partial v}{\partial \nu }=\frac{\partial w}{\partial \nu }=0, & t>0,~x\in \partial \Omega ,\\ u(0,x)=u_0(x),~~v(0,x)=v_0(x),~~w(0,x)=w_0(x), & x\in \Omega , \end{array}\right. } \end{aligned}\) where \(\Omega \subset \mathbb {R}^N(N\ge 1)\) is a bounded convex smooth domain, \(\chi _i,b_i,c_i,\alpha _i,\beta _i>0,\) \(a_i\in \mathbb {R}\) and \(1\le \frac{8(\chi _1+\chi _2)}{8(\chi _1+\chi _2)-(\chi _1-\chi _2)^2}\le \gamma _i(i=1,2).\) If \(\gamma _1<\alpha _1\) and \(\gamma _2<\min \{\alpha _2,\beta _2\}\) or \(\gamma _2<\alpha _2\) and \(\gamma _1<\min \{\alpha _1,\beta _1\}\) or \(\begin{aligned} \chi _1+\chi _2\ge {\left\{ \begin{array}{ll} \frac{N\chi _1^2\gamma ^2_1}{2}=\min \Big \{\frac{N\chi _1^2\gamma ^2_1}{2},\frac{N\chi _2^2\gamma ^2_2}{2}\Big \}, & \text {if}~\frac{\chi _1}{\chi _2}\gamma _1<\gamma _2<\min \{\alpha _2,\beta _2\},\\ \frac{N\chi _2^2\gamma ^2_2}{2}=\min \Big \{\frac{N\chi _1^2\gamma ^2_1}{2},\frac{N\chi _2^2\gamma ^2_2}{2}\Big \}, & \text {if}~\frac{\chi _2}{\chi _1}\gamma _2<\gamma _1<\min \{\alpha _1,\beta _1\},\\ \frac{N\chi _1\chi _2\gamma _1\gamma _2}{2}, & \text {if}~\gamma _1<\alpha _1~\text {and}~\gamma _2<\alpha _2,\\ \frac{N\chi _1\chi _2\gamma _1\gamma _2}{2}\in \Big (\frac{N\chi _1^2\gamma ^2_1}{2},\frac{N\chi _2^2\gamma ^2_2}{2}\Big ), & \text {if}~\frac{\chi _1}{\chi _2}\gamma _1<\gamma _2<\alpha _2,\\ \frac{N\chi _1\chi _2\gamma _1\gamma _2}{2}\in \Big (\frac{N\chi _2^2\gamma ^2_2}{2},\frac{N\chi _1^2\gamma ^2_1}{2}\Big ), & \text {if}~\frac{\chi _2}{\chi _1}\gamma _2<\gamma _1<\alpha _1,\\ \max \Big \{\frac{N\chi _1^2\gamma ^2_1}{2},\frac{N\chi _2^2\gamma ^2_2}{2}\Big \},\\ \end{array}\right. } \end{aligned}\) we prove that for any given nonnegative initial functions \(u_0,v_0\in C^0(\overline{\Omega })\) and \(w_0\in W^{1,\infty }(\Omega ),\) problem (0.1) admits a global classical solution. Moreover, if \(\begin{aligned} \chi _1+\chi _2\ge \max \Big \{\frac{N\chi _1^2\gamma ^2_1}{2},\frac{N\chi _2^2\gamma ^2_2}{2}\Big \}, a_1<-\frac{1}{\gamma _1}~\text {and}~a_2<-\frac{1}{\gamma _2}, \end{aligned}\) we prove that the problem (0.1) has a global bounded classical solution.