I have come up with the following solution to this integral, but is just incomplete to my standards $$f(n)=\int_0^\infty \frac{1}{e^{x^n}+1}=\left(1-2^{(n-1)/n}\right )\zeta(n^{-1})\Gamma(1+n^{-1})$$
Seems to only work for $x\in\Bbb{N},x\gt 2$
This identity, therefore does not apply to $n=1$, and we all know that $f(1)=\ln 2$ because $\zeta(1)$ diverges.
So my question is this: How can you generalize the integral solution I gave to fit the case $n=1$?
In fact the formula does apply if you perform the limiting procedure correctly. Notice that
$$ 2^{(n-1)/n}-1\sim_{n\rightarrow 1}(n-1)\log(2)+\mathcal{O}(n-1)^2 $$
and
$$ -\zeta(n^{-1})\sim_{n\rightarrow 1} \frac{1}{n-1}+\mathcal{O}(1) $$
Therefore we obtain
which yields together with $\Gamma(2)=1$ the desired result!