Let $$z=\cos\bigg(\frac{2\pi}{3^n+1}\bigg)+i\sin\bigg(\frac{2\pi}{3^n+1}\bigg)$$
Then $z^{3^n+1}=1$ and also $\displaystyle 2\cos \bigg(\frac{2\pi\cdot 3^k}{3^n+1}\bigg)=z^{3^k}+\frac{1}{z^{3^k}}$
Write $$\prod^{n}_{k=1}\bigg[1+2\cos\bigg(\frac{2\pi\cdot 3^k}{3^n+1}\bigg)\bigg]$$
$$=\bigg(1+z^3+\frac{1}{z^3}\bigg)\bigg(1+z^9+\frac{1}{z^9}\bigg)\cdots \cdots \bigg(1+z^{3n}+\frac{1}{z^{3n}}\bigg)$$
$$=\frac{(1+z^3+z^6)(1+z^9+z^{18})\cdots \cdots (1+z^{3^n}+(z^{3^n})^2)}{z^{3+9+\cdots \cdots +3^n}}$$
Multiply both Nr and Dr by $(1-z^3)$
$$=\frac{1-z^{3^{n+1}}}{(1-z^3)\cdot z^{3\frac{(3^n-1)}{2}}}= \frac{1-z^{-3}}{-(1-z^3)\cdot z^{-3}}=1.$$
$\text{Simplification}:\;\; $ From $z^{3n+1}=1\Rightarrow z^{3n}=z^{-1}$
And $$z^{3\frac{(3^n-1)}{2}}=z^{-3}\cdot z^{3\frac{(3^+1)}{2}}=z^{-3}\cdot \bigg(z^{\frac{(3^+1)}{2}}\bigg)^3=-z^3$$