I want to compute the discriminant of $x^{n+1}+x$.
It's easy to see the roots of this polynomial are $0$ and $e^{(2k+1)\pi i/n},k=0,1,\dots,n-1$. But it's still quite difficult, if we compute all $\sigma_k(c_1,\dots,c_n),k=1,2,\dots,n$ and use Newton's formula to compute all $s_k$.
Is there any easy way to do this?