Recreational math: If $f(f(x))=e^x$, bound the integral $\int_0^1 f(x)dx$ I've been studying functions $f:\mathbb R\to\mathbb R$ that satisfy $f(f(x))=e^x$ (or, half-iterates of the exponential function). I know that there's only one such analytic function, but it's really hard to study since it is almost certainly non-elementary and I only know how to find finitely many terms of its Maclaurin Series. 
Instead, I'm studying all continuous and increasing functions $f$ satisfying $f(f(x))=e^x$, and I've alighted on the following problem (which I came up with out of curiosity). I propose this question to all interested residents of MSE:

Given that $f$ is continuous and increasing and $f(f(x))=e^x$, find some bounds for the integral
  $$\int_0^1 f(x)dx$$

I've managed to come up with some pretty sweet bounds (in fact, they are the best possible bounds), and I'll post them after this question gets some answers.
I'll accept whatever answer has the tightest bounds, with a proof.
NOTE: Most people probably wouldn't think of this as recreational "fun" math, but hey, I did it for fun, and I'm proposing it as a problem to be done just for fun. So please try to enjoy it, and please don't try to close it.
 A: IMPORTANT EDIT:
I find that $f(x)$ can create some paradoxical result.
Assume $0<a<1$(which is the region of interest) and $f^{-1}(a)=b$. Then $0<b<a$.
Then, $$a=f(b)\implies \ln a=\ln f(b)=f^{-1}(b)$$
Clearly, $L.H.S.\le0$. However, we have proved that $f^{-1}(b)=R.H.S.>0$.
Can anyone resolve it?
END EDIT
Not a tight bound, but too long for a comment.
Assume the $f(x)$ is:


*

*Continuous

*A real function on the interval $(-\infty,\infty)$



*Strictly increasing



Since $f(x)$ is strictly increasing, it has a well-defined inverse $f^{-1}(x)$.
Let’s redefine the function as
$$f(x)=e^{f^{-1}(x)}$$
Firstly, since $e^x$ is positive for all real $x$, hence $f(x)$ is strictly positive.
Secondly, we are interested in the fixed point of $f(x)$. For fixed points $(t,f(t))$, $$t=f(t)= f^{-1}(t)$$
Thus, $$t=e^t$$ which has no real solution.
Therefore, we can see $y=f(x)$ has no intersection with $y=x$.
With these two observations, one can deduce $f(x)>x$ for all $x$.
$$f(x)>x\implies \int^1_0f(x)dx>\int^1_0xdx=\frac12$$
EDIT
From $$f(x)>x$$
we have $$f(f(x))>f(x)$$
$$e^x>f(x)$$
$$e-1>\int^1_0f(x)dx$$
EDIT 2:
I think that $f(0)$ cannot be defined(but maybe okay in the sense of limit).
Assume $f^{-1}(0)=C$.
It is easy to prove $f^{-1}(x)$ is strictly increasing.
Thus, from $f(x)>x$, we get $x>f^{-1}(x)$. So $C<0$.
We also observe that, due to $f(x)>0$, $f^{-1}(x)$ is undefined for $x<0$. Since $f(x)$ is defined wherever its inverse is defined, so $f(x)$ is undefined in $(0,-\infty)$.
$$f^{-1}(0)=C\implies f(C)=0$$
But $C<0$, $f(C)$ should be undefined! Also, I have proven that $f(C)>0$(above). That leads to a contradiction.
So a value of $C$ does not exist, and hence $f^{-1}(0)$ and $f(0)$ is undefined.
From the above observations, we can also see: $$0<f^{-1}(x)<x$$.
A: Suppose $f(x)$ is an increasing, continuous function satisfying $f(f(x))=e^x$. Then there must be some $\alpha\in(0,1)$ such that $f(\alpha)=1$. 
We can compute
$$
\begin{align*}
\int_\alpha^1 f(x)\,dx&=\int_\alpha^1 e^{f^{-1}(x)}\,dx\\
&=\int_0^\alpha e^u\,df(u)\\
&=e^\alpha f(\alpha)-e^0f(0)-\int_0^\alpha f(u)e^u\,du\\
&=e^\alpha-\alpha-\int_0^\alpha f(x)e^x\,dx.
\end{align*}
$$
So, we have
$$
\int_0^1 f(x)\,dx=\left(\int_0^\alpha+\int_\alpha^1\right) f(x)\,dx=e^\alpha-\alpha-\int_0^\alpha(e^x-1)f(x)\,dx.
$$
For $0< x< \alpha$, we have $\alpha< f(x)< 1$, and thus
$$
1=e^\alpha-\alpha-\int_0^\alpha e^x\,dx<\int_0^1 f(x)\,dx<e^\alpha-\alpha-\int_0^\alpha \alpha e^x\,dx=(1-\alpha)e^\alpha+\alpha^2.
$$
For $0< \alpha< 1$, the quantity $(1-\alpha)e^\alpha+\alpha^2$ takes a maximum value of $2+\log^2(2)-\log(4)$ at $\alpha=\log(2)$, so we have bounds
$$
1<\int_0^1 f(x)\,dx < 2+\log^2(2)-\log(4)\approx 1.0942.
$$
These bounds are sharp: for $0\leq x\leq \log(2)$, we can take $f(x)$ to be an arbitrary continuous, strictly increasing function satisfying $f(0)=\log(2)$ and $f(\log(2))=1$, then extend $f$ to a continuous increasing function on the entire real line satisfying $f(f(x))=e^x$ by iteratively using the relation $f(x)=e^{f^{-1}(x)}$. To make $\int_0^1 f(x)\,dx$ arbitrarily close to the lower bound, we can choose $f(x)$ to be close to $1$ for an arbitrarily large portion of the interval $[0,\log(2)]$, and to make the integral arbitrarily close to the upper bound, we can choose $f(x)$ to be close to $\log(2)$ for an arbitrarily large portion of the interval $[0,\log(2)]$.
