[guided]We now use the analytic division theorem that is designed for exactly this situation: Skoda's $L^2$ division theorem. In the form needed here, it says that if $D\subset\mathbb{C}^m$ is pseudoconvex, $g_1,\dots,g_p\in\mathcal O(D)$, $Q=\min\{m,p-1\}$, $\alpha>1$, $\psi$ is plurisubharmonic on $D$, and $H\in\mathcal O(D)$ satisfies
\begin{align*}
\int_D |H|^2 |g|^{-2(\alpha Q+1)} e^{-\psi}\,d\mathcal L^{2m}<\infty,
\end{align*}
then there are $v_1,\dots,v_p\in\mathcal O(D)$ such that
\begin{align*}
H=\sum_{j=1}^p g_jv_j
\end{align*}
and
\begin{align*}
\int_D \sum_{j=1}^p |v_j|^2 |g|^{-2\alpha Q} e^{-\psi}\,d\mathcal L^{2m}
\leq
\frac{\alpha}{\alpha-1}
\int_D |H|^2 |g|^{-2(\alpha Q+1)} e^{-\psi}\,d\mathcal L^{2m}.
\end{align*}
We apply this theorem on $D=\Omega'$ with the holomorphic tuple $g_1,\dots,g_p$, the holomorphic dividend
\begin{align*}
F:\Omega' &\to \mathbb{C} \\
z &\mapsto f(z)^{q+1},
\end{align*}
the parameter
\begin{align*}
q=\min\{m,p-1\},
\end{align*}
and the plurisubharmonic weight $\psi=0$.
Let us check the hypotheses. The domain $\Omega'$ is pseudoconvex by assumption. The functions $g_1,\dots,g_p$ are holomorphic on $\Omega'$ because they are holomorphic on the larger domain $\Omega$. The function $F=f^{q+1}$ is holomorphic on $\Omega'$ because powers of holomorphic functions are holomorphic. The weight $\psi=0$ is plurisubharmonic. Finally, the required weighted integral
\begin{align*}
\int_{\Omega'} |F|^2 |g|^{-2(\alpha q+1)}\,d\mathcal L^{2m}
\end{align*}
is finite by the preceding step.
Skoda's theorem therefore produces holomorphic functions
\begin{align*}
u_1,\dots,u_p:\Omega' \to \mathbb{C}
\end{align*}
such that
\begin{align*}
F=\sum_{j=1}^p g_j u_j
\qquad \text{on } \Omega',
\end{align*}
and such that the weighted $L^2$ estimate
\begin{align*}
\int_{\Omega'} \sum_{j=1}^p |u_j|^2 |g|^{-2\alpha q}\,d\mathcal L^{2m}
\le
\frac{\alpha}{\alpha-1}
\int_{\Omega'} |F|^2 |g|^{-2(\alpha q+1)}\,d\mathcal L^{2m}
\end{align*}
holds.
Substituting back $F=f^{q+1}$ gives the division identity
\begin{align*}
f^{q+1}=\sum_{j=1}^p g_j u_j
\qquad \text{on } \Omega',
\end{align*}
and the estimate becomes
\begin{align*}
\int_{\Omega'} \sum_{j=1}^p |u_j|^2 |g|^{-2\alpha q}\,d\mathcal L^{2m}
\le
\frac{\alpha}{\alpha-1}
\int_{\Omega'} |f|^{2(q+1)} |g|^{-2(\alpha q+1)}\,d\mathcal L^{2m}.
\end{align*}
This is the effective part of the theorem: the coefficient functions are not merely shown to exist, but are produced with an explicit weighted norm bound.[/guided]