The infinitesimal differential forms $dx_1$ and $dx_2$ span a two-dimensional vector space (cotangent space). The transformation $f$ acts on $(x_1,x_2)$ like a coordinate transformation. Thus,
$$ dx_i \to du_i = \frac{\partial u_i}{\partial x_j}dx_j,\qquad i = 1,2.$$
Let's look at the infinitesimal transformation of $x_i$, i.e. $$ u_i(x) = x_i +\epsilon_i(x).$$
Show that if $f$ is a conformal transformation:
$$ \omega(x)\delta_{ij} = -\frac{\partial \epsilon_j}{\partial x^i}-\frac{\partial \epsilon_i}{\partial x^j},$$
where $\omega(x) \in \mathbb{R}$ is a scale factor.
Hint: First, show that $$\delta_{ij}du^idu^j=\left(1+\omega(x)\right)\delta_{kl}dx^kdx^l.$$
I've already tried different things such as substituting the $du$'s in the expression given as a hint. All of this basically led to nothing, so it would be great if somebody could help me find the solution. How should I start?
Maybe I should mention that this is an exercise from a physics workbook, so this question has to be answered using the information given above.