J'ai une solution, mais c'est un poil bourrin.
Posons pour la suite Y = y^3 et réécrivons l'équation sous la forme Y = x^3-10^3.
On a alors (x-10)|Y, ie x = kY+10, k € Z. Avec ce changement de variable, après simplification, l'équation devient :
Y*(k^3*Y^2 + 30k²Y+300k-1)=0.
Si on exclut le cas Y=0 (dont on sait qu'il donne une solution), on est ramené à une équation du second degré ayant par hypothèse une racine réelle. Donc son delta doit être positif ou nul. Or, on a Delta = 4k^3(1-75k) < 0 si k est non nul. Mais k=0 redonne x=10 et y=0, ce qu'on avait exclu. Ainsi, l'équation admet bien pour seule solution le couple (10,0). 