[计算几何][找规律]迎风一刀斩 PTA L3-006

迎着一面矩形的大旗一刀斩下，如果你的刀够快的话，这笔直一刀可以切出两块多边形的残片。反过来说，如果有人拿着两块残片来吹牛，说这是自己迎风一刀斩落的，你能检查一下这是不是真的吗？

注意摆在你面前的两个多边形可不一定是端端正正摆好的，它们可能被平移、被旋转（逆时针90度、180度、或270度），或者被（镜像）翻面。

这里假设原始大旗的四边都与坐标轴是平行的。

输入格式：

输入第一行给出一个正整数N（≤20），随后给出N对多边形。每个多边形按下列格式给出：

kx1​y1​⋯xk​yk​

其中k（2<k≤10）是多边形顶点个数；(xi​,yi​)（0≤xi​,yi​≤10^8）是顶点坐标，按照顺时针或逆时针的顺序给出。

注意：题目保证没有多余顶点。即每个多边形的顶点都是不重复的，任意3个相邻顶点不共线。

输出格式：

对每一对多边形，输出YES或者NO。

输入样例：

8
3 0 0 1 0 1 1
3 0 0 1 1 0 1
3 0 0 1 0 1 1
3 0 0 1 1 0 2
4 0 4 1 4 1 0 0 0
4 4 0 4 1 0 1 0 0
3 0 0 1 1 0 1
4 2 3 1 4 1 7 2 7
5 10 10 10 12 12 12 14 11 14 10
3 28 35 29 35 29 37
3 7 9 8 11 8 9
5 87 26 92 26 92 23 90 22 87 22
5 0 0 2 0 1 1 1 2 0 2
4 0 0 1 1 2 1 2 0
4 0 0 0 1 1 1 2 0
4 0 0 0 1 1 1 2 0

输出样例：

YES
NO
YES
YES
YES
YES
NO
YES

题意： 给出两个多边形，询问是否能通过劈一刀大矩形得到这两个多边形。

分析： 切割一块大矩形能得到的情况一共只有五种，如下图所示： 

可以发现这两个小多边形的组合是一定是以下几种之一：四边形+四边形（再细分两种情况）、五边形+三角形、三角形+三角形、四边形+三角形。于是题目给出其它组合时就可以直接判否了，而对于这几种组合可以枚举出来能拼成矩形的所有情况，然后再找到共同的特征，如果两个n边型具有n-2个直角且那两个非直角夹边长度相同，那就可以用这两条长度相同的边作为公共边拼接成矩形。

有两种特殊情况要特判一下，第一种就是把大矩形且为两个小矩形时，这种情况下所有角都是直角，没法用上面的判断方法，所以需要枚举一下各边长，一旦出现了等长的两边就说明两矩形可以拼接。第二种是下图的情况：

这种情况虽然满足上面的判断方法，但却不能够拼成一个大矩形，所以需要求出两梯型的直角腰长度，若等长才可以拼成大矩形。

最后一个坑让我debug了好长时间，后来发现是题面描述错误，题面的意思是给你两个多边形，它可能被平移、旋转、翻面，然后检查其是否可以拼为矩形，但出题人的意思是给你一个矩形，切一刀产生两个多边形，并可能对两多边形平移、旋转、翻面，问这两个多边形是否来自该矩形。

这两种表述差别在哪？例如下图两多边形，如果按照题面表述明显返回YES，但按照出题人的意思就是NO，因为初始矩形各边是与坐标轴平行的，而旋转只能转一个直角，所以最后应该至多出现1条与坐标轴不平行的边。

 

具体代码如下：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <string>
#include <vector>
#define x first
#define y second
#define int long long 
#define pii pair<int, int>
using namespace std;

int px(pii a, pii b){//点乘 
	return a.x*b.x+a.y*b.y;	
}
pii operator-(pii a, pii b){
	return make_pair(a.x-b.x, a.y-b.y);
}

signed main()
{
	int n;
	cin >> n;
	for(int i = 1; i <= n; i++){
		int n1, n2;
		vector<pii> v1, v2;
		cin >> n1;
		for(int j = 1; j <= n1; j++){
			pii p;
			cin >> p.x >> p.y;
			v1.push_back(p);
		}
		cin >> n2;
		for(int j = 1; j <= n2; j++){
			pii p;
			cin >> p.x >> p.y;
			v2.push_back(p);
		}
		if(!(n1==4&&n2==4||n1==3&&n2==3||max(n1,n2)==5&&min(n1,n2)==3||max(n1,n2)==4&&min(n1,n2)==3)){
			puts("NO");
			continue;
		}
		int p1, p2, p3, p4, num1 = 0, num2 = 0;//直角个数 
		int cnt1 = 0, cnt2 = 0;
		p1 = p2 = p3 = p4 = -1;
		int _p1, _p2, _p3, _p4;//找到直角边 
		_p1 = _p2 = _p3 = _p4 = -1;
		for(int j = 0; j < n1; j++){
			pii a = v1[j], b = v1[(j+1)%n1], c = v1[(j+2)%n1];
			int t = px(b-a, c-b);
			if(t != 0){
				if(p1 == -1) p1 = j;
				else p2 = j;
			}
			else{
				num1++;
				if(_p1 == -1) _p1 = j;
				else _p2 = j;
			}
			if(v1[j].x != v1[(j+1)%n1].x && v1[j].y != v1[(j+1)%n1].y) cnt1++;
		}
		for(int j = 0; j < n2; j++){
			pii a = v2[j], b = v2[(j+1)%n2], c = v2[(j+2)%n2];
			int t = px(b-a, c-b);
			if(t != 0){
				if(p3 == -1) p3 = j;
				else p4 = j;
			}
			else{
				num2++;
				if(_p3 == -1) _p3 = j;
				else _p4 = j;
			}
			if(v2[j].x != v2[(j+1)%n2].x && v2[j].y != v2[(j+1)%n2].y) cnt2++;
		}
		if(cnt1 >= 2 || cnt2 >= 2){
			puts("NO");
			continue;
		}
		if(num1 == 4 && num2 == 4){//特判两矩形 
			int t1 = (v1[0].x-v1[1].x)*(v1[0].x-v1[1].x)+(v1[0].y-v1[1].y)*(v1[0].y-v1[1].y);
			int t2 = (v1[1].x-v1[2].x)*(v1[1].x-v1[2].x)+(v1[1].y-v1[2].y)*(v1[1].y-v1[2].y);
			int t3 = (v2[0].x-v2[1].x)*(v2[0].x-v2[1].x)+(v2[0].y-v2[1].y)*(v2[0].y-v2[1].y);
			int t4 = (v2[1].x-v2[2].x)*(v2[1].x-v2[2].x)+(v2[1].y-v2[2].y)*(v2[1].y-v2[2].y);
			if(t1 == t3 || t2 == t3 || t1 == t4 || t2 == t4)
				puts("YES");
			else
				puts("NO");
			continue;
		}
		if(num1 != n1-2 || num2 != n2-2){
			puts("NO");
			continue;
		}
		pii U1 = v1[(p1+1)%n1], V1 = v1[(p2+1)%n1];
		pii U2 = v2[(p3+1)%n2], V2 = v2[(p4+1)%n2];
		pii _U1 = v1[(_p1+1)%n1], _V1 = v1[(_p2+1)%n1];
		pii _U2 = v2[(_p3+1)%n2], _V2 = v2[(_p4+1)%n2];
		//判两斜边相等 
		if((U1.x-V1.x)*(U1.x-V1.x)+(U1.y-V1.y)*(U1.y-V1.y) == (U2.x-V2.x)*(U2.x-V2.x)+(U2.y-V2.y)*(U2.y-V2.y)){
			if(n1 == 4 && n2 == 4){//分割为两四边形时还需要特判 
				//判两直角边相等 
				if((_U1.x-_V1.x)*(_U1.x-_V1.x)+(_U1.y-_V1.y)*(_U1.y-_V1.y) == (_U2.x-_V2.x)*(_U2.x-_V2.x)+(_U2.y-_V2.y)*(_U2.y-_V2.y))
					puts("YES");
				else
					puts("NO");
			}
			else
				puts("YES");
		}
		else
			puts("NO");
	}
    return 0;
}