I understand the general Diophantine Equation proof when $ax + by = d,$ where $d = gcd(a,b),$ but I am stuck on how the proof changes for a multiple of the gcd
It is easy to transform (scale) one solution to the other, after recalling their linear genesis.
Recall that the general solution of a linear equation $\,L(x,y) := ax+by = d\,$ is the sum of any particular solution $\,(x_0,y_0)$ plus the general solution of the associated $\rm\color{#c00}{homogeneous}$ equation $\,L(x,y) = \color{#c00}0,\,$ which here is $\,(x,y) = k(b,-a)/d = (kb/d,-ka/d)$.
Now $L(x,y) = md\,$ has the same associated homogeneous equation so same homogeneous solution, and scaling yields a particular solution: $\, md = mL(x_0,y_0) = L(m(x_0,y_0)),\,$ so $m(x_0,y_0)$ is a particular solution. Therefore $\,L(x,y) = md\,$ has the general solution
$$(x,y)\, =\, \underbrace{m(x_0,y_0)}_{\rm particular} + \underbrace{k(b,-a)/d}_{\rm homogeneous}\, =\, (mx_0+kb/d,\, my_0-ka/d)\qquad $$
Summarizing $ $ the solution transformation process in general:
If $\,L\,\bar x_0 = d\ $ & $\ [L\,\bar x = 0 \!\iff\! \bar x = \bar x_h]\,$ then $\,L\,\bar x = d\,$ has general solution $\ \bar x\, =\, \bar x_0\, +\,\bar x_h$
so $\, L\,m\bar x_0 = md,\,$ therefore we deduce that: $\, L\bar x = md\,$ has general solution $\,\bar x = m\bar x_0 + \bar x_h$
e.g. from calculus: a simple linear differential equation
$\ D f\, =\, x\iff f = \frac{1}2 x^2 + c\,\ $ ($f= c\,$ is general solution of homogenous form $\,D f = 0)$
$D f = m x\!\iff\! f = \frac{m}2 x^2 + c,\,$ by simply scaling the particular solution by $\,m,\,$ as above, i.e. $\,\ D_x^{-1}(mx)\, =\, m\:\! D_x^{-1}(x),\ $ i.e. $\ \int_x mx = m\int_x x$.