手动计算 ACM

题目描述
LF是一个计算几何大佬，他每天都要做114514道几何题来锻炼自己的计算几何能力。但是今天LF做剩一道题的时候他竟然发现他这道题他不会做，这道题的题意是这样的：
平面上有两个椭圆，这两个椭圆的中心都是(0,0)，且一个椭圆的两个焦点在x轴上，一个椭圆的两个焦点在y轴上，给出这两个椭圆的方程，求他们面积的并的大小。
这道题实在是太难了，于是LF找到了你，你能帮他完成这个艰巨的任务吗？
输入描述:
输入第一行一个整数T\ (1 \leq T \leq 5)T (1≤T≤5)，表示数据组数。
接下来共T行，每行四个实数a,b,c,d，表示第一个椭圆的方程为\frac{x^2}{a2} +\frac{y^2}{b2}=1
a
2

x
2

​
+
b
2

y
2

​
=1,第二个椭圆的方程为\frac{x^2}{c2} +\frac{y^2}{d2}=1
c
2

x
2

​
+
d
2

y
2

​
=1, 数据保证a,b,c,d最多两位小数, 且1 \leq b<a\leq 8,1 \leq c<d\leq 81≤b<a≤8,1≤c<d≤8。
输出描述:
输出共T行，每行一个实数表示该组数据中这两个椭圆面积的并的大小，保留一位小数。
示例1
输入
复制
3
2 1 1 2
3 2 1 2
3 2 1 4
输出
复制
8.9
18.8
23.9

做法：定积分
求除两个椭圆焦点的横坐标(方程联立），计算第一象限内两个椭圆交集的面积（定积分（a，b，c，d，x0都被视为常数）），交集的面积等于两个椭圆分别在0到x0和x0到c的积分相加
总共的交集面积是第一象限的4倍
求出两个椭圆的面积，减去交集的面积

特殊情况，一个椭圆包含在另一个椭圆，独立判断

#include <iostream>
#include <cstdio>
#include <cmath>

using namespace std;

const double PI = 3.1415926535;
/*
3
2 1 1 2
3 2 1 2
3 2 1 4
*/
int main(void)  {
	int t;
	cin >> t;
	while (t--) {
		double a, b, c, d, ans = 0;
		cin >> a >> b >> c >> d;
		if (a > c && b > d) {
			ans = PI * a * b;
		} else if (c > a && d > b) {
			ans = PI * c * d;
		} else {
			double x0 = sqrt((a*a * c*c * (b*b - d*d))/(b*b * c*c - a*a * d*d));
			double s1 = 0.5 * a * b * (x0 / a) * sqrt(1 - x0*x0 / (a*a)) + 0.5 * a * b * (asin(x0 / a));
			double s2 = PI * c * d / 4 -  0.5 * c * d *  x0 / c * sqrt(1 - x0 * x0 / (c * c)) - 0.5 * c * d * asin(x0 / c);
			ans = PI * a * b + PI * c * d - 4 * (s1 + s2);
		}
		printf("%.1f\n", ans);
	}
}