I hope I'm just overlooking something. The Lagrangian is as follows: $$\mathcal{L}=\frac{1}{2}(\partial_\mu \phi)^2-(\frac{1}{2}m^2\phi^2+\frac{1}{3!}g\phi^3+\frac{1}{4!}\lambda\phi^4)$$ and I just want to calculate the saddle point approximation (SPA) for the 1PI effective action.
First, the counterterms from the "classical part" is $$\frac{1}{2}\delta_{m^2}\varphi^2+\frac{1}{3!}\delta_g\varphi^3+\frac{1}{4!}\delta_\lambda\varphi^4$$ and these are (up to the "source counterterm" which only enforces the one point function fluctuation to vanish, see Peskin.) the only ones available to counter anything from the SPA.
The divergent part of the SPA for a homogeneous background is given by (omitting the multiplicative constants here): $$(m^2+g\varphi+\frac{\lambda}{2}\varphi^2)^2/\epsilon.$$ This can be "seen" (if one chooses not to explicitly compute this) by dimensional analysis and the fact that the quadratic part of the lagrangian after expanding around $\varphi$ gives me this "mass term" (see Peskin, chapter 11.4 for an analogous calculation.)
Expanding this product gives me in particular this term here: $2m^2g\varphi/\epsilon.$
But this is weird: Which term (by counting $\varphi$) could counter this divergence? It does not seem to be the "source counterterm", as the one loop contribution for $\langle \eta \rangle$ (the one point fluctuation) does not have the same prefactors.
(They are given by something like $(g+\lambda\varphi)(m^2+g\varphi+\frac{\lambda}{2}\varphi^2)$ instead, so I cannot counter both using one $\delta J(x)$.)
Btw: I did compute everything.
EDIT: I did not forget the field renorm. It disappears due to assumption of homogeneity of $\varphi$. But even with this, i would still only get a quadratic contribution.
EDIT2: Regarding the source counterterm: It (i.e. the combination $\delta_J\phi$) vanishes when legendre transforming. So it does not appear anymore when evaluating the saddle point approximation (SPA)
EDIT3: Not even sure if this should be an edit or just answered directly: But a comment stated that the result is wrong. Well, the SPA is given by tr log $\frac{\mathcal{L}}{\delta \phi\delta\phi}$. And then you can just do the calculation to get the effective mass: $m^2+g\phi+\lambda/2 \phi^2$. And well, if you really want to think of the tadpole instead, remember we did a shift in the fields. So we generate another vertex from the original lambda interaction.
ALSO: If people claim the source counterterm does not vanish/ can be used: State this as an answer including derivation of the fact. This would deviate from standard literature such as P&S and Schwartz and would be very enlightening to see.