Alternatively, instead of expanding, one can obtain other linear dependences on the coefficients computing some derivatives at x=1 and at x=i in the above polynomial identity. (To this end, recall that the derivative at x=a of (x−a)mp(x) vanishes if m > 1 and it is just p(a) if m=1.)
Thus, for instance the first derivative at x=1 gives