由于要求10位小数的精度,可以想像模拟法会需要$10^{10}$以上数量级的试验,解析计算则非常复杂,一种手段是离散化,即先假定双方在$(0,1/n,2/n,\ldots,(n-1)/n,1)$之间选数,再令$n\to\infty$。
此处直接引用此文的结论得x的概率密度为$e-e^{1-x}$, y的概率密度为$e^2-e-e^{2-y}+(1-y)e^{1-y}$,从而y大于x的概率为
$$\int_{0}^{1}\left(\int_0^y(e-e^{1-x})dx\right)(e^2-e-e^{2-y}+(1-y)e^{1-y})dy=\int_0^1(ey+e^{1-y}-e)(e^2-e-e^{2-y}+(1-y)e^{1-y})dy=\frac{1}{4}(1+14e-5e^2)\approx0.5276662759.$$
本题无需编程。
|