In the case of quantum field theory:
First of all for a massless gauge field the most general form of the effective action will contain the renormalizable term
$ \mathcal{L} = -\frac{1}{4 g^2} F^{\mu \nu} F_{\mu \nu} $.
This follows simply from Lorentz Invariance (indices must be contracted properly essentially) and Gauge Invariance. I don't want to turn this into a post concerning Gauge Invariance simply because you didn't ask about that but I just want to say that while Gauge Invariance is a 'fake' symmetry it is a symmetry of the action none-the-less for our purposes so our effective Lagrangian should respect it. I also don't want to make this a post about renormalization so I just say that in general there will be higher order terms in the field strength ($ (F^{\mu \nu} F_{\mu \nu})^2$ and so on ) in the effective action but these are irrelevant both in the formal sense of the word in field theory and in the sense that we don't care about them here since you were really only asking about the quadratic part of the action.
If you just want the term quadratic in the gauge field $A^\mu$ it will take the form:
$ A_{\mu} \Pi^{\mu \nu} A_{\nu} $
which just follows from Lorentz invariance. In momentum space the only thing we have around that has indices on it is the 4-momentum, $p^\mu$ and the metric $\eta^{\mu \nu}$ so $\Pi^{\mu \nu}$ must take the form
$\Pi^{\mu \nu}(q^2) = \Pi (q^2) \left(\alpha \eta^{\mu \nu} + \beta q^\mu q^\nu \right) $.
Now from gauge invariance we know that $A^\mu \rightarrow A^\mu + q^\mu$ must be a symmetry or
$\alpha q^\nu + \beta q^2 q^\nu = 0 $
hence we can take $\alpha = q^2 $ and $\beta = -1$ and the remaining overall constant can be absorbed into $\Pi (q^2)$, giving
$\Pi^{\mu \nu}(q^2) = \Pi (q^2) \left( \eta^{\mu \nu} q^2 - q^\mu q^\nu \right) $.
For a $U(1)$ gauge field the $ A_{\mu} \Pi^{\mu \nu} A_{\nu} $ is the same things as $- \frac{1}{4 g_p^2} F^{\mu \nu} F_{\mu \nu} $ where $g_p $ is physical coupling constant. For non-abelian gauge fields there are 3 and 4 pt interactions as well in $ F^{\mu \nu} F_{\mu \nu} $ as well. So writing something like
$\mathcal{L} =- \frac{1}{4 g_p^2} F^{\mu\nu} F_{\mu \nu} + A_\mu \Pi^{\mu \nu} A_\nu $
is redundant, because $A_\mu \Pi^{\mu \nu} A_\nu$ is all already there in the $F^{\mu\nu} F_{\mu \nu}$.
For a massive vector field there is no gauge invariance and so my $\alpha$ and $\beta$ above are not constrained. In this case the effective action takes the form
$\mathcal{L}= - \frac{1}{4 g_p^2} F^{\mu\nu} F_{\mu \nu} - \frac{1}{2} m_p^2 A_\mu A^\mu $
and so part of the polarization tensor is not accounted for in the $F^{\mu\nu} F_{\mu \nu}$ term (where $m_p$ is the physical mass).
For a massless gauge field the polarization tensor takes the form
$\Pi_{ij} (q^2) = \Pi (q^2) \left( \delta_{ij} - \frac{q_i q_j}{q^2} \right)$
if one makes the gauge choice $A_0= 0 $ and $\ q_i A_i = 0 $.