For $n \geq 0$, let $a_n$ be the square of the Euclidean length of the vector of Littlewood-Richardson coefficients of $\sum_{\lambda \vdash n} s_\lambda^2$, where $s_\lambda$ are the symmetric Schur functions and the sum runs over all partitions $\lambda$ of $n$. These numbers can also be described via the generating function $\sum_{n \geq 0} a_n x^n = \prod_{i \geq 1} (1 - 4x^i)^{-1/2}$. This sequence appears in the OEIS (https://oeis.org/A067855) and begins 1, 2, 8, 26, 94, 326, 1196, 4358, $\dots$.
For $n > 1$, it appears empirically that $a_n$ is congruent to 0 (resp. 2) mod 4 when the exponent in the highest power of 2 dividing $n-1$ is odd (resp. even). Is this true?
(Additionally, it appears that $a_{4k}-a_{4k-2}$ is always 0 mod 8, while $a_{2^{i+1} k}-a_{2^{i+1} k- 2^i}$ is always 4 mod 8 when $i \geq 2$. There are similar but more complicated patterns modulo higher powers of 2. For instance, $a_{16k}-a_{16k-2}$ is 0 (resp. 32) mod 64 according to whether $k$ is 0, 4, or 5 (resp. 1, 2, or 3) mod 6.)