Solution:
This is equivalent to the existence of nonnegative integers c and d such that 3b−1=c(2a+1) and 3a−1=d(2b+1). Then
cd=(2a+1)(2b+1)(3b−1)(3a−1)=2a+13a−1⋅2b+13b−1<23⋅23=2.25.
Neither c nor d can equal 0 since that would give a=31 or b=31, so cd≤2.25 implies (c,d)∈{(1,1),(2,1),(1,2)}. Substituting (c,d) back in gives three systems of equations and the three solutions: (2,2),(12,17),(17,12).