Найти целые неотрицательные решения уравнения

В тех случаях, когда многочлен симметричен относительно $%x$% и $%y$%, удобно выражать его через новые переменные $%a=x+y$% и $%b=xy$%. В данном случае $%x^3+y^3=(x+y)(x^2-xy+y^2)=a(a^2-3b)$%, поэтому уравнение принимает вид $%a^3-3ab-3b-a=0$%, то есть $%a(a^2-1)=3b(a+1)$%. Ввиду неотрицательности $%x$% и $%y$%, число $%a+1$% положительно, и на него можно поделить обе части равенства. Получится $%a(a-1)=3b$%. При этом числа $%x$%, $%y$% станут корнями квадратного уравнения $%t^2-(x+y)t+xy=0$%, то есть $%t^2-at+b=0$%. Дискриминант должен быть неотрицателен, что означает $%D=a^2-4b\ge0$%. Тем самым, $%4a(a-1)=12b\le3a^2$%. При $%a=0$% имеем $%b=0$%, что соответствует нулевому решению $%x=y=0$%. Если $%a\ne0$%, то сокращаем на $%a$%. Знак неравенства при этом не меняется ввиду $%a > 0$%. Имеем $%4a-4\le3a$%, то есть $%a\le4$%. Дальнейшее просто: значение $%a=2$% не подходит, так как $%a(a-1)$% должно быть кратно трём, а для остальных значений получается $%a=1$%, $%b=0$%; $%a=3$%, $%b=2$%; $%a=4$%, $%b=4$%. Для каждой из этих пар значений решается система $%x+y=a$%, $%xy=b$%. Общее количество решений системы в целых неотрицательных числах равно шести: $%(x;y)\in\{(0;0),(1;0),(0;1),(1;2),(2;1),(2;2)\}$%.