2618: [Cqoi2006]凸多边形

2618: [Cqoi2006]凸多边形

Time Limit: 5 Sec  Memory Limit: 128 MB
Submit: 1011  Solved: 515
[ Submit][ Status][ Discuss]

Description

逆时针给出 n个凸多边形的顶点坐标，求它们交的面积。例如n=2时，两个凸多边形如下图：

则相交部分的面积为5.233。

Input

第一行有一个整数n，表示凸多边形的个数，以下依次描述各个多边形。第i个多边形的第一行包含一个整数m_i，表示多边形的边数，以下m_i行每行两个整数，逆时针给出各个顶点的坐标。

 

Output

    输出文件仅包含一个实数，表示相交部分的面积，保留三位小数。

 

Sample Input

2
6
-2 0
-1 -2
1 -2
2 0
1 2
-1 2
4
0 -3
1 -1
2 2
-1 0

Sample Output

5.233

HINT

100%的数据满足：2<=n<=10，3<=m_i<=50，每维坐标为[-1000,1000]内的整数

Source

[ Submit][ Status][ Discuss]
做一遍半平面交，然后求交出来的凸多边形的面积
为了防止忘记这个麻烦的算法，，，这里加一些注释
半平面表示方法：点+向量表示一条直线，一个半平面为这条直线左方的区域 算法流程： 1.对所有半平面求极角序atan2(v.y,v.x)，按极角序从小到大排序 2.先将第一个半平面加入双端队列，定义sec[i]为第i个半平面与第i+1个半平面的交点 3.从第二个到最后一个for一遍，考虑当前插入的半平面，若sec[tail-1]在其右边，说明第tail个半平面会被它所取代，tail-- 4.还需考虑队头，总之，sec[head]必须在当前的半平面左边(因为半平面可能绕过一圈使得队头的半平面多余) 5.插入半平面，考虑当前半平面与上一个是否平行，若是，舍去左边的一个 6.求sec[tail]
叉积与判断点在半平面左侧就不解释了 主要是GetIntersection操作，苟蒻搞了好久老是忘，，总算证了一遍
这里求半平面a,v1  与半平面b,v2的交点(用点+向量，表达方法随意啦) 设交点为c，令c = a+t*v1 重点在求出t 作出向量u = a - b Cross(u,v2) Cross(v2,v1)分别是那两个带红边的三角形的有向面积的两倍 图画成这样时，有|v2| = |v2'|  v2//v2' 所以说，两三角形的面积比等于其高比，又由v2//v2'，根据相似三角形，高比又可转换为t与|v1|之比
求出所有交点以后，多边形面积计算公式 随意取一个点o，Ans = ∑Cross(sec[i]-o,sec[i+1]-o)  (if (i == tail) i+1 = head) 最后输出的是Ans/2

#include<iostream>
#include<cstdio>
#include<queue>
#include<vector>
#include<bitset>
#include<algorithm>
#include<cstring>
#include<map>
#include<stack>
#include<set>
#include<cmath>
#include<ext/pb_ds/priority_queue.hpp>
using namespace std;

const int maxn = 1E3;
typedef double DB;
const DB eps = 1E-10;
const DB two = 2.00;

struct Point{
	DB x,y;
	Point(){}
	Point(DB x,DB y): x(x),y(y){}
	Point operator - (const Point &b) {return Point(x - b.x,y - b.y);}
	Point operator * (const DB &b) {return Point(x*b,y*b);} 
	Point operator + (const Point &b) {return Point(x + b.x,y + b.y);}
}poi[maxn],sec[maxn];
typedef Point Vector;

struct Line{
	Point p; Vector v; DB ang;
	Line(){}
	Line(Point p,Vector v): p(p),v(v){ang = atan2(v.y,v.x);}
	bool operator < (const Line &b) const {return ang < b.ang;}
}line[maxn],Q[maxn];

int n,m,head,tail,cnt;

DB Cross(Vector v1,Vector v2) {return v1.x*v2.y - v2.x*v1.y;}
bool Onleft(Line k,Point po) {return Cross(k.v,po - k.p) > 0;}

Point GetIntersection(Line l1,Line l2)
{
	Vector u = l1.p - l2.p;
	DB t = Cross(u,l2.v)/Cross(l2.v,l1.v);
	return l1.p + l1.v*t;
}

int main()
{
	#ifdef DMC
		freopen("DMC.txt","r",stdin);
	#endif
	
	cin >> n;
	while (n--) {
		scanf("%d",&m);
		for (int i = 1; i <= m; i++) {
			DB x,y; scanf("%lf%lf",&x,&y);
			poi[i] = Point(x,y);
		}
		poi[++m] = poi[1];
		for (int i = 2; i <= m; i++) 
			line[cnt++] = Line(poi[i-1],poi[i] - poi[i-1]);
	}
	sort(line,line + cnt);
	Q[head = tail = 0] = line[0];
	for (int i = 1; i < cnt; i++) {
		while (head < tail && !Onleft(line[i],sec[tail-1])) --tail;
		while (head < tail && !Onleft(line[i],sec[head])) ++head;
		Q[++tail] = line[i];
		if (fabs(Cross(Q[tail].v,Q[tail-1].v)) < eps) {
			--tail;
			if (Onleft(Q[tail],line[i].p)) Q[tail] = line[i];
		}
		if (head < tail) sec[tail-1] = GetIntersection(Q[tail],Q[tail-1]);
	}
	while (head < tail && !Onleft(Q[head],sec[tail-1])) --tail;
	sec[tail] = GetIntersection(Q[head],Q[tail]); sec[tail+1] = sec[head];
	DB Ans = 0; Point o = Point(233.233,322.322);
	for (int i = head; i <= tail; i++) Ans += Cross(sec[i],sec[i+1]);
	printf("%.3lf",Ans/two);
	return 0;
}