(3x+1)^3 - x^3
= (3x+1-x)[(3x+1)^2 + x^2 + x(3x+1)]
= (2x+1)[9x^2+6x+1+x^2+3x^2+x]
= (2x+1)[13x^2 + 7x + 1]