From SL Loney Trigonometry:
Prove that $$r^2+r_1^2+r_2^2+r_3^2=16R^2-a^2-b^2-c^2$$ for a triangle with sides $a$, $b$, $c$, circumradius $R$, inradius $r$, and exradii $r_1$, $r_2$, $r_3$.
This one is a really dirty question.
I have tried a lot of stuff like $r=4R\times \Pi(\sin\frac{A}{2}$) and $r_1=4R\times\sin\frac{A}{2}\cos\frac{B}2\cos\frac{C}2$ etc..
There are other ways like $-r+r_1+r_2+r_3=4R$, but all these are way too cumbersome or clumsy.
Anything more intuitive?
we Need the following formulas: $$r_1=\frac{\Delta}{s-a}$$ etc $$\Delta=rs$$, $$\Delta=\frac{abc}{4R}$$, $$s=\frac{a+b+c}{2}$$ so the right-hand side of your equation is given by $$\frac{8abc)^2}{s(s-a)(s-b)(s-c)}-a^2-b^2-c^2=-{\frac {{a}^{6}-{a}^{4}{b}^{2}-{a}^{4}{c}^{2}-{a}^{2}{b}^{4}+10\,{a}^ {2}{b}^{2}{c}^{2}-{a}^{2}{c}^{4}+{b}^{6}-{b}^{4}{c}^{2}-{b}^{2}{c}^{4} +{c}^{6}}{ \left( a+b+c \right) \left( a-b-c \right) \left( a-b+c \right) \left( a+b-c \right) }}={\frac {{x}^{4}{y}^{2}+{x}^{4}{z}^{2}+2\,{x}^{3}{y}^{3}+2\,{x}^{3}{y}^ {2}z+2\,{x}^{3}y{z}^{2}+2\,{x}^{3}{z}^{3}+{x}^{2}{y}^{4}+2\,{x}^{2}{y} ^{3}z+4\,{x}^{2}{y}^{2}{z}^{2}+2\,{x}^{2}y{z}^{3}+{x}^{2}{z}^{4}+2\,x{ y}^{3}{z}^{2}+2\,x{y}^{2}{z}^{3}+{y}^{4}{z}^{2}+2\,{y}^{3}{z}^{3}+{y}^ {2}{z}^{4}}{ \left( y+z+x \right) xyz}} $$ with $$a=y+z,b=x+z,c=x+y$$ the left-hand side is given by $$s(s-a)(s-b)(s-c)\left(\frac{1}{s^2}+\frac{1}{(s-a)^2}+\frac{1}{(s-b)^2}+\frac{1}{(s-c)^2}\right)$$ the result is the same as posted above