Soit $a,b,c\in\mathbb{R}$ de somme 1.
Soit $P$ le barycentre de $(M,a)$ et $(N,b)$ et $(O,c)$.
Il est facile de voir que les aires algébriques des triangles $NOP$, $OMP$ et $MNP$ sont proportionnelles aux nombres $a,b,c$.
Donc les 3 aires géométriques sont égales si et seulement si $|a|=|b|=|c|$.
Autrement dit il y a 4 points magiques
$$P_1[/latex] barycentre de [latex](M,1/3)[/latex], [latex](N,1/3)[/latex], [latex](O,1/3)$$
$$P_2[/latex] barycentre de [latex](M,1)[/latex], [latex](N,1)[/latex], [latex](O,-1)$$
$$P_3[/latex] barycentre de [latex](M,1)[/latex], [latex](N,-1)[/latex], [latex](O,1)$$
$$P_4[/latex] barycentre de [latex](M,-1)[/latex], [latex](N,1)[/latex], [latex](O,1)$$
Voilà !

Pour le 1), je pense que seul le centre de gravité (point d'intersection des médianes) du triangle convient. Mais je n'en suis pas sûr.

2) La 3D, c'est terre inconnue pour moi.

Coucou

Une médiane d'un triangle partage ce dernier en 2 triangles de même aire.
Les 3 médianes étant concourantes. Le point de concours G partage un triangle en 6 triangles de même aire.
En associant par 2 ces 6 triangles, on obtient 3 triangles de même aire.
Donc dans le plan, le centre de gravité est le point magique d'un triangle.

Soit M' le symétrique du point M par rapport au milieu du segment NO.
Le quadrilatère MNM'O est un parallélogramme (ses diagonales se coupent en leur milieu) donc ses diagonales le partagent en quatre triangles de même aire puisque
Les trois triangles  MOM' , NOM' et MPM' sont chacun formés de 2 de ces triangles donc A(MOM') = A(NOM') = A (MPM')

On obtient de le même façon 2 autres points magiques en prenant les symétriques des points N et O par rapport au milieu de leur côté opposé.

Dans l'espace, je n'ai rien vérifié mais intuitivement je dirais que la droite orthogonale au plan (MNO) passant par G nous donne l'ensemble des points magiques.

Question 2  : déterminer l'ensemble des points magiques dans l'espace 3D.

On se place dans $\mathbb{R}^3$ et l'on note  $Oxyz$ le repère orthonormé canonique. Pour simplifier on va traiter le cas particulier où les points $M,N,O$ sont donnés par
$$M=(1,0,0),\ N=(0,1,0),\ O=(0,0,0)\,.$$
On pose $P=(x,y,z)$.

On rappelle que si l'on a un triangle $ABC$ dans l'espace, son aire s'obtient à partir d'un produit vectoriel par la formule
$$\mathrm{Aire}(ABC)=\frac{\ \Vert\overrightarrow{AB}\wedge\overrightarrow{AC}\Vert\ }{2}$$
Cette formule permet de montrer que
$$4\,\mathrm{Aire}(MNP)^2=x^2+y^2+2z^2+2xy-2x-2y+1$$
$$4\,\mathrm{Aire}(NOP)^2=x^2+z^2$$
$$4\,\mathrm{Aire}(OMP)^2=y^2+z^2$$
Par suite l'ensemble des points magiques est donné par les solutions du système
$$y=\pm x$$
$$y^2+z^2+2xy-2x-2y+1=0$$
Pour $y=x$, on doit prendre l'intersection du plan $y=x$ avec une quadrique, ce qui donne une conique. De même pour $y=-x$. On voit aisément que ces 2 coniques sont disjointes ; de plus elles ne sont pas réduites à l'ensemble vide d'après la Question 1.

Finalement l'ensemble des points magiques est la réunion de 2 coniques disjointes. Avec un peu de patience, on peut chercher la nature de chacune de ces coniques.

On peut aussi étendre la méthode à un triangle $MNO$ quelconque. Dans ce cas l'ensemble des points magiques est une courbe obtenue comme intersection de 2 quadriques.

Voilà !

Un petit message pour indiquer que si le triangle $MNO$ est équilatéral, alors l'ensemble des points magiques est composé de 4 droites parallèles, perpendiculaires au plan $(MNO)$. La position des 4 droites résulte de la Question 1.
Il est facile de vérifier géométriquement que ces 4 droites font partie du lieu.