Which polynomial's roots are its coefficients? Start with any polynomial of degree $n$ with complex coefficients, e.g.,
$$z^3+z^2+2 z+3 \;.$$
Find its $n$ roots, and list them in order of their modulus:
$$-1.28, (0.14\pm 1.53 i)$$
Now form a new polynomial with leading coefficient $1$ and the other
coefficients the $n$ roots:
$$
z^3
-1.28\, z^2
+(0.14 -1.53 i)\, z
+(0.14 +1.53 i) \;.
$$
Here I used the sorting convention that $(a-b i) < (a + b i)$.
A bit more formally, let $P$ be
$$
P \;:\; z^n + a_{n-1}z^{n-1} + \cdots + a_0 \;,
$$
with roots
$$r_{n-1},  \ldots ,r_0 \;,$$
sorted so that $|r_i| \le |r_{i-1}|$.
Define $P_1$ as
$$
P_1 \;:\; z^n + r_{n-1}z^{n-1} + \cdots + r_0 \;.
$$
My question is:

Q1. Which are the polynomials whose roots are its coefficients,
  in the sense that $r_i=a_i$, $i=0,\ldots,n-1$, i.e., $P_1 = P$?

Here is one fixed point:
$$
z^3
-(0.782599 +0.521714 i)\, z^2
+(0.884646 -0.589743 i)\, z
+(0.680552 +1.63317 i)
$$
When this process is iterated, it appears to fall into cycles, not always of length $1$.

Q2. Does iteration of this process always fall into a cycle?

(Added. 3Sep15.) The question suggested by Igor Rivin is also interesting,
especially in light of Richard Stanley's remark concerning real polynomials:

Q3. Which are the polynomials whose roots are its coefficients, in any order? (This is an unordered version of Q1.)

(Added. 4Sep15.)
Concerning Q2, which I realize is less interesting than the other questions, here is
an example of a cycle of length $11$. I plot the magnitude of the vector
of $3$ roots versus iteration #, starting from a random cubic polynomial. The last
vector of roots is
$$
(-1.52+0.88 i,-0.48-0.90 i,1.00 +0.02 i)
$$
with magnitude $2.27$:

          


 A: This is a long comment rather than a complete answer.
This is very similar to a multi-variable Julia set question,
which has the same flavour as the No wandering domain theorem.
The last sentence on wikipedia is key:
"However, the result can be generalized to many situations where the functions naturally belong to a finite-dimensional parameter space, most notably to transcendental entire and meromorphic functions with a finite number of singular values." 
This hints that the answer is YES for most points, in some sense.
However, I also strongly suspect that there is a "small" (measure 0) "Julia-like", repelling set of parameters, so that numerics prevent you from actually finding a chaotic sequence of steps, but this set is on the boundary of the basin of attraction of your cycles. 
In conclusion, my intuition says that there is an exceptional set of starting parameters, which never converges to a cycle, but it is a small set, and you will not find it by computer experiments with your method.
However, if you manage to run your map "backwards", this exceptional set should be attracting instead. That might be a way to get an idea of what parameters to pick to find such a point.
A: There’s certainly a whole line of periodic points of order two, namely the $n$-tuples of roots $(t,0,\cdots,0)$, giving the polynomial $X^n-tX^{n-1}$. More generally if $(\mathbf r)$ is an $n$-tuple giving a periodic point of order $m$, then it seems that $((\mathbf r),0)$ is an $(n+1)$ tuple giving a periodic point of the same order for the $(n+1)$-situation.
A: It is known (though I don't have a reference except for the statement at http://ocw.mit.edu/courses/mathematics/18-s34-problem-solving-seminar-fall-2007/assignments/roots.pdf, #2) that the only real polynomials with nonzero constant term whose roots equal their coefficients are $x^2+x-2$, $x^3+x^2-x-1$, and (approximately) 
  $$  x^3 + .56519772x^2 - 1.76929234x + .63889690 $$
  $$ x^4 + x^3 - 1.7548782x^2 - .5698401x + .3247183. $$
Update. I found a reference: P. Stein, On polynomial equations with coefficients equal to their roots, Amer. Math. Monthly 73 (1966), 272-274.
A: This is a partial answer to Q1. I have found the polynomials of degree 1, 2, and 3 using Vieta's formulas.
Let $|\alpha| < |\beta| < |\gamma|$ be the roots of the polynomial.
For degree 1: We get that the last coefficient is $-\alpha$ by Vieta's formula, and it is $\alpha$ by the criterium so $\alpha=0$ and the only polynomial is $P(x)=x$, which indeeds statistifies.

For degree 2: We get $\alpha\beta=\beta$, $-\alpha-\beta=\alpha$. From the first equation we get that $\beta=0$ or $\alpha =1$. $\beta=0$ gives $\alpha=0$. $\alpha=1$ gives $\beta=-2$.
Therefore the polynomials are $P(x)=x^2$ and $P(x)=x^2+x-2.$ They indeed statistify.

For degree 3:  We get $-\alpha\beta\gamma=\gamma$, $\alpha\beta+\beta\gamma+\alpha\gamma=\beta$ and $-\alpha-\beta-\gamma=\alpha$.
From the first equation we deduce that $\gamma=0$ or $\alpha\beta=-1$. If $\gamma =0$, then $\alpha\beta=\beta$ and $-\alpha-\beta=\alpha$. Therefore we have the same equation as in the first part, thus we get the polynomials are $P(x)=x^3$ and $P(x)=x^3+x^2-2x$. Only the former statistifies $|\alpha| < |\beta| < |\gamma|$.
If $\alpha\beta=-1$, then $\alpha=-\frac{1}{\beta}$. We also have $\gamma=-2\alpha-\beta=2\frac{1}{\beta}-\beta$. Therefore $-1+2-\beta^2+1-\frac{2}{\beta^2}=\beta$. This rewrites to $\beta^4+\beta^3-2\beta^2+2=0$.
Solutions for $\beta$ are (W|A verification): $\beta=-1$, $\beta=-1.7693$, $\beta=0.88465-0.58974i$, $\beta=0.88465+0.58974i$. From this we can get the values of $\beta$, and thus the values of $\alpha$ and $\gamma$.

$$P(x)=x^3+x^2-x-1$$
$$P(x)=x^3+0.565195x^2-1.7693x+0.638909$$
$$P(x)=x^3+(-0.78260-0.52171i)x^2+(0.88465-0.58974i)x+(0.68055+1.63316i)$$
$$P(x)=x^3+(-0.78260+0.52171i)x^2+(0.88465+0.58974i)x+(0.68055-1.63316i)$$

The second violates the ordering condition, but does statistify for Q3. The other three all statistify.
