【BZOJ】1069: [SCOI2007]最大土地面积(凸包+旋转卡壳)

时间:2022-01-21 16:48:48

http://www.lydsy.com/JudgeOnline/problem.php?id=1069

显然这四个点在凸包上,然后枚举两个点找上下最大的三角形即可。

找三角形表示只想到三分QAQ。。。。。。。

看了题解发现。。。。这是单调的。。。。直接扫。。。然后,这货叫“旋转卡壳”?是qia还是ka。。。。

自己一开始写的wa了。。。。。。。。然后照标程写,,又wa了。。。。

后来发现是凸包写渣了QAQ。。。自己原来的找三角形是没问题的。。。我好sb。。。。。

可是,为嘛凸包这样写会错。

cross(b[cnt], a[i], b[cnt-1])<=0变成cross(a[i], b[cnt], b[cnt-1])>0为嘛会错!!!!!!!

//upd:^(*&%*^^我为什么那么逗。。。。不应该是>=0吗。。。。。。。。。。。。。

2015.4.20 upd:

那是因为数据中有重复点!去掉重复点就能a啦!

如果凸包求的是允许三点共线的,那么数据中不能出现重复点!(因为这样会使得求凸壳出错!


两个写法:

自己yy的(凸包那里改了)

#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>
#include <set>
#include <map>
using namespace std;
typedef long long ll;
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << (#x) << " = " << (x) << endl
#define error(x) (!(x)?puts("error"):0)
#define rdm(x, i) for(int i=ihead[x]; i; i=e[i].next)
inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<'0'||c>'9'; c=getchar()) if(c=='-') k=-1; for(; c>='0'&&c<='9'; c=getchar()) r=r*10+c-'0'; return k*r; } const int N=2005;
struct dat { double x, y; }a[N], b[N];
int n, cnt;
double suml[N][N], sumr[N][N];
bool cmp(const dat &a, const dat &b) { return a.x<b.x; }
double cross(const dat &a, const dat &b, const dat &c) {
static dat x, y;
x.x=a.x-c.x; x.y=a.y-c.y;
y.x=b.x-c.x; y.y=b.y-c.y;
return x.x*y.y-y.x*x.y;
}
void tb() {
sort(a+1, a+1+n, cmp);
for1(i, 1, n) {
while(cnt>1 && cross(b[cnt], a[i], b[cnt-1])<=0) --cnt;
b[++cnt]=a[i];
}
int k=cnt;
for3(i, n-1, 1) {
while(cnt>k && cross(b[cnt], a[i], b[cnt-1])<=0) --cnt;
b[++cnt]=a[i];
}
if(n>1) --cnt;
}
double S(const dat &a, const dat &b, const dat &c) { return abs(cross(a, b, c))/2.0; }
void work() {
int mid, nxt, j;
for1(i, 1, n) {
mid=(i%n)+1;
for1(t, 1, n-2) {
j=(i+t)%n+1;
nxt=(mid%n)+1;
while(nxt!=j && S(b[i], b[j], b[nxt])>=S(b[i], b[j], b[mid])) { mid=nxt; nxt=(mid%n)+1; }
suml[i][j]=S(b[i], b[j], b[mid]); //printf("%d %d chose:%d: S:%.3f\n", i, j, mid, suml[i][j]);
}
mid=(i-2+n)%n+1;
for1(t, 1, n-2) {
j=(i-2+n-t+n)%n+1;
nxt=(mid-2+n)%n+1;
while(nxt!=j && S(b[i], b[j], b[nxt])>=S(b[i], b[j], b[mid])) { mid=nxt; nxt=((mid-2+n)%n)+1; }
sumr[i][j]=S(b[i], b[j], b[mid]);
}
}
} int main() {
read(n);
for1(i, 1, n) scanf("%lf%lf", &a[i].x, &a[i].y);
tb(); n=cnt;
//for1(i, 1, n) printf("%.2f %.2f\n", b[i].x, b[i].y);
work(); double ans=0;
for1(i, 1, n) for1(j, 1, n) if(suml[i][j] && sumr[i][j]) ans=max(ans, suml[i][j]+sumr[i][j]);
printf("%.3f", ans);
return 0;
}

照hzwer神犇的:

#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>
#include <set>
#include <map>
using namespace std;
typedef long long ll;
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << (#x) << " = " << (x) << endl
#define error(x) (!(x)?puts("error"):0)
#define rdm(x, i) for(int i=ihead[x]; i; i=e[i].next)
inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<'0'||c>'9'; c=getchar()) if(c=='-') k=-1; for(; c>='0'&&c<='9'; c=getchar()) r=r*10+c-'0'; return k*r; } const int N=2005;
struct dat { double x, y; }a[N], b[N];
int n, cnt;
bool cmp(const dat &a, const dat &b) { return a.x<b.x; }
double cross(const dat &a, const dat &b, const dat &c) {
static dat x, y;
x.x=a.x-c.x; x.y=a.y-c.y;
y.x=b.x-c.x; y.y=b.y-c.y;
return x.x*y.y-y.x*x.y;
}
void tb() {
sort(a+1, a+1+n, cmp);
for1(i, 1, n) {
while(cnt>1 && cross(b[cnt], a[i], b[cnt-1])<=0) --cnt;
b[++cnt]=a[i];
}
int k=cnt;
for3(i, n-1, 1) {
while(cnt>k && cross(b[cnt], a[i], b[cnt-1])<=0) --cnt;
b[++cnt]=a[i];
}
if(n>1) --cnt;
}
double S(const dat &a, const dat &b, const dat &c) { return abs(cross(a, b, c)); }
double work() {
b[cnt+1]=b[1];
int dmid, umid;
double ret=0;
for1(i, 1, n) {
dmid=(i%n)+1, umid=((i+2)%n)+1;
for1(j, i+2, n) {
while((dmid%n)+1!=j && S(b[i], b[j], b[dmid+1])>=S(b[i], b[j], b[dmid])) dmid=(dmid%n)+1;
while((umid%n)+1!=i && S(b[i], b[j], b[umid+1])>=S(b[i], b[j], b[umid])) umid=(umid%n)+1;
ret=max(ret, S(b[i], b[j], b[dmid])+S(b[i], b[j], b[umid]));//printf("%d %d chose:%d: S:%.3f\n", i, j, mid, suml[i][j]);
}
}
return ret/2.0;
} int main() {
read(n);
for1(i, 1, n) scanf("%lf%lf", &a[i].x, &a[i].y);
tb(); n=cnt;
//for1(i, 1, n) printf("%.2f %.2f\n", b[i].x, b[i].y);
printf("%.3f", work());
return 0;
}

  


Description

在某块平面土地上有N个点,你可以选择其中的任意四个点,将这片土地围起来,当然,你希望这四个点围成的多边形面积最大。

Input

第1行一个正整数N,接下来N行,每行2个数x,y,表示该点的横坐标和纵坐标。

Output

最大的多边形面积,答案精确到小数点后3位。

Sample Input

5
0 0
1 0
1 1
0 1
0.5 0.5

Sample Output

1.000

HINT

数据范围 n<=2000, |x|,|y|<=100000

Source