First suppose that a=1. Then φ(1)=1. For all positive integers b we have gcd(a,b)=1. Therefore in this case the equation is 1+b=1+φ(b)+1, or equivalently, φ(b)=b−1. This is equivalent to the statement that there exists a unique integer from {1,2,…,b}; which then has to be b itself (since unless b=1, we have gcd(b,b)>1, but if b=1 we have φ(b)=b). In other words, this is equivalent to b being a prime number. Hence the solutions for a=1 are precisely the pairs (1,p) with p a prime number. Similarly, the solutions for b=1 are precisely the pairs (p,1) with p a prime number.
Now assume that a,b≥2. As gcd(b,b)>1 we have φ(b)≤b−1. Therefore
gcd(a,b)=a+b−φ(a)−φ(b)≥a−φ(a)+1.
Let p be the minimal prime divisor of a (which exists as a≥2). Since for all multiples tp≤a of p, we have gcd(tp,a)>1, it follows that a−φ(a)≥p1⋅a. Therefore we have
gcd(a,b)≥a−φ(a)+1≥pa+1.
The two largest divisors of a are a and pa. Since gcd(a,b) is a divisor of a that is at least pa+1, it must equal a. Hence gcd(a,b)=a. In the same way we prove that gcd(a,b)=b. So a=b.
The equation now is equivalent to 2a=2φ(a)+a, so also to a=2φ(a). Note that 2∣a. Therefore write a=2k⋅m with k≥1 and m odd. By a well-known property of the φ-function, we have φ(a)=φ(2k)⋅φ(m)=2k−1⋅φ(m), and the equation becomes 2k⋅m=2⋅2k−1⋅φ(m), or equivalently
m=φ(m).
Therefore m=1, and a=b=2k. Indeed, the equation holds for all pairs (2k,2k) with k≥1.