Dirichlet convolution of two arithmetic expressions (divisor-sum expansion plus symbolic identities)
DirichletConvolve[f, g, n, m]
DirichletConvolve[n, n, n, m]
→ m*DivisorSigma[0, m]DirichletConvolve[n + 2, n ^ 2 + 1, n, m]
→ 2*DivisorSigma[0, m] + DivisorSigma[1, m] + m*DivisorSigma[1, m] + 2*DivisorSigma[2, m]{DirichletConvolve[n, n ^ 2, n, m], DirichletConvolve[n ^ 2, n, n, m]}
→ {m*DivisorSigma[1, m], m*DivisorSigma[1, m]}