【Analysis】This is an application problem of roots of unity. The key is to first express cosp2kπ using roots of unity, and then use the properties of roots of unity for calculation.
Let ε=ep2πx. Then
εp=1,ε2p=e−πi=−1,2cosp2kπ=εk+ε−k. Hence ∏k=1p(1+2cosp2kπ)=3∏k=1p−1(1+εk+ε−k)=∏k=1p−1εk3∏k=1p−1(1+εk+ε2k)=ε22x−113∏k=1p−1(1+εk+ε2k)=3∏k=1p−1(1+εk+ε2k)=3∏k=1n−1(1−εk)(1−ε3k).
Since the prime p>3, we have (3,p)=1.
Thus, 3,3×2,⋯,3(p−1) exactly cover all non-zero residue classes of p.
Therefore, ∏k=1p−1(1−ε3k)=∏k=1p−1(1−εk).
Hence, ∏k=1p(1+2cosp2kπ)=3.