For every positive integer $d$, we let $\tau\left(d\right)$ be the number of positive divisors of $d$.
Prove that \begin{align} \sum_{d|n} \tau^3(d) = \left(\sum_{d|n} \tau (d)\right)^2 \end{align} for each positive integer $n$, where the sums range over all positive divisors $d$ of $n$.
Now I only know that both sides are multiplicative arithmetic functions in $n$. Could you tell me what I need to do next?