Given Ramanujan tau function $\tau(n)$, which is the nth Fourier coefficient of the modular discriminant
$\displaystyle \Delta(q)=q\prod_{m=1}^\infty (1-q^m)^{24} = \sum_{n=1}^\infty \tau(n)\,q^n\tag{1a}$
Nothing forbids us from generalizing the generating function to
$\displaystyle \Delta_{k}(q)= \sum_{n=1}^\infty \tau((2k-1)n)\,q^n\tag{1b}$
for integers $k\gt0$
After some experimentation with wolfram mathematica, one discovers the following intriguing patterns
$\Delta_{1}(q)=q^{\color{blue}{1}}-24q^2+252q^3-1472q^4+4830q^5-6048q^6-16744q^7+84480q^8-113643q^{\color{blue}{9}}-115920q^{10}+543612q^{11}-\cdots$
$\Delta_2(q)=252q-6048q^2-113643q^{\color{brown}{3}}-370944q^4+1217160q^5+2727432q^6-4219488q^7+21288960q^8-73279080q^9-29211840q^{10}+134722224q^{11}+167282496q^{12}-145589976q^{13}+101267712q^{14}-548895690q^{15}+248758272q^{16}-1740295368q^{17}+1758697920q^{18}+2686677840q^{19}-1791659520q^{20}+1902838392q^{21}-3233333376q^{22}+4698104544q^{23}-9600560640q^{24}-6425804700q^{25}+3494159424q^26+1665188361q^{\color{brown}{27}}+\cdots$
$\vdots$
From which we deduce the following conjecture
The nth coefficient $\tau((2k-1)n)$ is odd if and only if the powers of q in the partial sum of $(1b)$ depend on whether $(2k-1)$ is square-free or non square-free as follows
$n\stackrel{\mathrm{def}}{=} \begin{cases}{\color{brown}{(2k-1),9(2k-1),25(2k-1),49(2k-1),...}} & \mbox{ where $2k-1=1,3,5,7,11,13,...$ is any odd square-free integer}\\{\color{blue}{1,9,25,49,81,...}} & \mbox{ where $2k-1=9,25,27,49,81,...$ is any odd non square-free integer} \end{cases}$
Question: How do we prove the conjectured criterion?
Motivation: The Möbius function $\mu(n)$ is either $1$ or $-1$ if $n$ is square-free depending on whether it has a number of even prime factors or odd respectively and $\mu(n)=0$ if it is non square-free. It happens to be a very important arithmetic function in number theory. For example, the Dirichlet series that generates the Möbius function is the reciprocal of the Riemann zeta function $\displaystyle \sum_{n=1}^{\infty} \frac{\mu(n)}{n^s}=\frac{1}{\zeta(s)}$