Yo no me centraría en la parametrización como tal. En lugar de ello, observaría que $4 a^3 = d^2 + 27 b^2.$ Esto da una restricción fácil de establecer en la factorización de primos de $a.$ Es necesario y suficiente que $a \geq 0$ y, siempre que cualquier primo $p | a$ y $p \equiv 2 \pmod 3,$ entonces el exponente de $p$ debe ser uniforme. Así que $a = 2$ o $a = 5$ o $a=10$ son imposibles. Sin el cubo, también habría una restricción en el primo 3, pero resulta que no importa porque se tiene $a^3.$ Sin el $4$ en $4a^3,$ habría una competencia entre $d^2 + 27 b^2$ y las demás formas del género, $4 u^2 \pm 2 u v + 7 v^2.$ Pero $$ \left( \begin{array}{cc} 4 & 0 \\ 1 & 1 \end{array} \right) \; \cdot \; \left( \begin{array}{cc} 1 & 0 \\ 0 & 27 \end{array} \right) \; \cdot \; \left( \begin{array}{cc} 4 & 1 \\ 0 & 1 \end{array} \right) = \; \;\; \; \left( \begin{array}{cc} 16 & 4 \\ 4 & 28 \end{array} \right) $$
Ahora, intente $a = 4 = 2^2,$ para que el exponente de 2 sea par, se obtiene $4a^3 = 256,$ y se obtiene $d = 16, b = 0,$ que parece demasiado fácil. pondré algo más abajo, tal vez salte cuadros... Pero, dada una legalidad $a$ como se ha descrito, el número de pares $d,b$ es finito.
$$a = 1,4a^3 = 4, ( d = \pm 2, b = 0), $$ $$a = 3, 4a^3 = 108, ( d = \pm 9, b = \pm 1), ( d = 0, b = \pm 2), $$ $$a = 7, 4a^3 = 1372, ( d = \pm 20, b = \pm 6), ( d = \pm 7, b = \pm 7), $$ $$a = 12, 4a^3 = 6912, ( d = \pm 72, b = \pm 8), ( d = 0, b = \pm 16), $$ $$a = 13, 4a^3 = 8788, ( d = \pm 70, b = \pm 12), ( d = \pm 65, b = \pm 13). $$
Bien, dadas todas las representaciones posibles $a = s^2 + s t + t^2,$ se puede construir todo $d,b,$ una tarea molesta a menos que $a$ es primo con $a \equiv 1 \pmod 3.$ Esto es, esencialmente, lo que la factorización de los enteros de Eisenstein le dará.
En realidad, esta última parte era innecesariamente pesimista. Con un poco de tratamiento especial de los primos 2,3, puedo ver cómo crear todas las representaciones posibles de $a = j^2 + 3 k^2,$ crear todas las representaciones posibles de $4a^3 = m^2 + 3 n^2,$ entonces sólo guarda los cuando $3 | n.$ Bastante complicado pero bastante fácil de programar.