旋转卡壳简介(POJ2187)(洛谷P1452)

标签: POJ 旋转卡壳 洛谷
16人阅读 评论(0) 收藏 举报
分类:

旋转卡壳简介(POJ2187)

读音

word上是这么读的:
这里写图片描述

前置技能

凸包

算法用途

旋转卡壳可以在O(n)的时间内确定一对对踵点对,它的用途包括但不限于:计算距离(凸多边形直径)、计算外接矩形(最小面积/周长)、三角剖分(洋葱三角剖分)等其他奇奇怪怪的东西。

算法实现(模板)

POJ2187洛谷P1452)为例。

这道题要我们求所有点之间的最大距离。我们先求出这些点的凸包,最大距离肯定在凸包的对踵点对上。那么我们枚举所有的对踵点对就好了,而旋转卡壳就是实现这个的。

对于每个点,我们枚举它的对踵点对。假设当前点的为i,上一个点的对踵点对为x,那么当(x+1x)×(ix)(x+1x)×(i+1x)时,说明x仍不是i的对踵点对(因为i+1i更符合)。那么把x+1并继续验证。当找到对踵点对时更新答案就好了。

代码实现(旋转卡壳为RC()):

#include<cmath>
#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 50005
#define F friend
#define eps 1e-7
#define sqr(x) ((x)*(x))
#define abs(x) ((x)>0?(x):-(x))
using namespace std;
struct P{
    int x,y;
    F bool operator < (P a,P b){
        return abs(a.y-b.y)<eps?a.x<b.x:a.y<b.y;
    }
    F int cs(P a,P b,P c){//叉积
        return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);
    }
    F int dis(P a,P b){
        return sqr(a.x-b.x)+sqr(a.y-b.y);
    }
}t[N],s[N];
int n,tp,ans;
inline char readc(){
    static char buf[100000],*l=buf,*r=buf;
    if (l==r) r=(l=buf)+fread(buf,1,100000,stdin);
    if (l==r) return EOF; return *l++;
}
inline int _read(){
    int x=0,f=1; char ch=readc();
    while (!isdigit(ch)) { if (ch=='-') f=-1; ch=readc(); }
    while (isdigit(ch)) x=(x<<3)+(x<<1)+(ch^48),ch=readc();
    return x*f;
}
inline bool cmp(P a,P b){
    int f=cs(t[1],a,b); if (f) return f>0;
    return dis(t[1],a)-dis(t[1],b)<eps*eps;
}
inline void GR(){//凸包
    sort(t+1,t+n+1),s[++tp]=t[1],sort(t+2,t+n+1,cmp);
    for (int i=2;i<=n;s[++tp]=t[i++])
        while (tp>1&&cs(s[tp-1],s[tp],t[i])<eps) tp--;
}
inline void RC(){//旋转卡壳
    s[0]=s[tp]; int x=1;
    for (int i=0;i<tp;ans=max(ans,dis(s[i++],s[x])))
        while (cs(s[x],s[x+1],s[i])-cs(s[x],s[x+1],s[i+1])>eps)
            x=(x+1)%tp;
}
int main(){
    n=_read();
    for (int i=1;i<=n;i++)
        t[i].x=_read(),t[i].y=_read();
    return GR(),RC(),printf("%d\n",ans),0;
}
查看评论

poj2187 最远点对(旋转卡壳)

Bessie, Farmer John’s prize cow, has just won first place in a bovine beauty contest, earning the ti...
  • zjy2015302395
  • zjy2015302395
  • 2017-03-29 02:15:08
  • 368

计算几何——旋转卡壳【POJ2187】

裸题,只有代码
  • wyh0410
  • wyh0410
  • 2016-09-18 11:46:05
  • 121

旋转卡壳 模板 【poj2187】Beauty Contest

旋转卡壳,kake?qiake?kaqiao?qiaqiao? 一番研究过后,发现旋转卡(qia)壳(ke)是比较正确的读法(也是我最难受的读法QAQ) 旋转卡壳是在凸包上求两个距离最远的点。 ...
  • Todobe
  • Todobe
  • 2017-02-25 21:42:27
  • 394

POJ2187【旋转卡壳】

用向量的外积来代替斜率.可以有效减小误差. /* I will wait for you */ #include #include #include #include #include #incl...
  • Lethelody
  • Lethelody
  • 2015-04-01 01:09:58
  • 234

【OI康复训练】POJ2187旋转卡壳

随便写一个旋转卡壳,结果有重复点实在是日了狗了 #include #include #include #include using namespace std; struct Point { i...
  • johann_oier
  • johann_oier
  • 2015-07-10 17:46:47
  • 584

【POJ2187】Beauty Contest【旋转卡壳】

【题目链接】 学习一发旋转卡壳。 /* Telekinetic Forest Guard */ #include #include #include using namespace...
  • BraketBN
  • BraketBN
  • 2016-05-26 15:16:22
  • 171

POJ2187 Beauty Contest(旋转卡壳)

题目大意:凸包最远点对 思路:模板套,初写很粗糙。好几个相同函数拼凑成的,写的很丑陋。。 program: #include #include #include #include  ...
  • kg_second
  • kg_second
  • 2012-09-21 12:27:48
  • 424

【旋转卡壳】Beauty Contest POJ2187

Beauty Contest Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 23413 ...
  • jiangzh7
  • jiangzh7
  • 2013-04-14 12:55:30
  • 448

POJ2187 Andrew凸包+旋转卡壳

浏览(水了水)网上各种凸包和旋转卡壳的讲解,然后来水这道模板题,结果被坑死了。。。。。。。 我的模板里是double类型的,放到POJ上狂WA(输出用的%.lf,末尾没0啊),于是很愤慨地改成了整型...
  • nlj1999
  • nlj1999
  • 2016-01-09 16:44:17
  • 242

poj2187 Beauty Contest(旋转卡壳)

模板题 注意:避免凸包上三点共线 代码 #include #include #define eps 1e-10 struct Point { double x,y; Poin...
  • cjk_cjk
  • cjk_cjk
  • 2015-06-15 15:16:50
  • 263
    个人资料
    专栏达人 持之以恒
    等级:
    访问量: 5万+
    积分: 3539
    排名: 1万+
    博客专栏
    文章分类