先说几句批话:
某个菜鸡在干计算几何的时候被精度整出怄火症,两道题调了三天。
今天本来是要学半平面交的,但是想想还是写篇小说总结一下吧。
希望大家不要跟我一样犯下低级错误了!
由于我太菜了,很可能会出锅……
旋转卡壳入门
前置知识
向量
啊是吧。
凸多边形
凸多边形的切线
如果一条直线与凸多边形有交点,并且整个凸多边形都在这条直线的一侧,那么这条直线就是该凸多边形的一条切线。
对踵点
如果过凸多边形上两点作一对平行线,使得整个多边形都在这两条线之间,那么这两个点被称为一对对踵点。
凸多边形的直径
即凸多边形上任意两个点之间距离的最大值。直径一定会在对踵点中产生,如果两个点不是对踵点,那么两个点中一定可以让一个点向另一个点的对踵点方向移动使得距离更大。并且点与点之间的距离可以体现为线与线之间的距离,在非对踵点之间构造平行线,一定没有在对踵点构造平行线优。
感性理解一下吧。
旋转卡壳
这里直接放上一幅图:
对,字面意思。
以求凸包直径为例,在凸包上维护双指针, t t t 指针到 A i A i + 1 A_iA_{i+1} AiAi+1 这条边的距离最大,就是对踵点。此时直接用 ∣ A i A t ∣ |A_iA_t| ∣AiAt∣ 和 ∣ A i + 1 A t ∣ |A_{i+1}A_t| ∣Ai+1At∣ 更新答案就行了。
需要注意的是,旋转卡壳中, t t t 指针可能会转一圈再回来,毛估估应该最多转 ⌈ 3 t o p 2 ⌉ \lceil \frac{3top}{2} \rceil ⌈23top⌉ 次吧。
代码
这里用 s t st st 储存凸包上的点,其中 s t 1 = s t t o p = 1 st_1=st_{top}=1 st1=sttop=1,所以实际上我们只要枚举 i ∈ [ 1 , t o p ) i \in [1,top) i∈[1,top)。
struct node{
long long x,y;
bool operator<(const node&cwx)const{return x!=cwx.x?x<cwx.x:y<cwx.y;}//求凸包用的排序
node operator-(const node&cwx)const{return {x-cwx.x,y-cwx.y};}
long long operator*(const node&cwx)const{return x*cwx.y-y*cwx.x;}//叉积
long long dis2(){return x*x+y*y;}//模长的平方
}p[50001];
long long dis(node a,node b,node c){
long long tmp=(a-c)*(b-c);
return abs(tmp*tmp/(a-b).dis2());
}//c点到ab边的距离的平方
for(int i=1;i<top;++i){
while(dis(p[st[i]],p[st[i+1]],p[st[t]])<dis(p[st[i]],p[st[i+1]],p[st[t+1]]))t=t%(top-1)+1;//绝对不要写<=
ans=max(ans,max((p[st[i]]-p[st[t]]).dis2(),(p[st[i+1]]-p[st[t]]).dis2()));//dis表示t到当前边的距离
}
还有注意: t t t 的初值一定要设为 2 2 2,并且代码里面的 < < < 不能改成 ≤ \leq ≤。
为什么这样做不会遗漏所有有效情况?
首先,由于是凸包,所以铁定不会出现三点共线的情况。
最容易寄的也就这种情况。
按刚刚的做法看上去好像有点荒谬:那个紫色的点到当前边的距离和绿色点是一样的,但是很明显它到蓝色点的距离更远啊,你怎么能忽略他呢?
哎,这会在枚举下一条边的时候被算到啊……
所以说,前方的点,以后再来探索吧。
但是这样还有个问题,看下面这种情况:
当凸包退化为线段的时候,可以发现
t
t
t 指针永远也移动不了。
这只是特殊情况,放在一般情况下就是一开始选的边不会记录到答案里面。
解决方案就是令 t = 2 t=2 t=2,这样不会遗漏任何一种情况。
还有,如果你写了 ≤ \leq ≤,那就会一直死循环,嗨嗨嗨。
完整代码
对应题目:P2742。
#include<bits/stdc++.h>
using namespace std;
int n,st[50002],top,xia,t=1;
bool vis[50001];
struct node{
long long x,y;
bool operator<(const node&cwx)const{return x!=cwx.x?x<cwx.x:y<cwx.y;}
node operator-(const node&cwx)const{return {x-cwx.x,y-cwx.y};}
long long operator*(const node&cwx)const{return x*cwx.y-y*cwx.x;}
long long dis2(){return x*x+y*y;}
}p[50001];
long long ans;
inline long long dis(node a,node b,node c){
long long tmp=(a-c)*(b-c);
return abs(tmp*tmp/(a-b).dis2());
}
int main(){
scanf("%d",&n),t=2;
for(int i=1;i<=n;++i)scanf("%lld%lld",&p[i].x,&p[i].y);
sort(p+1,p+n+1),st[++top]=1;
for(int i=2;i<=n;++i){
while(top>1&&(p[st[top]]-p[st[top-1]])*(p[i]-p[st[top]])<=0)vis[st[top--]]=0;
vis[st[++top]=i]=1;
}xia=top;
for(int i=n-1;i;--i)
if(!vis[i]){
while(top>xia&&(p[st[top]]-p[st[top-1]])*(p[i]-p[st[top]])<=0)vis[st[top--]]=0;
vis[st[++top]=i]=1;
}
for(int i=1;i<top;++i){
while(dis(p[st[i]],p[st[i+1]],p[st[t]])<dis(p[st[i]],p[st[i+1]],p[st[t+1]]))t=t%(top-1)+1;
ans=max(ans,max((p[st[i]]-p[st[t]]).dis2(),(p[st[i+1]]-p[st[t]]).dis2()));
}return 0&printf("%lld",ans);
}
继续批话
本来想再放一道典中典例题的,但是已经九点二十了,马上回寝室,只能草草收尾。