zoukankan      html  css  js  c++  java
  • bzoj1069[SCOI2007]最大土地面积

    传送门

    这个题我改了好久呢。(最终居然是挂在常识上,我的算法完全没写错)
    算是计算几何的第一步啦
    这个题很容易想到四边形的四个点在凸包上,但是暴力枚举复杂度依然感人
    可以考虑先枚举对角线,那么另外两个点就满足单峰,可以用三分优化,据说这样写卡卡常能过
    然后发现枚举对角线的过程中,另外两个点也是跟着在动的,维护一下就好了
    时间复杂度(O(n^2))
    代码:(懒得判叉积方向了,直接用fabs,反正是面积,符号不影响)

    #include<cstdio>
    #include<iostream>
    #include<algorithm>
    #include<cmath>
    using namespace std;
    #define rg register
    const int maxn=2010;
    int n,top,top1;double ans;
    struct oo{double x,y;}a[maxn],st[maxn],stt[maxn];
    double operator*(oo a,oo b){return a.x*b.y-a.y*b.x;}
    oo operator-(oo a,oo b){return (oo){b.x-a.x,b.y-a.y};}
    bool cmp(oo a,oo b){return a.x<b.x;}
    int main()
    {
    	scanf("%d",&n);
    	for(rg int i=1;i<=n;i++)scanf("%lf%lf",&a[i].x,&a[i].y);
    	sort(a+1,a+n+1,cmp);
    	++top,st[top]=a[top],++top,st[top]=a[top];
    	for(rg int i=3;i<=n;i++)
    	{
    		while(top-1&&(st[top-1]-st[top])*(st[top-1]-a[i])<=0)top--;
    		st[++top]=a[i];
    	}
    	++top1,stt[top1]=a[top1],++top1,stt[top1]=a[top1];
    	for(rg int i=3;i<=n;i++)
    	{
    		while(top1-1&&(stt[top1-1]-stt[top1])*(stt[top1-1]-a[i])>=0)top1--;
    		stt[++top1]=a[i];
    	}
    	for(rg int i=top1-1;i>=2;i--)st[++top]=stt[i];st[top+1]=st[1];
    	for(rg int i=1;i<=top;i++)
    	{
    		int a=i%top+1,b=(i+2)%top+1;
    		for(rg int j=i+2;j<=top;j++)
    		{
    			while(a%top+1!=j&&fabs((st[a+1]-st[i])*(st[a+1]-st[j]))>fabs((st[a]-st[i])*(st[a]-st[j])))a=a%top+1;
    			while(b%top+1!=i&&fabs((st[b+1]-st[i])*(st[b+1]-st[j]))>fabs((st[b]-st[i])*(st[b]-st[j])))b=b%top+1;
    			ans=max(ans,fabs((st[a]-st[j])*(st[a]-st[i]))+fabs((st[b]-st[i])*(st[b]-st[j])));
    		}
    	}
    	printf("%.3lf",ans/2.0);
    }
    
  • 相关阅读:
    随机色块
    JQ命令汇总
    JQ选择器
    cookie
    tab切换
    Ajax跨域
    RocksDB介绍:一个比LevelDB更彪悍的引擎
    谷歌的诀窍:如何取消验证码
    Ruby on Rails创始人DHH谈如何进行混合移动APP开发
    SequoiaDB 架构指南
  • 原文地址:https://www.cnblogs.com/lcxer/p/10324343.html
Copyright © 2011-2022 走看看