Let $0<x<1$ then we have :
$$x^{{\operatorname{W}(2ex)}^{2x}}+(1-x)^{{\operatorname{W}(2e(1-x))}^{2(1-x)}}\leq 1$$
The equality case is $x=0.5$.
To show it I have tried to follow the lemma 7.1 and 7.2 of this paper by Vasile Cirtoaje.The problem is that the resulting expression is awful !
I have tried also Bernoulli's inequality with any effect because it's not sharp enough.
Update 18/12/2020 :
It's a another tried .We can build an approximation like this :
Let $0<\beta<x\leq 0.5$ then we have to determine the constants such that:
$$x^{{\operatorname{W}(2ex)}^{2x}}\leq \frac{1}{2}\operatorname{W}^{\alpha}(2ex)$$
We have numerically speaking $\frac{115}{100}<\alpha<\frac{125}{100}$
To reduces the gap I have tried to introduce a linear function : $$x^{{\operatorname{W}(2ex)}^{2x}}\leq \frac{1}{2}\operatorname{W}^{\alpha}(2ex)+ax+b$$
But again it's not enough to work so we may consider a general polynomial like :
$$x^{{\operatorname{W}(2ex)}^{2x}}\leq \frac{1}{2}\operatorname{W}^{\alpha}(2ex)+\sum_{k=0}^{n}a_nx^n$$
Well it's a first step and in the future I shall tried to find the coefficients of this general polynomial .
Update 20/12/2020 :
We can reformulate the problem as :
Let $x,y>0$ such that $ye^y+xe^x=2e$ then we have :
$$\left(\frac{xe^x}{2e}\right)^{(x)^{\frac{xe^x}{e}}}+\left(\frac{ye^y}{2e}\right)^{(y)^{\frac{ye^y}{e}}}\leq 1$$
Where I use the inverse function of the Lambert's function .
Well using the form $f(x)=\left(\frac{xe^x}{2e}\right)^{(x)^{\frac{xe^x}{e}}}=g(x)^{h(x)}$ I can show that the function $f(x)$ is convex on $(0,W(2e))$ so (I have tried) we can use Slater's inequality to find an upper bound .Like this it doesn't work . On the other hand we can use Karamata's inequality but I have not tried !
Well If we use Karamata's I have a strategy :
We have by Karamata's inequality and $0\leq\varepsilon_n'\leq\varepsilon_n<y<x$:
$$f(x)+f(y)\leq f(x+\varepsilon_n)+f(y-\varepsilon_n')$$
With $(y-\varepsilon_n')e^{y-\varepsilon_n'}+(x+\varepsilon_n)e^{x+\varepsilon_n}\geq 2e$
Now we want to repeat the process to get a series of inequalities of the kind :
$$f(x)+f(y)\leq f(x+\varepsilon_n)+f(y-\varepsilon_n')\leq f(x+\varepsilon_{n-1})+f(y-\varepsilon_{n-1}')< 1$$
But it's very complicated.
It doesn't work for all the value but I think we have the inequality $y> 0.5 \geq x$ :
$$p(x)=(1-x^{xe^{x-1}})^2+x^{xe^{x-1}} \frac{xe^{x-1}}{2} (2-x^{xe^{x-1}})-x^{xe^{x-1}} \frac{xe^{x-1}}{2} (1-x^{xe^{x-1}}) \ln\left(\frac{xe^{x-1}}{2}\right)$$ We have : $$f(x)+f(y)\leq p(y)+2^{-\varepsilon}p^{1+\varepsilon}(x)< 1$$
With $0\leq \varepsilon \leq\frac{1}{10}$
Where I use the Lemma 7.2 of the paper above .
The last idea :
Using the majorization theorem :
Let $a\geq b>0$ and $c\geq d >0$ and $n$ a natural number large enough such that :
$$a\geq c$$
And :
$$\left(a\frac{n}{n+1}+c\frac{1}{n+1}\right)\left(b\frac{n}{n+1}+d\frac{1}{n+1}\right)\geq cd$$
Then we have :
$$a+b\geq c+d$$
Proof:It's a direct consequence of the Karamata's inequality.
We have another theorem :
Let $2>x,y>0$ ,$n$ a natural number large enough and $\varepsilon>0 $
If we have :
$$xy<1-\varepsilon $$ $$x+y<2-\varepsilon$$ then we have :
$$\ln\left(\frac{n}{n+1}+x\frac{1}{n+1}\right)+\ln\left(\frac{n}{n+1}+y\frac{1}{n+1}\right)\leq 0$$
Example :
Using the theorem of majorization we have ($x=0.4$):
$$(1-x)^{{\operatorname{W}(2e(1-x))}^{2(1-x)}}< 1-\operatorname{W}^{1.25}(2ex)0.5$$
And :
$$\left(\frac{1}{4000}x^{{\operatorname{W}(2ex)}^{2x}}+\frac{3999}{4000}\operatorname{W}^{1.25}(2ex)0.5\right)\left(\frac{1}{4000}(1-x)^{{\operatorname{W}(2e(1-x))}^{2(1-x)}}+\frac{3999}{4000}(1-\operatorname{W}^{1.25}(2ex)0.5)\right)< (1-\operatorname{W}^{1.25}(2ex)0.5)\operatorname{W}^{1.25}(2ex)0.5$$
Dividing both side by the RHS and using the second theorem remarking that :
$$\frac{x^{{\operatorname{W}(2ex)}^{2x}}(1-x)^{{\operatorname{W}(2e(1-x))}^{2(1-x)}}}{\operatorname{W}^{1.25}(2ex)0.5(1-\operatorname{W}^{1.25}(2ex)0.5)}<1-\varepsilon$$
And :
$$\frac{x^{{\operatorname{W}(2ex)}^{2x}}}{\operatorname{W}^{1.25}(2ex)0.5}+\frac{(1-x)^{{\operatorname{W}(2e(1-x))}^{2(1-x)}}}{(1-\operatorname{W}^{1.25}(2ex)0.5)}<2-\varepsilon\quad (I)$$
Now I think it's easier because we can take the logarithm and study the behavior of the function .
To prove the $(I)$ we can use the bound :
Let $0<x<\frac{1}{100}$ :
$$e^x<(1+x)^2-x$$
Obviously if we study separately the differents elements of the LHS .
Then to study $(I)$ we have a quite good approximation :
Let $0< x \leq \frac{1}{2}$ then we have :
$${\operatorname{W}(2ex)}^{2x}\geq (2x)^{\frac{915}{1000}\left(x\right)^{\left(\frac{87}{100}\right)}}$$
In fact we have the following refinement on $(0,0.5]$ :
$$x^{{\operatorname{W}(2ex)}^{2x}}+(1-x)^{{\operatorname{W}(2e(1-x))}^{2(1-x)}}\leq x^{(2x)^{\frac{915}{1000}\left(x\right)^{\left(\frac{87}{100}\right)}}}+ (1-x)^{(2(1-x))^{\frac{915}{1000}\left((1-x)\right)^{\left(\frac{87}{100}\right)}}}\leq 1$$
Remarks : The method using the majorization theorem have two advantages . We need to choose two values of the same order with respect to the values in the LHS . One can be inferior (and the other necessary superior).On the other hand the bound with the exponential ,his accuracy depends of the initial approximation in $(I)$ . Finally if we split in two the LHS in $(I)$ and if for a one we prove a stronger result then the other element is slighlty easier to show .
I build an approximation on $(0,1)$ wich have the form :
$$x^{(2x)^{\frac{915}{1000}\left(x\right)^{\left(\frac{87}{100}\right)}}}\simeq \left(\left(2^{(2x)^{x^{1.25}}} \frac{x}{2}\right)^{0.5}0.5^{0.5}*0.5^{{(2 (1-x))}^{x^{-0.25}}}\right)^{0.5}\quad (S)$$
You can play with the coefficients $-0.25$ and $1.25$ wich are not the best (make me a comment if you have better please:-))
We can slightly improve $(S)$ in using the logarithm we have on $[0.5,1)$:
$$x^{(2x)^{\frac{915}{1000}\left(x\right)^{\left(\frac{87}{100}\right)}}}\simeq \left(\left(2^{(2x)^{x^{3}}} \frac{x}{2}\right)^{0.5}0.5^{0.5}*0.5^{{(2 (1-x))}^{x^{-0.2}}}\right)^{0.5}-0.5\ln\left(\left(\left(2^{(2x)^{x^{3}}} \frac{x}{2}\right)^{0.5}0.5^{0.5}*0.5^{{(2 (1-x))}^{x^{-0.2}}}\right)^{0.5}\right)+0.5\ln\left(x^{(2x)^{\frac{915}{1000}\left(x\right)^{\left(\frac{87}{100}\right)}}}\right)\quad (S')$$
We can replace the coefficient $\frac{915}{1000}$ by $\frac{912}{1000}$,and $3$ by $3.5$ and finally $-0.2$ by $-0.19$ and I think it's the same order so we can apply the majorization theorem .Ouf!
Any idea to solve it ?
Thanks