Suppose we have two Dirichlet polynomials:
$$ f_1(s) = \sum_{n=1}^{m} \frac{a_n}{n^s} \\ f_2(s) = \sum_{n=1}^{m} \frac{b_n}{n^s} $$
Their product will also be a Dirichlet polynomial:
$$ f_1(s)f_2(s) = \sum_{n=1}^{m}\sum_{k=1}^{m} \frac{a_n b_k}{(nk)^s} = \sum_{n=1}^{m^2} \frac{c_n}{n^s}. $$
Is there a general way to express $c_n$ through $a_n$ and $b_n$? It appears that a variant of Dirichlet convolution can be used, that is,
$$ c_n = \sum_{d:\: d|n,\: d \leq m, \: n/d \leq m} a_d b_{\frac{n}{d}} $$
if $n \in \{lk : l=\overline{1,m}, k=\overline{1,m}\}$, and $c_n = 0$ otherwise, but such description looks too unwieldy and complex. Is there a simpler way?
Note: There is no simpler way compared to what you have already indicated. But, in fact it's not more complicated than the presumably more familiar Cauchy product of polynomials. Let's make a comparison head-to-head: