I am trying to solve the system of equations. My code

```
{x, y, z, a, b, m, n} /. Solve[{1/x + 1/y == 1/z, a/x + b/y == m/n, 2 <= x <= 30, 2 <= y <= 30, 1 <= z <= 30, 2 <= a <= 10, 2 <= b <= 10, a > b, 2 <= m <= 10, 2 <= n <= 10, x > y, GCD[m, n] == 1, GCD[a, b, m] == 1, m <= n}, {x, y, z, a, b, m, n}, Integers]
```

For a long time, I can not get the result. How can I reduce the time?