En primer lugar,
$$I=\int_0^1x(1-x)f(x)f'(x)\text{d}x \,\overset{IBP}{=\!=\!=}\,\frac12x(1-x)f^2(x)\bigg|_0^1+\frac12\int_0^1 (2x-1)f^2(x)\text{ d}x,$$
$$I=\frac12\int_0^1 (2x-1)f^2(x)\text{ d}x.\tag1$$
Las condiciones dadas significa que $f(x)$ tiene una raíz de multiplicidad $2$
en $x=0$ y una raíz simple de a $x=1.$
Teniendo en cuenta la Besou y teorema de la condición de $|f'''(x)|\le 1,$ uno puede conseguir
$$|f(x)| \le a(x-0)^2(1-x) = a(x^2-x^3),\quad a\in\left(0,\frac16\right),\tag2$$
$$|I|\le\dfrac1{72}\int_0^1(2x-1)(x^2-x^3)^2\text{ d}x = \dfrac1{30240}.$$
$\color{green}{\textbf{EDIT of 12.08.20.}}$
Permítanos integrar el dado de la desigualdad
$$-1\le f'''(x) \le 1\tag{1n}$$
bajo las condiciones dadas en el intervalo de $(0,x),$luego
$$
\begin{cases}
-x\le f''(x) - f''(0)\le x\\[4pt]
-\frac12x^2 \le f'(x) - xf''\left(0\right) \le \frac12x^2\\[4pt]
-\frac16x^3 \le f(x) - \frac12x^2f''\left(0\right) \le \frac16x^3._{\Large\mathstrut}
\end{casos}\etiqueta{2n}$$
De $(2n.3)$debe
$$|6f(x)-3f''(0)x^2| \le x^3,\quad -1\le -3f''(0)\le1,$$
$$6|f(x)|\le |x^3+ax^2|,\quad |a|\le 1,\quad f(1)=0,$$
$$6|f(x)|\le x^2-x^3 = h(x).$$
Esto confirma las fórmulas de $(2).$
Del mismo modo, de $(2n.2)$
$$|f'(x)| \le \frac12 x^2 -\frac13x = h'(x),$$
en $h'(x)$ está sincronizado con $h(x).$
Por lo tanto,
$$|I| \le \left|\int_0^1 x(1-x) h(x) h'(x) \text{ d}x\right| = \frac1{30240}.$$