We have two conjectured generalizations of the question asked at a Hankel matrix of involution numbers by Tewodros Amdeberhan. Let $n!!=1!\,2!\cdots n!$.
Conjecture 1. Let $I_k$ denote the number of involutions in the symmetric group $\mathfrak{S}_k$. Then the Smith normal form of the matrix $[I_{i+j}]_{i,j=0}^n$ has diagonal entries $0!, 1!, \dots, n!$,
Conjecture 2. Let $J_k=k!\sum_{\lambda\vdash k}s_\lambda$, where $s_\lambda$ is a Schur function. Set $J_0=1$. When the symmetric function $\det[J_{i+j}]_{i,j=0}^n$ is expanded in terms of power sums, then every coefficient is an integer divisible by $n!!$. (Tewodros' question is equivalent to the coefficient of $p_1^{n(n+1)}$ being equal to $n!!$.)
Is there a nice formula or combinatorial interpretation of the coefficients in Conjecture 2?