Proof: Let f(x)=g(x)h(x), where g,h are polynomials with integer coefficients, not constants, and the leading coefficients are 1, then xn≡g(x)h(x)(mod2). Therefore, we can set
g(x)=xr+2(br−1xr−1+⋯+b0),h(x)=xn−r+2(cn−r−1xn−r−1+⋯+c0),1⩽r<n.
By comparing the constant terms, we get
b0+xn−r=c0=±1,xr⋅2(cn−r−1xn−r−1+⋯+c0)⋅2(br−1xr−1+⋯+b0)≡0(mod4).
By comparing the coefficients of the lowest degree terms, we know r=n−r, i.e., n=2r. Thus,
(cr−1+br−1)xr−1+⋯+(b1+c1)x≡
0(mod2),b1+c1≡0(mod2). The coefficient of the linear term in g(x)h(x) is 4(b1c0+c1b0)=4c0(b1+c1)≡ 0(mod8), which contradicts the fact that the coefficient of the linear term in f is 4.