At a fixed time (set $t = 0$), we can write the free complex scalar field as \begin{equation*} \begin{split} \phi(x) &= \int\frac{d^3p}{(2\pi)^3}\frac{1}{\sqrt{2E_p}} \left[b_{-p}^\dagger+a_p\right] e^{ipx}\\ &=\int\frac{d^3p}{(2\pi)^3}\frac{1}{\sqrt{2E_p}} \left[a_pe^{ipx}+b_p^\dagger e^{-ipx}\right]\\ \phi^{\dagger}(x) &= \int\frac{d^3p}{(2\pi)^3}\frac{1}{\sqrt{2E_p}} \left[a^{\dagger}_{-p}+b_p\right] e^{ipx} \end{split} \end{equation*}
The parity reversal operator $P$ flips the sign of the space coordinate $x$. Its transformation rules (on a Dirac field) are \begin{equation*} Pa_p^sP = a_{-p}^s\quad\text{and}\quad Pb_p^sP = -b_{-p}^s\qquad\Rightarrow\qquad P\psi(t,x)P = \gamma^0\psi(t,-x) \end{equation*}
My question is how does the operator affect the $\phi$ field? Is it right we have
\begin{equation*} \begin{split} P\phi(x)P &= \int\frac{d^3p}{(2\pi)^3}\frac{1}{\sqrt{2E_p}} \left[Pb_{-p}^\dagger P+Pa_pP\right] e^{-ipx}\\ &= \int\frac{d^3p}{(2\pi)^3}\frac{1}{\sqrt{2E_p}} \left[-b_{-p}^\dagger +a_{-p}\right] e^{-ipx} \end{split} \end{equation*} I'm not quite sure how to further simplify this.
Now suppose we have the pseudoscalar Yukawa Lagrangian: $$ L = \frac{1}{2}\partial_\mu\phi\partial^\mu\phi-\frac{1}{2}m^2\phi^2+\bar\psi(i\not\partial-m)\psi-g\gamma^5\phi\bar\psi\psi. $$
It is invariant under the parity transformation: $$P\psi(t,x) P = \gamma^0\psi(t,-x), \qquad P\phi(t,x) P = -\phi(t,-x).$$ If we replace the pseudoscalar term with $g\phi\bar\psi\psi$ instead, it seems like the parity transformation on $\phi$ should be $P\phi(t,x) P = \phi(t,-x)$ to make $L$ invariant. Is that correct? If so, how can I understand the difference of $P\phi(t,x) P$ in two cases, based on the expansion of $\phi$ above?