The maximum of a polynomial on the unit circle
Let's create a proof a la Koosis. All the techniques used below can be found in his book "The Logarithmic Integral".
Take $a>1$ and put $f(z)=a\prod_j(1+z/z_j)^{p_j}$. That is a nice analytic function and, while its absolute value is somewhat hard to understand, its argument is very simple: it is an $n$-piece piecewise linear function on the circle with slope $\frac 12 n$ (with respect to the usual circle length) and jumps at $z_j$. Let $I$ be the image of one of the arcs between two adjacent points $z_j$ and $z_{j+1}$ under the mapping $z\mapsto \operatorname{arg}f(z)$. We can transplant all functions defined on the circle arc $[z_j,z_{j+1}]$ to $I$ using this mapping. Note that the integral of any function over the arc with respect to the circle length is just $2/n$ times the integral of its transplant over $I$ with respect to the line length.
Let $\Phi$ be the transplant of $|f|$. Assume that $\Phi<2$ on $I$. The transplant of $f$ is then just $F(t)=\Phi(t)e^{it}$. The key observation is the following: $$ \int_I \log|2-F(t)|dt\ge \log 2(|I|-\pi). $$ Assuming that it is true, we conclude that the full integral of $\log|2-f|$ over the unit circle is at least $2/n$ times the sum of the right hand sides over the intervals corresponding to all arcs, which is $0$. On the other hand, if $a>1$, then $\log|2-f(0)|=\log|2-a|<0$, so $2-f$ must have a root inside the disk and the maximum principle finishes the story.
Now let us prove the observation claim. The only thing we really know about $\Phi$ is that it is log-concave and, thereby, unimodal. Fortunately, that's all we need. So, in what follows, $\Phi$ will be just any unimodal function on $I$ with values in $[0,2]$. Since we can always extend $\Phi$ by $0$ outside $I$, we can switch to any larger interval we want without making the inequality easier. So, WLOG, $I=[-2\pi n-\frac\pi 2,2\pi n+\frac\pi 2]$
Now, let us observe that for every fixed $t$, the integrand is minimized for $\Phi(t)=2\max(0,\cos t)$ (that is just the nearest point on the line) and that the farther we go away from this optimal value, the larger the integrand is. Therefore, to minimize the left hand side, we need to stay as close to the black regime on the picture (the graph of $2\cos_+ t$) as we can.
Suppose that the actual $\Phi$ is given by the blue line. Then, replacing $\Phi$ by the red line $\Psi$, we come closer to the optimum at every point. But the red line consists of several full periods (horizontal pieces) and several pieces that together constitute one full positive arc of $\cos t$. Now, each full period means running over some circle around the origin, so the average value of $\log|2-\Psi(t)e^{it}|$ over each full period is exactly $\log 2$. At last, the $2\cos t$ part gives $\int_0^\pi\log (2|\sin t|)dt=0$, which is exactly the loss of $\pi \log 2$ compared to $\log 2$ times its length $\pi$.
That's it. Feel free to comment and/or ask questions.
Here is a second proof. Let $$H(z)=\prod_j|z-z_j|^{p_j}.$$ Let $$F(z)=\prod_j\left((z-z_j)(1/z-\bar{z}_j)\right)^{p_j/n},\quad |z|\geq 1.$$ It is easy to see that: $$|F(z)|=H^{2/n}(z),\quad |z|=1,$$ and that $F$ maps conformally the exterior of the unit disc on the exterior of a "star" consisting of $n$ straight segments $[0,a_j]$ with some complex $a_j$. Moreover, $F(z)\sim z$ as $z$ tends to infinity.
Now we see that the result is equivalent to the following theorem of Dubinin: Let $K$ be the union of some intervals of the form $[0,a_j]$, and suppose that capacity of $K$ is $1$. Then $\max_j|a_j|$ takes its minimal value when the star is symmetric: all |a_j| are equal and the angles between adjacent intervals are equal.
The reference is Dubinin, Uspekhi Mat. Nauk (=Russian Math. Surveys), 49 (1994). Fedja's proof is much simpler, so it is a new proof of Dubinin's theorem.