Aquí está la prueba anticipada por Stefan4024, con base en esta pregunta vinculada a su respuesta.
En primer lugar demostrar que si $p \equiv 1 \pmod 6$, luego
$$
p(a^2 + ab + b^2)^2 \,\mid\, (a+b)^p - ^p - b^p .
$$
Considere algunos fijos $b$, y deje $f(x)$ ser el polinomio
$$
f(x) = (x + b)^p - x^p - b^p.
$$
Deje $\omega = e^{2\pi i / 3}$ ser una primitiva de tercera raíz de la unidad, y se nota que desde $1 + \omega + \omega^2 = 0$, lo que tenemos que $1 + \omega = -\omega^2$ es un sexto de la raíz de la unidad.
Vamos a mostrar que el $\omega b$ es una raíz de ambos $f(x)$, y su derivado $f^\prime (x)$.
Tenemos que $f(\omega b)$ es igual a
$$
(\omega b + b)^p (\omega b)^p - b^p = b^p \left( (1 + \omega)^p - \omega^p - 1 \right)
$$
Desde $\omega^3 = (1 + \omega)^6 = 1$, e $p-1$ es divisible por $6$,$\omega^p = \omega$, e $(1 + \omega)^p = (1 + \omega)$. Así tenemos que
$$
f(\omega b) = b^p \left( (1 + \omega) - \omega - 1 \right) = 0.
$$
Por lo tanto $\omega b$ es una raíz de $f$. Del mismo modo,
$$
f^\prime (\omega b) = p (\omega b + b)^{p-1} - p(\omega b)^{p-1} = p b^{p-1} \left( (1 + \omega)^{p-1} - \omega^{p-1} \right) = p b^{p-1} (1 - 1) = 0.
$$
Por lo tanto $\omega b$ es una raíz de ambos $f$$f^\prime$, del que se desprende que $(x - \omega b)^2$ es un factor de $f$. Desde $f$ tiene coeficientes reales, $(x - \bar{\omega})^2$ es también un factor de $f$, y vemos que $(x^2 + xb + b^2)^2$ es un factor de $f$.
Ahora es un hecho bien conocido que para $1 \leq k \leq (p-1)$, que el coeficiente binomial $\binom{p}{k}$ es divisible por $p$, y así vemos que por el teorema del binomio de que todos los coeficientes de $f$ son divisibles por $p$. Por lo tanto $\frac{1}{p} f$ es un polinomio con coeficientes enteros, y es divisible por $(x^2 + xb + b^2)^2$, que es un monic polinomio con coeficientes enteros. De ello se sigue que podemos escribir
$$
\frac{1}{p} f(x) = (x^2 + xb + b^2)^2 \cdot g(x)
$$
donde $g(x)$ es de algún polinomio con coeficientes enteros. De esto se deduce fácilmente que el $f(a) = (a + b)^p - a^p - b^p$ es divisible por $p(a^2 + ab + b^2)^2$ si $p \equiv 1 \pmod 6$.
Ahora supongamos que $p \,\mid\, a^2 + ab + b^2$. Vamos a mostrar que el $p \equiv 1 \pmod 6$, por lo que
$$
p^3 \,\mid\, p(a^2 + ab + b^2)^2 \,\mid\, (a + b)^p - ^p - b^p .
$$
Tomamos nota de que
$$
p \,\mid\, 4a^2 + 4ab + 4b^2 = (2a + b)^2 + 3b^2.
$$
Si $p \,\mid\, b$, entonces podemos ver que también debemos tener ese $p \,\mid\, a$, y por lo $(a + b)^p - a^p - b^p$ es divisible por $p^p$, y así es, sin duda divisible por $p^3$. Supongamos ahora que $b$ no es divisible por $p$. Entonces tenemos que
$$
\left((2a + b) \cdot b^{-1} \right)^2 \equiv -3 \pmod p
$$
donde $b^{-1}$ es el inverso multiplicativo de a $b$ modulo $p$, y así vemos que el $-3$ es un residuo cuadrático módulo $p$.
Así
$$
\left( \frac{-1}{p} \right)\left( \frac{3}{p} \right) = \left( \frac{-3}{p} \right) = 1
$$
donde $\left( \frac{\;}{} \right)$ es el símbolo de Jacobi. Si $p \equiv 1 \pmod 4$, luego
$$
\left( \frac{-1}{p} \right) = 1
$$
y por la reciprocidad cuadrática,
$$
\left( \frac{3}{p} \right) = \left( \frac{p}{3} \right) =
\begin{cases}
1 & \text{ if } p \equiv 1 \pmod 3 \\
-1 & \text{ if } p \equiv 2 \pmod 3
\end{casos}.
$$
Vemos que en este caso, debemos tener la $p \equiv 1 \pmod 3$.
Por otro lado, si $p \equiv 3 \pmod 4$, entonces sabemos que
$$
\left( \frac{-1}{p} \right) = -1
$$
por lo tanto, debemos tener
$$
\left( \frac{3}{p} \right) = -1.
$$
Desde $3$ $p$ ambos $3$ mod $4$, la reciprocidad cuadrática en este caso nos da que
$$
-1 = \left( \frac{3}{p} \right) = -\left( \frac{p}{3} \right),
$$
y así volvemos a tener ese $p \equiv 1 \pmod 3$.
En cualquier caso, vemos que $p \equiv 1 \pmod 3$, y por lo $p \equiv 1 \pmod 6$, y el resultado de la siguiente manera.