J'y suis allé totalement à la bourrin :
Si tu écrit que

alors
x^2)
donc
^2=u^2+uv+v^2\ (E_1))
De même
\vec{AB}+v\vec{AC})
implique que
^2=(u\!-\!1)^2+(u\!-\!1)v+v^2\ (E_2))
De même
\vec{AC})
implique que
^2=u^2+u(v\!-\!1)+(v\!-\!1)^2\ (E_3))
Ensuite,

et

te donnent deux équations linéaires en

et

(à coefficients "réellement" constant) qui te permettent de les déterminer en fonction de

et il n'y a plus qu'à réinjecter ça dans n'importe laquelle des 3 pour obtenir une équation bicarré en

dont la solution est celle ci dessus et dans laquelle apparait "magiquement" la formule de Héron donnant la surface du triangle de cotés

(le terme avec la racine).
EDIT :En fait, si on note

la surface du triangle équilatéral (ABC),

les surfaces des trois triangles équilatéraux de cotés respectifs

et

celle du triangle de cotés

, alors on a :

Et ça se "lit" directement sur une figure sur laquelle on fait 3 rotations comme celle de chan79 pour faire apparaitre à la fois

et

(le triangle ACA' est isométrique à ABD, idem pour les deux autres et ça explique que la surface totale coloriée soit le double de celle du triangle ABC)
