Solve this Diophantine Equation: $a^3-b^3 = ab+61$ "solve this diophantine equation: $a^3-b^3 = ab+61$"
For Diophantine Equations I have always had the question in the form $ax+by = c$ and solved it via Extended Euclidean Algorithm. However, this question is not formatted in such a way. 
What I am thinking is that we move $ab$ to the LHS so that we can have a single constant on the RHS. Maybe then we could expand? I'm not sure. Hints/Advice is what I am looking for. 
Maybe we can substitute $a$ and $b$ with something else of a form I am more comfortable with? I'm not really sure. 

Finding all integer solutions to the equation $a^3-b^3=ab+61$ is equal to solving the diophantine equation $a^3-ab-b^3=61$, and placing $-b$ instead of $b$ gives $a^3+ab+b^3=61$.
Since the equation above is symmetric ( that is, the value doesn't changes even if we swap the position of the variables a and b), it can be expressed with the sum and product of a and b. Let c=a+b and d=ab. Then the quadratic equation $t^2-ct+d$ has two zeros, namely a and b. Thus, the discriminant $D=c^2-4d≥0$, which yields $d\le \frac {c^2}{4}$ ( and of course $c, d \in\mathbb Z$). In addition, $a^3+ab+b^3=61 \iff a^3+b^3+ab=61$
$\iff (a+b)^3-3(a+b)ab+ab=61$
$\iff c^3-3cd+d=61$
$\iff d=(c^3-61)/(3c-1)$
Since $d \in\mathbb Z$,  $27d = \frac{27c^3-1647}{3c-1} \in\mathbb Z$.
But then, $\frac{1646}{3c-1} \in\mathbb Z$, because $\frac{27c^3-1647}{3c-1} = 9c^2+3c+1\in\mathbb Z$. This means that $3c-1$ must divide $1646=2\cdot 823$ (823 is a prime). So $3c-1$ can be 1, 2, 823, 1646, -1, -2, -823, or -1646. Then $3c$ can be 2, 3, 824, 1647, 0, -1, -822, or -1645. Since c is an integer, the only possible values for c are c=0, 1, 549. This gives us the pair $(c,d) = (0,61), (1, -30), (549, 100528)$, and the only pair that $d\le \frac{c^2}{4}$ holds is $(1, -30)$. Thus, $c=1$ and $d=-30$, and it follows that a and b are the two zeros of the quadratic equation $x^2-x+30=0$. We now can conclude that the only solutions to the Diophantine equation $a^3+ab+b^3=61$ are $(a, b) = (6, -5), (-5, 6)$.
