In this paper we study some aspects of the integrability problem for polynomial vector fields x = P(x, y), = Q(x, y). We analyze the possible existence of first integrals of the form I(x, y) = (y - g(1) (x))(alpha1) (y - g(2)(x))(alpha2...) (y -g(e)(x))(alphae) h(x), where g(1)(x),..., g(e)(x) are unknown particular solutions of dy/dx = Q(x, y)/P(x, y), alpha(i) are unknown constants and h(x) is an unknown function. We show that for certain systems some of the particular solutions remain arbitrary and the other ones are explicitly determined or are functionally related to the arbitrary particular solutions. We obtain in this way a nonlinear superposition principle that generalize the classical nonlinear superposition principle of the Lie theory. In general, the first integral contains some arbitrary solutions of the system but also quadratures of these solutions and an explicit dependence on the independent variable. In the case when all the particular solutions are determined, they are algebraic functions and our algorithm gives an alternative method for determining such type of solutions.