As I said, it is the chain rule. We have [itex]r= \sqrt{x^2+ y^2+ z^2}[/itex] so that [itex]\partial r/\partial x= x(x^2+ y^2+ z^2)^{-1/2}= x/r[/itex], [itex]\partial r/\partial y= y(x^2+ y^2+ z^2)^{-1/2}= y/r[/itex], [itex]\partial r/\partial z= z(x^2+ y^2+ z^2)^{-1/2}= z/r[/itex]. So [itex]grad r= (xi+ yj+ zk)/r[/itex].
If we write [itex]\vec{A(r)}= A_1(r)i+ A_2(r)j+ A_3(r)k[/itex] then [itex]d\vec{A(r)}/dr= (dA_1/dr) i+ (dA_2/dr)j+ (dA_3/dr)k[/itex] and [itex](d\vec{A})/dr\cdot grad r= (x(dA_1/dr)+ y(dA_2/dr)+ z(dA_3/dr))/r[/itex]
On the left, [itex]div \vec{A(r)}= dA_1/dr+ dA_2/dr+ dA_3/dr= [(\partial A_1/\partial x)(\partial r/\partial x)+ (\partial A_1/\partial y)(\partial r/\partial y)+ (\partial A_1/\partial z)(\partial z/\partial r)]+ [(\partial A_2/\partial x)(\partial r/\partial x)+ (\partial A_2/\partial y)(\partial r/\partial y)+ (\partial A_2/\partial z)(\partial z/\partial r)]+ [(\partial A_3/\partial x)(\partial r/\partial x)+ (\partial A_3/\partial y)(\partial r/\partial y)+ (\partial A_3/\partial z)(\partial z/\partial r)][/itex]
Now use the fact that [itex]\partial x/\partial r= 1/(\partial r/\partial x)= r/x[/itex], [itex]\partial y/\partial r= r/y[/itex], and [itex]\partial z/\partial r= r/z[/itex].