Suppose that a polynomial $p(x_1\ldots x_n, y_1\ldots y_n)$ in $2n$ variables is invariant under the following operations:
1) $p(x_1\ldots x_n, y_1\ldots y_n)=p(y_1\ldots y_n, x_1\ldots x_n)$
2) $\forall \sigma\in S_n, p(x_1\ldots x_n, y_1\ldots y_n)=p(y_{\sigma(1)}\ldots y_{\sigma(n)}, x_{\sigma(1)}\ldots x_{\sigma(n)})$
In other words, the polynomial has a symmetry group $S_n \times \mathbb{Z}_2$.
My question is: is there a simple polynomial basis for polynomials of degree $\leq d$ with this symmetry?
Clearly one can find a linear basis for such polynomials by taking monomials and applying all elements of the symmetry group. For example, a linear basis for this space is given by polynomials of the form
$\displaystyle\sum_{i_1\neq i_2 \neq \ldots i_d =1}^{n} \left(x_{i_1}^{\alpha_1} x_{i_2}^{\alpha_2} \ldots x_{i_d}^{\alpha_d} y_{i_1}^{\beta_1} y_{i_2}^{\beta_2} \ldots y_{i_d}^{\beta_d}+ y_{i_1}^{\alpha_1} y_{i_2}^{\alpha_2} \ldots y_{i_d}^{\alpha_d} x_{i_1}^{\beta_1} x_{i_2}^{\beta_2} \ldots x_{i_d}^{\beta_d} \right)$
where $\alpha_1,\ldots, \alpha_d, \beta_1, \ldots \beta_d$ are a string of integers summing to $\leq d$. The size of this basis scales exponentially with $d$.
I'm asking if there's a much simpler polynomial basis. (In other words, if every such polynomial $p$ can be written as a polynomial $q$ in some simple basis elements, where $q$ is a generic polynomial). Ideally the number of elements in the basis would grow only polynomially with $d$. For example, in the case of polynomials on $n$ variables with symmetry group $S_n$, the elementary symmetric polynomials of degree $\leq d$ are a simple polynomial basis for the space of $S_n$-symmetric polynomials of degree $\leq d$, with merely $d$ elements in the basis. (In contrast, a linear basis for this space has many elements, namely the number of partitions of $d$.) I'm asking if there is an analogous polynomial basis known for the case of $2n$ variables and symmetry group $S_n\times \mathbb{Z}_2$. Any suggestions are much appreciated.
Denote $G=S_n \times \mathbb{Z}_2$. Let the group $G$ acts on the polynomial ring $k[X_{i,j}], i \leq n, l =\{0,1\}$ ($k$ be a field of chracteristic $0$) by $(\sigma, l)=X_{\sigma(i),{j+l \mod 2}} $. To find a basis of the algebra $k[X_{i,j}]^G$ of $G$-invariants you should use the $G$-homomorphism $k[X_{i,j}] \to k[X_{i,j}]^G$, $f \mapsto R(f)$ where $R$ is the Reinolds average operator $$ R=\frac{1}{|G|}\sum_{g \in G}g. $$ Then algebra of invariants $k[X_{i,j}]^G$ is generated by the elements $R(f)$ where $f$ runs all polynomial of $k[X_{i,j}]$ up to degree $2 n!$. But of course, it is not a minimal generating set. I hope the upper bound for the degree of invariants is $n.$
Some calculation for $n=4.$
Degree 1. There is only one invariant ( in terms of $x,y$): $$ y_{{2}}+y_{{3}}+y_{{4}}+y_{{1}}+x_{{2}}+x_{{3}}+x_{{4}}+x_{{1}}. $$ Degree 2. There is $ 3$ linearly independed invariants $$ {x_{{1}}}^{2}+{x_{{2}}}^{2}+{x_{{3}}}^{2}+{x_{{4}}}^{2}+{y_{{1}}}^{2}+ {y_{{2}}}^{2}+{y_{{3}}}^{2}+{y_{{4}}}^{2},\\ x_{{1}}y_{{1}}+x_{{2}}y_{{2}}+x_{{3}}y_{{3}}+x_{{4}}y_{{4}},\\ x_{{1}}y_{{2}}+x_{{1}}y_{{3}}+x_{{1}}y_{{4}}+x_{{2}}y_{{1}}+x_{{2}}y_{ {3}}+x_{{2}}y_{{4}}+x_{{3}}y_{{1}}+x_{{3}}y_{{2}}+x_{{3}}y_{{4}}+y_{{1 }}x_{{4}}+x_{{4}}y_{{2}}+x_{{4}}y_{{3}} $$ Degree 3. I have found $6$ invariants $$ {x_{{1}}}^{3}+{x_{{2}}}^{3}+{x_{{3}}}^{3}+{x_{{4}}}^{3}+{y_{{1}}}^{3}+ {y_{{2}}}^{3}+{y_{{3}}}^{3}+{y_{{4}}}^{3} ,\\ {x_{{1}}}^{2}y_{{1}}+{x_{{2}}}^{2}y_{{2}}+{x_{{3}}}^{2}y_{{3}}+{x_{{4} }}^{2}y_{{4}}+{y_{{1}}}^{2}x_{{1}}+{y_{{2}}}^{2}x_{{2}}+{y_{{3}}}^{2}x _{{3}}+{y_{{4}}}^{2}x_{{4}},\\ \left( x_{{1}}x_{{3}}+x_{{1}}x_{{4}}+x_{{1}}x_{{2}} \right) y_{{1}}+ \left( x_{{2}}x_{{3}}+x_{{1}}x_{{2}}+x_{{2}}x_{{4}} \right) y_{{2}}+ \left( x_{{3}}x_{{4}}+x_{{1}}x_{{3}}+x_{{2}}x_{{3}} \right) y_{{3}}+ \left( x_{{1}}x_{{4}}+x_{{2}}x_{{4}}+x_{{3}}x_{{4}} \right) y_{{4}}+ \left( x_{{1}}+x_{{2}} \right) y_{{2}}y_{{1}}+ \left( x_{{1}}+x_{{3}} \right) y_{{3}}y_{{1}}+ \left( x_{{4}}+x_{{1}} \right) y_{{4}}y_{{1}} + \left( x_{{2}}+x_{{3}} \right) y_{{3}}y_{{2}}+ \left( x_{{2}}+x_{{4} } \right) y_{{4}}y_{{2}}+ \left( x_{{3}}+x_{{4}} \right) y_{{4}}y_{{3} },\\ \left( {x_{{4}}}^{2}+{x_{{3}}}^{2}+{x_{{2}}}^{2} \right) y_{{1}}+ \left( {x_{{1}}}^{2}+{x_{{4}}}^{2}+{x_{{3}}}^{2} \right) y_{{2}}+ \left( {x_{{2}}}^{2}+{x_{{1}}}^{2}+{x_{{4}}}^{2} \right) y_{{3}}+ \left( {x_{{3}}}^{2}+{x_{{1}}}^{2}+{x_{{2}}}^{2} \right) y_{{4}}+ \left( x_{{4}}+x_{{2}}+x_{{3}} \right) {y_{{1}}}^{2}+ \left( x_{{1}}+ x_{{3}}+x_{{4}} \right) {y_{{2}}}^{2}+ \left( x_{{1}}+x_{{2}}+x_{{4}} \right) {y_{{3}}}^{2}+ \left( x_{{3}}+x_{{1}}+x_{{2}} \right) {y_{{4} }}^{2},\\ x_{{1}}x_{{2}}x_{{3}}+x_{{1}}x_{{2}}x_{{4}}+x_{{1}}x_{{3}}x_{{4}}+x_{{ 2}}x_{{3}}x_{{4}}+y_{{1}}y_{{2}}y_{{3}}+y_{{1}}y_{{2}}y_{{4}}+y_{{1}}y _{{3}}y_{{4}}+y_{{2}}y_{{3}}y_{{4}},\\ \left( x_{{2}}x_{{4}}+x_{{2}}x_{{3}}+x_{{3}}x_{{4}} \right) y_{{1}}+ \left( x_{{3}}x_{{4}}+x_{{1}}x_{{4}}+x_{{1}}x_{{3}} \right) y_{{2}}+ \left( x_{{1}}x_{{4}}+x_{{2}}x_{{4}}+x_{{1}}x_{{2}} \right) y_{{3}}+ \left( x_{{1}}x_{{3}}+x_{{1}}x_{{2}}+x_{{2}}x_{{3}} \right) y_{{4}}+ \left( x_{{3}}+x_{{4}} \right) y_{{2}}y_{{1}}+ \left( x_{{2}}+x_{{4}} \right) y_{{3}}y_{{1}}+ \left( x_{{2}}+x_{{3}} \right) y_{{4}}y_{{1}} + \left( x_{{4}}+x_{{1}} \right) y_{{3}}y_{{2}}+ \left( x_{{1}}+x_{{3} } \right) y_{{4}}y_{{2}}+ \left( x_{{1}}+x_{{2}} \right) y_{{4}}y_{{3} }, $$
and so on..