Hola
Bien. Tu idea se puede refinar para conseguir una muy buena cota.
1) Como bien dices basta maximizar la función \( f(x)=g(x)+g(1-x) \) siendo \( g(x)=x(1-x)^p \) en el intervalo \( [0,0.5) \).
2) Como bien has justificado:
\( g(x)\leq \dfrac{1}{p+1}\left(\dfrac{p}{1+p}\right)^p \) para todo \( x\in [0,0.5) \)
3) Ahora nos queda acotar \( h(x)=g(1-x)=(1-x)x^p \) en \( [0,0.5) \). Pero:
\( h'(x)=-x^p+p(1-x)x^{p-1}=x^{p-1}(p-px-x)\geq x^{p-1}(p-p/2-1/2)\geq 0 \)
Es decir la función es creciente en ese intervalo y así:
\( h(x)\leq h(0.5)=\dfrac{1}{2^{p+1}} \)
4) Combinando dos y tres obtenemos:
\( f(x)\leq \dfrac{1}{p+1}\left(\dfrac{p}{1+p}\right)^p+\dfrac{1}{2^{p+1}} \)
5) Es interesante observar que esa cota aunque más "fea" es mejor que la que yo proponía \( 1/(pe) \) para \( p\geq 8 \) y a media que crece \( p \)
se aproxima cada vez más al máximo exacto.
Esto se ve muy bien en la siguiente tabla: la primera columna es el valor de \( p \), la segunda el valor exacto del máximo, la tercer la cota obtenida en (4) y la cuarta la cota \( 1/(pe) \).

Saludos.