I am now playing with permutation matrices, http://mathworld.wolfram.com/PermutationMatrix.html.
Also, there is a similar discussion: Symmetric Permutation Matrix. I want to ask more details than this one.
As we know, a permutation matrix is orthogonal, i.e., $E^T=E^{-1}$. I am interested in when it is symmetric, i.e., $E^T=E^{-1} = E$
Suppose
- Start from an identity matrix $I_n$.
- $n$ can be even or odd number.
- Pick $(i,j)$, where $0<i,j\leq n$ and $i, j$ are integer. Exchange $i$-th and $j$-th columns of $I_n$ (identity matrix) and get $E$. Then $E$ is symmetric. This is because $E_{ii}=E_{jj}=0$ and $E_{ij}=E_{ji}=1$.
- Based on 3., if I pick a number of sets $(i,j)$, $(k,l)$, $(q,r), \ldots$, without repeated index in each $(\cdot,\cdot)$, and permute columns of $I_n$ according to these sets, then the resulting permutation matrix $E$ is symmetric.
One key thing here is "without repeated index in each $(\cdot,\cdot)$". This is because if I do $(1,2)$ and $(2,3)$ for $I_3$ for example, I get
$$\begin{bmatrix}0 & 0 & 1 \\ 1 & 0 & 0 \\ 0 & 1 & 0 \end{bmatrix},$$
which is not symmetric. In this case, I repeat $2$ in each suit.
Is the above correct? Or I miss some key assumptions?