ZOJ 4111 Square on the Plane (数值分析+分类讨论)

Source The 16th Zhejiang Provincial Collegiate Programming Contest Sponsored by TuSimple - Problem L

组团狙击 大 失 败，还是 Claris 比较快，
不过讲道理其他题都不是我做的，拼不过罚时也不能怪我现场没过这个没什么希望现场过的 L 题吧（正论）。

总之先画个图，

所求即为点 $P$ 从 $A(a,0)$ 出发沿着椭圆 $\frac{x^2}{a^2}+\frac{y^2}{b^2}=1$ 逆时针走 $t$ 个单位长度之后，以 $P$ 为中心 $2|y_P|$ 为边长的边与坐标轴平行的正方形 $S$ 扫过的面积（这里无论正方形是否包含内部，结果都是一样的），

不难发现，上半平面和下半平面是完全独立的，记椭圆的周长为 $c$，那么只需要考虑 $0\le t\le \frac{c}{2}$ 的情况，

先考虑给定 $P$ 计算 $t$ 也就是从 $A$ 出发沿着椭圆逆时针走到 $P$ 的长度，令
$$\{\begin{array}{l}{x}_P=a\cos \theta ,\\ y_P=b\sin \theta ,\end{array}$$ 那么有
$$t={\int }_0^{\alpha }\sqrt{a^2?(a^2?b^2){\cos }^2\theta }d\theta ,$$ 其中 $\alpha =\arccos (\frac{x_P}{a})$，就可以转化为第二类不完全椭圆积分，使用 std::tr1::ellint_2 函数进行计算，顺便也可以算出椭圆的周长 $c$，
现在对于给定的 $\alpha$，已经可以计算出 $t(\alpha )$，且易知 $t^{\prime }(\alpha )=\sqrt{a^2?(a^2?b^2){\cos }^2\alpha }$，就可以通过牛顿迭代在给定 $t$ 时解出 $\alpha$，进而求出点 $P$ 的坐标，实测取初值 $\frac{\pi }{2}$ 迭代 $4$ 次精度就够了，

接下来考虑正方形 $S$ 扫过的面积，分四种情况讨论，

第一种情况，椭圆在点 $P$ 处的切线斜率 $<?1$，即 $x_P>\frac{a^2}{\sqrt{a^2+b^2}}$，
此时扫过的面积就是 $t$ 时刻正方形的面积 $4y_P^2$，

第二种情况，椭圆在点 $P$ 处的切线斜率 $\in [?1,0]$，即 $0\le x_P\le \frac{a^2}{\sqrt{a^2+b^2}}$，
此时扫过的面积可以分为一个正方形的面积 $4y_P^2$ 和一个曲边梯形的面积，而曲边梯形的面积可以继续分为一个矩形的面积 $\frac{2b^2}{\sqrt{a^2+b^2}}?(\sqrt{a^2+b^2}?(x_P+y_P))$ 和一个曲边三角形的面积，
不难求出正方形的右上角 $Q$ 的轨迹方程为
$$\frac{(x?y/2{)}^2}{a^2}+\frac{y^2}{4b^2}=1,$$ 那么曲边三角形的面积为
$$s={\int }_{\frac{2b^2}{\sqrt{a^2+b^2}}}^{2y_P}(x(y)?(x_P+y_P))dy={\int }_{\frac{2b^2}{\sqrt{a^2+b^2}}}^{2y_P}(\sqrt{1?\frac{y^2}{4b^2}}+\frac{y}{2}?(x_P+y_P))dy,$$ 这个积分比较基础，这里就不展开了，

第三种情况，椭圆在点 $P$ 处的切线斜率 $>1$，即 $x_P<?\frac{a^2}{\sqrt{a^2+b^2}}$，
此时扫过的面积就是一整块的面积，可以转化成第二种情况代入横坐标 $0$ 进行计算，

第四种情况，椭圆在点 $P$ 处的切线斜率 $\in (0,1]$，即 $?\frac{a^2}{\sqrt{a^2+b^2}}\le x_P<0$，
此时扫过的面积是第三种情况中一整块的面积减去一个曲边梯形的面积，而曲边梯形的面积可以转化成第二种情况代入横坐标 $?x_P$ 进行计算，

至此所有情况都讨论完了，由于数据组数较多，还需要精细地实现以避免超时。