From this equation, $a·f_x + b·f_y = c$ we can get this set of coupled ODE's:
$$\frac{dx}{a}=\frac{dy}{b}=\frac{df}{c}$$
In your case:
$$\frac{dx}{y}=\frac{dy}{x}=\frac{df}{1}$$
Considering that whenever $\frac{m}{n}=\frac{o}{p}$ then $\frac{m+o}{n+p}=\frac{m}{n}$, we get:
$$\frac{dx+dy}{x+y}=\frac{df}{1};\;u=x+y;\;du=dx+dy$$
$$\frac{du}{u}=\frac{df}{1};\;\ln{u} = f + c_1;\;\ln{(x+y)}=f+c_1$$
Now, with $\frac{dx}{y}=\frac{dy}{x}$, $xdx=ydy$, so:
$$x^2=y^2+c_2\;or\;x^2-y^2=c_2$$
We have now the equations of two surfaces, the intersection of which form the characteristics. It must exist a relation beetween them. The equality $c_1=g(c_2)$ is such a relation, with g an arbitrary single variable differentiable function. Then, $c_1=g(x^2-y^2)$
At last, we get the general solution:
$$\ln{(x+y)}=f+g(x^2-y^2)\;or$$
$$f(x,y)=\ln{(x+y)}+g(x^2-y^2)$$
We are ready to impose the boundary conditions. We need at least the values of f along some curve in order to get an unique solution. Unfortunately yours only offers the value k at a single point $(x_0,y_0)$, so we have a whole set of surfaces solution of the PDE.