Let $\Omega\subseteq \mathbb{R}^2$ be a bounded domain with smooth boundary. Prove that, for all $p>1$ and $1\leq q<\infty$ for all $f\in L^{p}(\Omega)$. Then there exists a unique $u\in H_0^1(\Omega)$ such that $\Delta u=u^q+f$ in $\Omega.$
This is not a linear pde, so I'm looking for a standard proof of the existence of such a weak solution. About the uniqueness, I think energy estimate may work, but I also find it hard to establish such an estimate.
Any solution or hint is highly appreciated.