The original sum can be written as
$$T(\alpha,\beta,\gamma,n)=\sum_{x,y\le n}x^\alpha y^\beta(x,y)^\gamma,$$
where $(x,y)=\mathrm{gcd}(x,y)$. One can find asymptotic formula for this sum using standart approach. Let $d=(x,y)$. Then
$$T(\alpha,\beta,\gamma,n)=\sum_{d\le n}d^\gamma\sum_{{x,y\le n\atop (x,y)=d}}x^\alpha y^\beta=\sum_{d\le n}d^{\alpha+\beta+\gamma}\sum_{{x,y\le n/d\atop (x,y)=1}}x^\alpha y^\beta.$$
The condition $(x,y)=1$ can be removed using Möbius function:
$$T(\alpha,\beta,\gamma,n)=\sum_{d\le n}d^{\alpha+\beta+\gamma}\sum_{\delta\le n/d}\mu(\delta)\sum_{{x,y\le n/d\atop \delta\mid(x,y)}}x^\alpha y^\beta=\sum_{d\le n}d^{\alpha+\beta+\gamma}\sum_{\delta\le n/d}\mu(\delta)\delta^{\alpha+\beta}\sum_{x,y\le n/(d\delta)}x^\alpha y^\beta.$$
The last sum (for $\alpha,\beta>-1$) is $\sim\frac{n^{\alpha+\beta+2}}{(\alpha+1)(\beta+1)(d\delta)^{\alpha+\beta+2}},$ so (for $\gamma>1$)
$$T(\alpha,\beta,\gamma,n)\sim \frac{n^{\alpha+\beta+\gamma+1}}{\zeta(2)(\alpha+1)(\beta+1)}.$$

The special case $\gamma=1$,  $\alpha=0$, $\beta=−1$ is more tricky. We can write the given sum as
$$T(n)=\sum_{x,y\le n}\frac{(x,y)}{y}=T_1(n)+T_2(n),$$
where for some $U>1$
$$T_1(n)=\sum_{{x,y\le n \atop (x,y)\le U}}\frac{(x,y)}{y},\quad T_2(n)=\sum_{{x,y\le n \atop (x,y)> U}}\frac{(x,y)}{y}.$$ The second sum will be in the error term ($y=dy_1$, $x=dx_1$):
$$T_2(n)=\sum_{d> U}d\sum_{{x,y\le n \atop (x,y)=d}}\frac{1}{y}\ll \sum_{d> U}\sum_{x_1,y_1\le n/d }\frac{1}{y_1}\ll \sum_{d> U}\frac{n}{d}\log\frac{n}{d}\ll n\log^2 \frac{n}{U}.$$
Here it is clear that for $U=n\log^{-2}n$ we get error term $O(R(n))$ with $R(n)=n\log^2\log n.$

The first sum gives main term:
$$T_1(n)=\sum_{y\le n}\frac{1}{y}\sum_{{d\le U \atop d\mid y}}d\sum_{{x\le n \atop (x,y)=d}}1=\sum_{y\le n}\frac{1}{y}\sum_{{d\le U \atop d\mid y}}d\sum_{{x_1\le n/d \atop (x_1,y/d)=1}}1.$$ The last sum is [known][1]:
$$\sum_{{x_1\le n/d \atop (x_1,y/d)=1}}1=\frac{\varphi(y/d)}{y/d}\frac nd+O(\tau(y/d)).$$
Hence
$$T_1(n)=n\sum_{y\le n}\frac{1}{y}\sum_{{d\le U \atop d\mid y}}\frac{\varphi(y/d)}{y/d}+O(R_1(n)),$$
where
$$R_1(n)=\sum_{d\le U}d\sum_{{y\le n \atop d\mid y}}\frac{\tau(y/d)}{y}=\sum_{d\le U}\sum_{y_1\le n/d }\frac{\tau(y_1)}{y_1}\ll U\log^2n\ll R(n).$$
So $$T_1(n)=n\sum_{d\le U}\frac{1}{d}\sum_{y_1\le n/d}\frac{\varphi(y_1)}{y_1^2}+O(R(n))=n\sum_{d\le U}\frac{1}{d}\left(\frac{1}{\zeta(2)}\left(\log (n/d)+\gamma-\frac{\zeta'(2)}{\zeta(2)}\right)+O\left(\frac{\log n}{n/d}\right)\right)+O(R(n))=\frac{n}{\zeta(2)}\sum_{d\le U}\frac{1}{d}\left(\log (n/d)+\gamma-\frac{\zeta'(2)}{\zeta(2)}\right)+O(R(n))=\frac{n}{\zeta(2)}\left(\left(\log n+\gamma-\frac{\zeta'(2)}{\zeta(2)}\right)\left(\log U+\gamma+O(1/U)\right)-\sum_{d\le U}\frac{\log d}{d}\right)+O(R(n)).$$
We also  [know][2] that
$$\sum_{d\le U}\frac{\log d}{d}=\frac{\log^2 U}{2}+\gamma_1+O(U^{-1}\log U).$$
Collecting all together we'll have
$$T(n)=\frac{n}{\zeta(2)}\left(\frac{\log^2n}{2}+\log n\left(2\gamma-\frac{\zeta'(2)}{\zeta(2)}\right)\right)+O(n\log^2\log n).$$


  [1]: https://mathoverflow.net/q/22953/5712
  [2]: https://en.wikipedia.org/wiki/Stieltjes_constants