Given the setup of a VAE, such as that outlined in Kingma and Welling (2014), there is a conditional probability distribution $p(z|x)$ describing the distribution of generated data $z \in \mathcal{Z}$ given sample data $x \in \mathcal{X}$, and there is likewise a conditional probability distribution $p(x|z)$ for the reverse. We do not know either of these distributions, so we approximate them using neural networks $q_{\phi}(z|x)$ and $p_{\theta}(x|z)$ with parameters $\phi$ and $\theta$. We then find $\phi$ and $\theta$ such that the ELBO (Evidence Lower BOund) is minimized.
Here is what confuses me:
$p(z|x)$ and $p(x|z)$ are conditional distributions. This means that for a fixed $x \in \mathcal{X}$, $\displaystyle \int_{\mathcal{Z}} p(z|x) \; dz = 1$, and similarly, for a fixed $z \in \mathcal{Z}$, $\displaystyle \int_{\mathcal{X}} p(x|z) \; dx = 1$. But we are approximating them using neural networks $q_{\phi}(z|x)$ and $p_{\theta}(x|z)$, which are in general just continuous functions parametrized by $\phi$ and $\theta$. So why is this a valid approximation when there is no requirement for a neural network to integrate to $1$ like $p(z|x)$ and $p(x|z)$ do?