Here's an alternate solution.
Lemma: Let $0\le j,k < n$. Then
$$
\sum_{j=0}^{n-1} \cos(2\pi k j/n) = \begin{cases}n,&k=0\\0,&\text{else}\end{cases}
$$
This is clear geometrically: the non-unit roots of unity sum to zero, and when $k=0$ there are $n$ terms of $1$.
Take $f(n)$ and multiply through by $2^n$:
$$2^n f(n) = \sum_{j=0}^{2n-1}{\left(2\cos\left( \frac{j \pi}{2n}\right)\right)^n \left( 2\cos \left( \frac{2 j \pi}{n} \right) + 1\right) \cos \left( \frac{j \pi}{2} - \frac{2 j \pi}{n} \right)}$$
$$= \sum_{j=0}^{2n-1}{\left(2\cos\left( \frac{j \pi}{2n}\right)\right)^n \left( 2\cos \left( \frac{4 j \pi}{2n} \right) + 1\right) \cos \left( \frac{nj \pi}{2n} - \frac{4 j \pi}{2n} \right)}$$
By symmetry, it suffices to sum over $0\le j \le n-1$ and then double:
$$= 2\sum_{j=0}^{n-1}{\left(2\cos\left( \frac{j \pi}{2n}\right)\right)^n \left( 2\cos \left( \frac{4 j \pi}{2n} \right) + 1\right) \cos \left( \frac{nj \pi-4j\pi}{2n} \right)}$$Denote $X=\frac{j\pi}{2n}$:
$$= 2\sum_{j=0}^{n-1}\left(2\cos( X)\right)^n \left( 2\cos (4X) + 1\right) \cos \left( (n-4)X\right)$$Now we are going to use the product to sum rules to turn these products into a sum of frequencies. This is going to be great for two reasons:
- By the lemma, we need only worry about the constant term, as all the other frequencies will die off.
- In fact, the constant term is a polynomial in $n$. Once we figure out its degree and sample that many points plus one, we can reverse-engineer it.
Claim: the constant term is a polynomial of degree four in $n$. Indeed, this is because of the factors of $4$ that appear in the arguments of the cosines.
Writing out the constants produces the sequence
$$\{1,2,4,8,16,31,57,99,\ldots\}$$
This is a well-known sequence, whose terms represent the number of regions when $n$ points in general position on a circle are connected, as seen here.
In summary,
$$ 2^n f(n) = 2n \left(\binom{n}{4}+\binom{n}{2}+\binom{n}{0}\right),$$ and the conclusion follows.