Once again, I failed to make a concise post so feel free to skip to the emphasized parts.
In the context of symplectic and contact geometry, I would like to establish the following linear algebra fact:
Proposition. Let $E$ be a real vector space of dimension $2n$ and let $\omega\in\Lambda^2 E^*$, then $\omega$ is non-degenerate if and only if $\omega^n$ is a volume form of $E$.
Proof. There exists $(\omega_{i,j})_{1\leqslant i,j\leqslant 2n}$ a skew-symmetric $2n\times 2n$ matrix with real entries such that, one has: $$\omega=\sum_{i,j=1}^{2n}\omega_{i,j}\,\mathrm{d}x_i\wedge\mathrm{d}x_j.$$ A straightforward computation then leads to the following equality: $$\omega^n=\left(\sum_{\sigma\in\mathfrak{S}_{2n}}\varepsilon(\sigma)\prod_{i=1}^n\omega_{\sigma(2i-1),\sigma(2i)}\right)\mathrm{d}x_1\wedge\cdots\wedge\mathrm{d}x_{2n},$$ where $\mathfrak{S}_{2n}$ stands for the permutation group on the set $\{1,\ldots,2n\}$ letters and $\varepsilon(\sigma)$ is the signature of $\sigma$. Therefore, let us define the following quantity associated to $\omega$:
$$\operatorname{vol}(\omega):=\sum_{\sigma\in\mathfrak{S}_{2n}}\varepsilon(\sigma)\prod_{i=1}^n\omega_{\sigma(2i-1),\sigma(2i)}.$$
The proof boils down to show that $(\omega_{i,j})_{1\leqslant i,j\leqslant 2n}$ is invertible if and only if $\operatorname{vol}(\omega)$ is non-zero. $\Box$
Hence, my question arises:
Question. How can I show that $\operatorname{vol}(\omega)$ is non-zero if and only if $(\omega_{i,j})_{1\leqslant i,j\leqslant 2n}$ is invertible?
I am aware that $\operatorname{vol}(\omega)=2^nn!\operatorname{pf}(\omega)$, where $\operatorname{pf}(\omega)$ is the Pfaffian of $\omega$ and that one has: $$\det(\omega)=\operatorname{pf}(\omega)^2,\tag{1}$$ which immediately leads to the desired result, but $(1)$ is a pain to establish and I would like to avoid using it. If you are aware of a clever and rather elementary way to prove $(1)$, I will also be happy with it.
Any enlightenment will be greatly appreciated.