Let me state my motivation. I believe that the Polya Enumeration Theorem and Multiple Zeta Values (the classic being the Basel problem and the values of the Riemann zeta function at the even integers) are among the most intriguing objects in mathematics.
With this message I present a conjecture that connects these two (PET and MZVs). Most likely there is a proof somewhere in the literature, which I invite readers to submit in text form if it admits a compact presentation or as a reference if it does not.
To understand this conjecture you need to learn about multiset cycle indices, which I presented in a different context at this MSE link.
Multisets and their cycle indices form a combinatorial species like any other, e.g. cycles, sequences and sets. They represent multisets and are identified by the partition that corresponds to the multiplicities of the elements of the multiset, e.g. $\mathfrak{M}_{1,2,3}$ is a multiset that contains three elements, one of which in two copies and another in three. Substitution of an OGF into a multiset cycle index yields the generating function of the multiset.
There is a special multiset cycle index which is written $\mathfrak{M}_{1,1,1,\ldots,1}$ which gives the multiset where all elements are unique, i.e. the species of sets. This cycle index has been known for quite some time and is used to compute the set operator $\mathfrak{P}.$ The cycle index $Z(P_n)$ of the set operator $\mathfrak{P}_{=n}$ is the difference between the cycle index $Z(A_n)$ of the alternating group and the cycle index $Z(S_n)$ of the symmetric group on $n$ elements. It admits the following simple recursive definition:
$$Z(P_0) = 1 \quad\text{and}\quad Z(P_n) = \frac{1}{n} \sum_{l=1}^n (-1)^{l+1} a_l Z(P_{n-l}).$$
Here are the cycle indices $Z(P_3), Z(P_4)$ and $Z(P_5):$ $$\begin{array}{|l|l|} \hline Z(P_3) & \frac{1}{6}\,{a_{{1}}}^{3}-1/2\,a_{{2}}a_{{1}}+1/3\,a_{{3}}\\ \hline Z(P_4) & \frac{1}{24}\,{a_{{1}}}^{4}-1/4\,a_{{2}}{a_{{1}}}^{2} +1/3\,a_{{3}}a_{{1}}+1/8\,{a_{{2}}}^{2}-1/4\,a_{{ 4}}\\ \hline Z(P_5) & {\frac {1}{120}}\,{a_{{1}}}^{5}-\frac{1}{12}\,a_{{2}}{a_{{1}}}^{3} +1/6\,a_{{3}}{a_{{1}}}^{2}+1/8\,a_{{ 1}}{a_{{2}}}^{2}-1/4\,a_{{4}}a_{{1}} -1/6\,a_{{2}}a_{{3}}+1/5\,a_{{5}}\\ \hline\end{array}$$
With these definitions we are now ready to state the conjecture which is quite simply that $$\large \color{blue}{\zeta(s, s, s, \ldots, s) = Z(P_n)(\zeta(s), \zeta(2s), \ldots, \zeta(ns))}$$ i.e. the MZV of a unique argument $s$ repeated $n$ times is equal to the substituted cycle index $Z(P_n)$ of the set operator $\mathfrak{P}_{=n}$ according to the rule $a_l = \zeta(ls).$
Note that the case $n=3$ is found on the Wikipedia page for MZVs that I linked to in the introduction.
Example. The conjecture says e.g. that $$\zeta(5,5,5,5) = 1/24\, \zeta \left( 5 \right) ^{4}-1/4\,\zeta \left( 10 \right) \zeta \left( 5 \right) ^{2}+1/3\,\zeta \left( 15 \right) \zeta \left( 5 \right) +1/8\, \zeta \left( 10 \right) ^{2}-1/4\,\zeta \left( 20 \right).$$
Important observation. The reader may well ponder the statement of the conjecture and say that there is nothing to prove here, which I will accept as proof, if accompanied by a brief explanation why this should be so.
Addendum. An interesting application of the above result is found at this MSE link.