Suppose $P(x)$ is a quadratic polynomial with integer coefficients satisfying the identity \[ P(P(x)) - P(x)^2 = x^2 + x + 2016 \] for all real $x$. What is $P(1)$?

Let's think step by step and output the final answer within \boxed{}.
user