首页 > 代码库 > 【BZOJ】1069: [SCOI2007]最大土地面积(凸包+旋转卡壳)
【BZOJ】1069: [SCOI2007]最大土地面积(凸包+旋转卡壳)
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吗。。。。。。。。。。。。。
两个写法:
自己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
0 0
1 0
1 1
0 1
0.5 0.5
Sample Output
1.000
HINT
数据范围 n<=2000, |x|,|y|<=100000
Source
【BZOJ】1069: [SCOI2007]最大土地面积(凸包+旋转卡壳)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。