here $d(n)$ counts the number of positive divisors of $n$.
I've tried 2 things:
Using Bell series. But then again it just showed me that the bell series of the square of a function is not the square of the bell series of the function.
Expand it all and use combinatorics. If: $$n = \prod_{i=1}^r p_i^{\alpha_i} $$
then if $d \vert n$:
$$d = \prod_{i=1}^rp_i^{\beta_i}, 0\leq \beta_i \leq \alpha_i$$
so
$$\sum_{t \vert n}d(t) = \sum \prod_{i=1}^r (\beta_i + 1)$$ for all possible combitions of $(\beta_1,\beta_2,...,\beta_r)$ but I couldn't develop it very much