nn = 32;
b = Table[
Exp[MangoldtLambda[Divisors[n]]]^-MoebiusMu[Divisors[n]], {n, 1,
nn}];
j = 1;
MatrixForm[
Table[Table[
Product[(b[[n]][[m]]*
DirichletCharacter[b[[n]][[m]] + 2, j + 1, k] - (b[[n]][[m]] -
1)), {m, 1, Length[Divisors[n]]}], {n, 1, 32}], {k, 1, 32}]]
Above is as a Mathematica program a generalization of the Dirichlet inverse of the Euler totient function to Dirichlet characters.