(The following is inspired by Integral inequality with a function twice differentiable: Integrating by parts transforms the integral with $f$ to an integral with $f'$. The condition $\int_0^1 f(x) \, dx = 0$ is used to add a term to the first integral so that the $u(b)v(b)-u(a)v(a)$ term vanishes. Cauchy-Schwarz then helps to estimate the integral containing $f'$ by an integral containing $f'^2$.)
Integrating by parts we get
$$
\int_0^1 xf(x) \, dx = \frac 12 \int_0^1 (2x-1)f(x) \, dx = \frac 12 \int_0^1 x(1-x) f'(x) \,dx \\
= \frac 12 \int_0^1 \frac{x \sqrt{1-x}}{\sqrt{1+x}} \sqrt{1-x^2} f'(x) \, dx \, .
$$
Now apply Cauchy-Schwarz:
$$
\left( \int_0^1 xf(x) \, dx \right)^2 \le \frac 14 \int_0^1 \frac{x^2(1-x)}{1+x} \, dx \int_0^1 (1-x^2) (f'(x))^2 \, dx \\
\le \frac 14 \int_0^1 x^2(1-x) \, dx \int_0^1 (1-x^2) (f'(x))^2 \, dx \\
= \frac{1}{48 }\int_0^1 (1-x^2) (f'(x))^2 \, dx
$$
which is better than the desired estimate by a factor of $2$.
Using the exact value $\int_0^1 \frac{x^2(1-x)}{1+x}\, dx = 2 \ln(2) - 4/3$ we get the sharp estimate
$$
\int_0^1 (1-x^2) (f'(x))^2 \, dx \ge C \left( \int_0^1 xf(x) \, dx \right)^2
$$
with
$$
C = \frac{2}{\ln(2)-2/3} \approx 75.53 \, .
$$
Equality holds if equality holds in the Cauchy-Schwarz inequality, and that is if
$$
f'(x) = \text{const} \cdot \frac{x}{x+1}
$$
so that the integrands are linearly dependent. Together with the condition $\int_0^1 f(x) \, dx = 0$ this gives (up to a multiplicative constant)
$$
f(x) = x - \ln(x) + 2 \ln(2) - \frac 23 \, .
$$