I) We will assume that the real n\times n matrix A is not necessary symmetric. Therefore the Bosonic Gaussian integral in OP's last formula should read
\int_{\mathbb{R}^n} \! d^nx~ \exp \left[-\frac{1}{2}x^i A_{ij} ~x^j\right]
~=~\int_{\mathbb{R}^n} \! d^nx~ \exp \left[-\frac{1}{4}x^i (A+A^T)_{ij} ~x^j\right]
~=~\sqrt{\frac{(2\pi)^n}{\det(\frac{A+A^T}{2})}},\tag{1}
where we assume that the symmetric part of the matrix A is positive definite. Hence OP's last formula is basically an integral representation for
\frac{\det(A)}{\sqrt{\det(\frac{A+A^T}{2})}}.\tag{2}
II) OP may also be interested in the Grassmann integral representation of the Pfaffian
\tag{3} {\rm pf}(A) ~\propto ~\int \! d^n\theta~ \exp \left[\frac{1}{2}\theta^i A_{ij} ~\theta^j\right]~=~\int \! d^n\theta ~\exp \left[\frac{1}{4}\theta^i (A-A^T)_{ij} ~\theta^j\right],
where we leave it as an exercise to the reader to fix the correct normalization factor in eq. (3). In the second equality we have used that the Grassmann-numbers anticommute
\theta^i\theta^j~=~-\theta^j\theta^i.
If we use the integral representation (3) as a definition, then the Pfaffian {\rm pf}(A) only depends on the antisymmetric part of the n\times n matrix A. Let us therefore assume that the matrix A is antisymmetric from now on. One may then prove that the square of the Pfaffian is the determinant:
\tag{4} {\rm pf}(A)^2~=~\det(A).
This post imported from StackExchange Mathematics at 2015-02-07 08:27 (UTC), posted by SE-user Qmechanic