Un boceto de una prueba de que $c(n)=p_{DO}(n)$ se puede encontrar aquí; se trata de una sencilla manipulación de diagramas de Ferrers. El resto parece más fácil que ver con la generación de funciones.
Primera nota de que en
$$\prod_{k\ge 1}\frac1{1+x^k}=\prod_{k\ge 1}(1-x^k+x^{2k}-+\dots)\;,$$
el individuo $x^n=x^{k_1}x^{k_2}\dots x^{k_m}$ términos son positivos o negativos de acuerdo como $m$ es par o impar, y por lo tanto
$$\prod_{k\ge 1}\frac1{1+x^k}=\sum_{k\ge 0}\big(p_E(k)-p_O(k)\big)x^k\,:$$
una partición de $(k_1,k_2,\dots,k_m)$ se cuenta de forma positiva al $m$ es incluso negativa y al $m$ es impar.
Del mismo modo, en
$$\prod_{k\ge 0}(1-x^{2k+1})$$
el $x^n$ términos son todos de la forma $x^{2k_1+1}x^{2k_2+1}\dots x^{2k_m+1}$, donde el $k_i$ son distintos, y el término es positivo incluso para $m$ y negativo para los impares $m$, por lo que
$$\prod_{k\ge 0}(1-x^{2k+1})=\sum_{k\ge 0}(-1)^kp_{DO}(k)x^k\;.$$
Si podemos demostrar que
$$\prod_{k\ge 0}(1-x^{2k+1})=\prod_{k\ge 1}\frac1{1+x^k}\;,\tag{1}$$
vamos a ser capaces de concluir que el $p_E(k)-p_O(k)=(-1)^kp_{DO}(k)$ y que, por ende,$p_{DO}(k)=(-1)^k(p_E(k)-p_O(k))$.
Para demostrar $(1)$, se observa que la
$$\begin{align*}
\prod_{k\ge 1}\frac1{1+x^k}&=\prod_{k\ge 1}\frac{1-x^k}{1-x^{2k}}\\\\
&=\frac{\prod\limits_{k\ge 1}(1-x^k)}{\prod\limits_{k\ge 1}(1-x^{2k})}\\\\
&=\frac{\prod\limits_{k\ge 1}(1-x^{2k})\prod\limits_{k\ge 0}(1-x^{2k+1})}{\prod\limits_{k\ge 1}(1-x^{2k})}\\\\
&=\prod_{k\ge 0}(1-x^{2k+1})\;.
\end{align*}$$