You can not (yet?) define polynomials over NN in Sage, because it is not a ring. What you can do is to define your polynomial over ZZ, find its factors (as elements of ZZ[x]), and then recombine them to find polynomials with non-negative entries.