Find the all positive integer solutions $(a,b)$ to $\frac{a^3+b^3}{ab+4}=2020$. 
Find the all positive integer solutions of given equation
  $$\frac{a^3+b^3}{ab+4}=2020.$$

I find two possible solutions, namely $(1011,1009)$ and $(1009,1011)$, but the way I solve  the equation was messy and I don't know if there are any other solutions.
Source: Turkey $1.$ TST for IMO $2020$

Note that if $p\mid a^2-ab+b^2$, where $p$ is a prime natural number s.t. $p\equiv 2\pmod{3}$, then $p\mid a$ and $p\mid b$.  For $p=2$, the claim is easily seen by inspection.  Let now $p>2$.  We prove by contradiction.  Suppose that $p\nmid a$ or $p\nmid b$.  It follows immediately that $p\nmid a$ and $p\nmid b$. Since
$$4(a^2-ab+b^2)=(2a-b)^2+3b^2\equiv0 \pmod{p},$$
we have $x^2\equiv-3\pmod{p}$, where $x=(2a-b)c$ if $c$ is an inverse of $b$ modulo $p$.  Consequently,
$$\left(\frac{-3}{p}\right)=1.$$
By quadratic reciprocity,
$$\left(\frac{p}{3}\right)=\left(\frac{p}{-3}\right)=1.$$
Hence $p\equiv 1\pmod{3}$, which is a contradiction.
From $\frac{a^3+b^3}{ab+4}=2020$, we get
$$2020(ab+4)=a^3+b^3=(a+b)(a^2-ab+b^2).$$
If $101\mid a^2-ab+b^2$, then $101\mid a$ and $101\mid b$ by the paragraph above.  Thus $101^3\mid (a+b)(a^2-ab+b^2)$, but clearly $101^3\nmid 2020(ab+4)$.  Hence, $101\mid a+b$.
Similarly, $5\mid a+b$.
If $2\mid a^2-ab+b^2$, then $2\mid a$ and $2\mid b$.  Write $a=2u$ and $b=2v$, then
$$1010(uv+1)=(u+v)(u^2-uv+v^2).$$
If $2\mid u^2-uv+v^2$, then $2\mid u$ and $2\mid v$.  Therefore $8\mid (u+v)(u^2-uv+v^2)$, but clearly $8\nmid 1010(uv+1)$.  Hence $2\mid u+v$.  Consequently $4\mid a+b$.  On the other hand, if $2\nmid a^2-ab+b^2$, then $4\mid a+b$ as $4\mid 2020$.  In any case, $4\mid a+b$.  This means $$a+b=4\cdot 5\cdot 101 \cdot k=2020 k$$ for some positive integer $k$.
Hence
$$k=\frac{a+b}{2020}=\frac{ab+4}{a^2-ab+b^2}.$$
Clearly $a\ne b$.  Since $a+b=2020k$, $a\equiv b\pmod{2}$ so $(a-b)^2\ge 4$.  If $(a-b)^2>4$, then $a^2-2ab+b^2>4$, making $ab+4<a^2-ab+b^2$, so
$$k=\frac{ab+4}{a^2-ab+b^2}<1,$$
which is a contradiction.  Hence, $(a-b)^2=4$, so that $ab+4=a^2-ab+b^2$, making $k=1$.  Thus we have $a+b=2020k=2020$ and $a-b=\pm\sqrt4=\pm2$.  This gives 
$$(a,b)=(1009,1011)\vee (a,b)=(1011,1009),$$
and both are solutions.  In fact, these two solutions are the only integer solutions (positive or negative) to the required equation.
