Here's a less trivial solution. Without loss of generality, translate your integration variable so the bounds of the integral are $[-a,a]$, where $a=(x_1-x_0)/2$.
$m$ and $n$ have opposite parity:
When $m$ and $n$ have opposite parity, a solution of the form $f(x,y) = \sum_{i} \alpha_i(x) \beta_i(y)$, where each of the $\alpha_i(x)$ is an odd function is a solution. This means that both $f(x,y)$ and $\frac{\partial f}{\partial y}$, are odd functions of $x$ and that consequently $f^{-n} (\frac{\partial f}{\partial y})^m$ is as well, meaning the integral will vanish for any value of $y$.
$m$ and $n$ are both even:
In this case $f^n$ and $(\frac{\partial f}{\partial y})^m$ are both positive, which requires that either $f(x,y)=0$ or $\frac{\partial f(x,y)}{\partial y}=0$ for each value of $(x, y)$. I think this pretty much means that the only solutions are trivial, ie $f(x,y)=g(x)$ for some function $g(x)$.
$m$ and $n$ are both odd:
Need to think about this a bit more.