Referring to the formula for how much bulging occurs for a given speed (up to on the first order but it's good enough), we have that the difference in polar and equatorial radii relative to the mean radius $a = \frac{a_e + a_e + a_p}{3}$ is that $$\frac{a_e - a_p}{a} = \frac{5}{4} \frac{\omega^2 a^3}{GM}$$.
The only other constraint on $a_e, a_p$ we have is the volume of the planet derived from the mass and density, $\rho \cdot \frac{4 \pi}{3} a_e^2 a_p = M$. And what we want from the equatorial acceleration and polar acceleration is that $$3 g_e = 3 \left( \frac{GM}{a_e^2} - a_e \omega^2 \right) = q_p = \frac{GM}{a^2_p}$$
ignoring the non-spherical components of gravity. We can't decide the exact values as
We know from the volume constraint that $a_p = \frac{a_V^3}{a_e^2}$, where $a_V$ is the radius of the spherical planet. This turns the other two into an equation on $a_e$ and $\omega$. For the gravity magnitude relation, we get a quadratic on $a_e^3$, so $$a_e = \sqrt[3]{\frac{-3\omega^2 + \sqrt{9\omega^4 + 12 \frac{G^2M^2}{a^6_V}}}{2 \frac{GM}{a^6_V}}} \text{ or } \omega^2 = \frac{GM}{a_e^3} \left( 1 - \frac{a_e^6}{3a_V^6}\right)$$
Let $a_e = f a_V$ and we get the final monster formula when subbing into the bulge relation, $$\frac{3f^3 - 3}{2 f^3 + 1} = \left(1 - \frac{f^3}{3} \right) \frac{5}{4} \frac{\left(\frac{2}{3}f^3 + \frac{1}{3}\right)^3}{f^9}$$
This is an absolute monster of a polynomial, a quintic of $f^3$ and definitely not solvable by analytical methods. By graph, the answer is roughly $1.16017$. We could get more decimals, but we have enough error from our approximations.

We thus get our variables as
$$ \begin{align}
a_v &= \sqrt[3]{\frac{3M}{4 \pi \rho} } = 7458 \text{ km} \\
a_e &= f \cdot a_V = 8652 \text{ km} \\
a_p &= f^{-2} \cdot a_V = 5541 \text{ km} \\
\omega &= \sqrt{\frac{4 \pi \rho G}{3 f^3} \left( 1 - \frac{1}{3} f^6\right) } = 4.293 \cdot 10^{-4} \text{Hz} \\
\implies T &= \frac{2 \pi }{\omega} = 4 \text{ hours } 3 \text{ minutes } 56 \text{ seconds}
\end{align}$$
Yeah that's fast. Definitely beyond the range of the approximations taken here.
Edit : BTW the gravitational acceleration at equator is 0.707 G and at the poles is 2.12G. If you scale your radius by $s$ while keeping density fixed, mass changes by $s^3$, while acceleration changes by $s$, but the time period remains fixed. So to get 0.4G at equator and 1.2 at poles at the same density, the planet must have a mass of $0.3 M_\oplus$, radii of $a_V = 4262 \text{ km}, a_e = 3166 \text{ km}, a_p = 4944 \text{ km}$.