As unit normal on the surface is $$\hat{n} = \frac{ax\hat{i}+by\hat{j}+cz\hat{k}}{\sqrt{a^2x^2+b^2y^2+c^2z^2}}$$ I took $$\overrightarrow{F} = \frac{\hat{i}}{ax}$$
so that $$\overrightarrow{F}\cdot \hat{n} = \frac{1}{\sqrt{a^2x^2+b^2y^2+c^2z^2}}$$
Hence the given integral is equal to
$$ \int \overrightarrow{F}\cdot dS = \int div \overrightarrow{F} dV = \int -\frac{1}{ax^2} dV $$ over ellipsoid.
I can't think of any method to compute this last integral. Please give some ideas.
You should choose suitable $\mathbf{F}$,
\begin{align} \mathbf{F} &= \frac{(x,y,z)}{ax^2+by^2+cz^2} \\ \nabla \cdot \mathbf{F} &= \frac{1}{ax^2+by^2+cz^2} \\ \end{align}
Using reduced coordinates:
\begin{align} (x,y,z) &= \left( \frac{x'}{\sqrt{a}}, \frac{y'}{\sqrt{b}}, \frac{z'}{\sqrt{c}} \right) \\ dV &= \frac{dx' \, dy' \, dz'}{\sqrt{abc}} \\ \iiint \nabla \cdot \mathbf{F} \, dV &= \iiint \color{red}{\frac{1}{x'^2+y'^2+z'^2}} \times \frac{dx' \, dy' \, dz'}{\sqrt{abc}} \\ &= \frac{1}{\sqrt{abc}} \int_0^{2\pi} d\phi \int_0^{\pi} \sin \theta \, d\theta \int_0^1 \color{red}{\frac{1}{r^2}} \times r^2 \, dr \\ &= \frac{4\pi}{\sqrt{abc}} \end{align}
Compare the results of another answer of mine here.