It is very straightforward to use complex numbers for this issue.
Proof of the first formula (though not asked) : Let $n=2m$.
The roots of $x^{2m}-1=0$ are the 2m-th roots of unity ; they can be written in the following way :
$$\omega_k=e^{ik \pi/m}, \ \ \ \ \ k=0,1,... (2m-1) $$
Thus, we can decompose:
$$x^{2m}-1=(x-\omega_0)(x-\omega_1) \cdots (x-\omega_{2m-1})\tag{1}$$
Let us group these factors 2 by 2 in this way :
$$(x-\omega_k)(x-\omega_{2m-k})=x^2-(e^{ik\pi/m }+e^{i(2m-k)\pi/m })x+ e^{i \pi k/m}.e^{i(2m-k)\pi/m}$$
which gives, knowing that $e^{2i\pi}=1,$
$$x^2-\underbrace{(e^{ik\pi/m }+e^{-ik\pi/m})}_{2 \cos(k\pi/m)}x+ 1 \tag{2}$$
using Euler formula.
Then it suffices to use (2) in (1) giving identity
$$x^{2m}-1=(x^2-1)\prod_{k=1}^{m-1}(x^2-2x\cos \frac{k\pi}{m}+1)\tag{3}$$
i.e., your identity, because the product of factors $(x-1)(x-(-1))=(w-\omega_0)(x-\omega_m)$ has been set apart.
2nd formula :
Factorizing $x$ in each 2nd degree factor on the RHS of (3) gives :
$$x^m(x^{m}-\tfrac{1}{x^m})=x^m(x-\tfrac{1}{x})\prod_{k=1}^{m-1}(x-2\cos \frac{k\pi}{m}+\tfrac{1}{x})\tag{4}$$
Simplifying both sides by $x^{m}$, (4) becomes :
$$\frac{x^{m}-\tfrac{1}{x^m}}{x-\tfrac{1}{x}}=\prod_{k=1}^{m-1}(x-2\cos \frac{k\pi}{m}+\tfrac{1}{x})\tag{5}$$
identity $\frac{a^m-b^m}{a-b}=a^{m-1}+a^{m-2}b+...+b^{m-1}$ ($m$ terms) allows us to write the LHS of (5) under the form
$$x^{m-1} + x^{m-3} + ... \tfrac{1}{x^{m-3}}+\tfrac{1}{x^{m-1}} \ \ (m \ \text{factors})$$
Let us now take $m=16$. Let $x$ tends to $1$,
the LHS of (5) tends to $16=4^2$, whereas
the RHS of (5) tends to the product of $m$ factors. Using classical formula $1-\cos2a=2 \sin^2a$, this RHS is equal to :
$$4^{15}\Pi_{k=1}^{15}(\sin^2(k\pi/32))$$
Equating RHS and LHS, and taking the square root of both sides, we get the desired formula.