DSolve[x^2 - 2*a*y[x] + (a + x*Derivative[1][y][x])^2 == 0,y[x],x]