For $x,y>0$ there is a unique solution $y(x)$ to $x = \frac{1}{3}y(y+1)(2y+1)^2(2y^2+2y+1)$ given by $$y=\tfrac{1}{2} 3^{-1/3} \sqrt{\frac{\left(\sqrt{11664 x^2-3}+108 x\right)^{2/3}+\sqrt[3]{3}}{(\sqrt{11664 x^2-3}+108 x)^{1/3}}}- \tfrac{1}{2}.$$$$y=\tfrac{1}{2} 3^{-1/3} \sqrt{\frac{\left(\sqrt{11664 x^2-3}+108 x\right)^{2/3}+3^{1/3}}{\bigl(\sqrt{11664 x^2-3}+108 x\bigr)^{1/3}}}- \tfrac{1}{2}.$$ Here is a plot of $y$ versus $x$.