I'm trying to simplify this combinatorial looking sum: $$\sum_{ax+by=n}_{(a,x,b,y)\in \mathbb{N^4}}\max{\{a,b\}}$$ In terms of possibly some scaled divisor functions plus a Cauchy product/convolution of divisor functions.
I thought maybe I could break it up into the cases where $(a<b)\wedge (a>b) \wedge (a=b)$
Then separate the cases and use the representation:
$$\sum_{ax+by=n}_{(a,x,b,y)\in \mathbb{N^4}}a=\sum_{ax+by=n}_{(a,x,b,y)\in \mathbb{N^4}}b=\sum_{k=1}^{n-1}d(k)\sigma(n-k)$$
To evaluate the entire sum, but it's not really working out.
Can someone more experienced with this sort of thing help me simplify the first sum in terms of divisor function convolutions?
Or if it doesn't look like it can be simplified in such a closed form, give me some insight as to why it can't?