Aquí es una manera de llegar al resultado. Utilizamos la simple regularización de la Delta de Dirac dada por
$$\delta_\epsilon(x)=\begin{cases}\frac1{2\epsilon}&,x\in[-\epsilon,\epsilon]\\0&,\text{elsewhere} \end{cases}$$
donde $\lim{\epsilon\to0}\delta\epsilon(x) \sim \delta(x)$ en el sentido de distribuciones. Es decir, para cada función lisa con la ayuda compacta (es decir, función de prueba), tenemos
$$\lim{\epsilon\to0}\int{-\infty}^\infty f(x)\delta_\epsilon(x)\,dx=f(0)$$
Entonces, tenemos
$$\begin{align} \int{-1}^1 x^2\delta(x^3)\,dx&=\lim{\epsilon\to0}\int{-1}^1 x^2 \delta\epsilon(x^3)\,dx\\ &=\lim{\epsilon\to0}\int{-\sqrt[3]{\epsilon}}^\sqrt[3]{{\epsilon}}x^2\left(\frac1{2\epsilon}\right)\,dx\\ &=\frac13 \end {Alinee el} $$
como era de esperar!
Para la observación, la regularización que utilizamos aquí en no único.