Solution:
Let z1,…,z5 be the roots of Q(z)=z5+2004z−1. We can check these are distinct (by using the fact that there's one in a small neighborhood of each root of z5+2004z, or by noting that Q(z) is relatively prime to its derivative). And certainly none of the roots of Q is the negative of another, since z5+2004z=1 implies (−z)5+2004(−z)=−1, so their squares are distinct as well. Then, z12,…,z52 are the roots of P, so if we write C for the leading coefficient of P, we have
P(−1)P(1)=C(−1−z12)⋯(−1−z52)C(1−z12)⋯(1−z52)=[(i−z1)⋯(i−z5)]⋅[(i+z1)⋯(i+z5)][(1−z1)⋯(1−z5)]⋅[(1+z1)⋯(1+z5)]=[(i−z1)⋯(i−z5)]⋅[(−i−z1)⋯(−i−z5)][(1−z1)⋯(1−z5)]⋅[(−1−z1)⋯(−1−z5)]=(i5+2004⋅i−1)(−i5+2004⋅(−i)−1)(15+2004⋅1−1)(−15+2004⋅(−1)−1)=(−1+2005i)(−1−2005i)(2004)(−2006)=−20052+120052−1=−4020024/4020026=−2010012/2010013.
Alternative Solution: In fact, we can construct the polynomial P explicitly (up to multiplication by a constant). We write P(z2) as a polynomial in z; it must use only even powers of z and be divisible by z5+2004z−1, so we are inspired to try a difference of squares,
P(z2)=(z5+2004z−1)(z5+2004z+1)=(z5+2004z)2−12=z2(z4+2004)2−1,
giving
P(z)=z(z2+2004)2−1.
Now plugging in z=1 and z=−1 rapidly gives (20052−1)/(−20052−1) as before.