You have to solve two equations simultaneously:
$$\begin{align*}
I_{sat_1}\left(e^{\frac{V_{D_1}}{n V_T}}-1\right) &= I_{sat_2}\left(e^{\frac{V_{D_2}}{n V_T}}-1\right)\\\\
V_0&=V_{D_1} +V_{D_2}= -3\:\textrm{V}
\end{align*}$$
This follows from the fact that the currents in both diodes must be equal and that the sum of their voltages must match with the supply voltage. Pretty obvious, really.
Solving these simultaneously is a bit tricky. You could attempt it, iteratively. Or you could attempt it with the Lambert-W (aka ProductLog) function. (Doable for a closed solution, but still takes some work.)
But we can use a symmetry argument to claim that -1 term can be ignored. This allows for a very simple solution:
$$\begin{align*}
V_{D_1} &\approx -\frac{1}{2} \left[V_T\: \operatorname{ln}\left(\frac{I_{sat_2}}{I_{sat_1}}\right) + V_0\right]\\\\
V_{D_2} &= V_0 - V_{D_1}
\end{align*}$$
Using \$V_T=26\:\textrm{mV}\$ and \$I_{sat_1}=0.1\:\textrm{pA}\$ and \$I_{sat_2}=5\:\textrm{pA}\$, this instantly gives the correct answer to quite a few places: \$V_{D_1}\approx -1.5508563\:\textrm{V}\$. This should check out okay.
You can "read" the above solution equation to say:
Start out with the assumption that the voltage is divided in half.
Then apply a correction which will be one half of \$V_T\$ times the
logarithm of the ratio. (The sign of the correction will, of course,
depend upon which saturation currents are used in the numerator and
denominator.)
It's my suspicion that this is the approach you were supposed to take because it focuses on what's important (the saturation current ratios) and avoids getting hung up on numerical solutions or overly mathy discussions which distract from rather than shed light on the subject.
A note about a huge assumption in all this is about the nominal temperature. It's enough to simply say "room temperature" and use a value for \$V_T\$ that is commonly used (somewhere from about \$25\:\textrm{mV}\$ to perhaps \$26\:\textrm{mV}\$ is often picked.) However, to imagine recomputing that for any temperature to get the right behavior from the above equations by just substituting in the new value for \$V_T\$ is wrong. It turns out that the saturation currents are a function of \$T^3\$ to \$T^4\$ and so they vary as well. In fact, they vary so much that they overwhelm the effect of \$V_T\$ enough to reverse the sign of the effect!
So it is probably fine to just assume that the saturation currents one is given are meant at "room temperature." But that's all. If the model is to apply over a wide temperature range, the variation of the saturation currents must also be incorporated into the model. And that is a whole other issue.