Let $x(t)$ be a bandlimited signal such that $X(j\omega) =0 $ when $|\omega|>M$. Also $p(t) = p_1(t) - p_1(t-\Delta)$ is a nonuniformly spaced periodic pulse train where $$p_1(t) = \sum_{k = -\infty}^{+\infty}\delta\left(t - \frac{2\pi k}{M}\right), \quad \text{with}\quad \Delta= \frac{\pi}{2M}$$
Let $x_p(t) = x(t)p(t)$ and apply an ideal low-pass filter with cutoff frequency $\omega_c = M$ to $X_p(j\omega)$. The result is $z(t)$. Design a system which recovers $x(t)$ from $z(t)$.
My try:
It's easy to see that $$P(j\omega) = \left(1 - e^{-j\Delta \omega}\right)\left(M\sum_{k = -\infty}^{+\infty}\delta(\omega - kM)\right) = M\sum_{k = -\infty}^{+\infty}\left(1 - e^{-j\Delta kM}\right)\delta(\omega - kM) $$
Also we have
$$ X_p(j\omega) = \frac{1}{2\pi}X(j\omega)\star P(j\omega) = \frac{M}{2\pi}\sum_{k = -\infty}^{+\infty}\left(1 - e^{-j\Delta kM}\right)X\big(j(\omega - kM)\big) $$
Since $\Delta= \frac{\pi}{2M}$ we have $$ X_p(j\omega) = \frac{M}{2\pi}\sum_{k = -\infty}^{+\infty}\left(1 - e^{\frac{-jk\pi}{2}}\right)X\big(j(\omega - kM)\big) $$
After low-pass filter we get
$$ Z(j\omega) = \begin{cases} \frac{M}{2\pi}\bigg[(1+j)X\big(j(\omega - M)\big) + (1-j)X\big(j(\omega + M)\big)\bigg], & \lvert \omega\rvert < M\\ 0, & \text{O.W} \end{cases} $$
I've got stuck here. How can we recover $x(t)$? I tried to use phase shifter and $z(t)\cos(Mt)$ but it didn't work.