[计算几何]旋转卡壳

欢迎大家访问我的老师的OJ———caioj.cn

题面描述

传送门

思路

凸包可以打SCY的。

去看wjyyy的卡壳吧

SCY的std有点问题的。

对于卡壳这一部分

while(mul(h[i+1],h[now],h[i])<mul(h[i+1],h[now+1],h[i]))

可以随意变号，就是不要加等于号！因为等于的东西实际上是等效的，所以不用加。加了反而会T。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
using namespace std;
const int N=5e4+10;
struct node{double x,y;}p[N],h[N];
double mul(node p1,node p2,node p0)
{
    double x1=p1.x-p0.x,y1=p1.y-p0.y;
    double x2=p2.x-p0.x,y2=p2.y-p0.y;
    return x1*y2-x2*y1; 
}
bool cmp(node a,node b){return a.x==b.x?a.y<b.y:a.x<b.x;}
double dis(node a,node b){return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);}
int sta[N],tp,n;
bool used[N];
void tb()
{
    sort(p+1,p+n+1,cmp);
    tp=0;
    sta[++tp]=1;
    for(int i=2;i<=n;i++)//下凸壳 
    {
        while(tp>1&&mul(p[sta[tp]],p[i],p[sta[tp-1]])<=0)used[sta[tp--]]=0;
        used[i]=1;
        sta[++tp]=i;
    }
    int top=tp;
    for(int i=n-1;i;i--)//上凸壳 
        if(!used[i])
        {
            while(tp>top&&mul(p[sta[tp]],p[i],p[sta[tp-1]])<=0)used[sta[tp--]]=0;
            sta[++tp]=i;
        }
    for(int i=1;i<=tp;i++)h[i]=p[sta[i]];
}
void kq()
{
    double ans=0;
    int now=2;
    for(int i=1;i<tp;i++)
    {
        while(mul(h[i+1],h[now],h[i])<mul(h[i+1],h[now+1],h[i]))
            {now++;if(now==tp)now=1;}
        ans=max(ans,max(dis(h[i+1],h[now]),dis(h[i],h[now])));
    }
    printf("%.0lf\n",ans);
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%lf%lf",&p[i].x,&p[i].y);
    tb();
    kq();
    return 0;
}