• HDU 3934


    /*这是用的有旋转卡壳的思想。
    首先确定i,j,对k进行循环,知道找到第一个k使得cross(i,j,k)>cross(i,j,k+1),如果k==i进入下一次循环。
    对j,k进行旋转,每次循环之前更新最大值,然后固定一个j,同样找到一个k使得cross(i,j,k)>cross(i,j,k+1)。对j进行++操作,继续进行下一次,
    知道j==k为止。
    */ 
    
    #include <iostream>
    #include <cstdio>
    #include <algorithm>
    #include <cmath>
    
    using namespace std;
    
    struct point{
    	double x,y;
    }p[1000100];
    int n;
    
    int ans[1000100],st[1000100],cnt,stop;
    
    bool cmp(point A, point B){
    	if(A.y<B.y) return true;
    	else if(A.y==B.y){
    		if(A.x<B.x) return true;
    	}
    	return false;
    }
    
    double multi(point a,point b,point c){
    	point p1; p1.x=a.x-c.x; p1.y=a.y-c.y;
    	point p2; p2.x=b.x-c.x; p2.y=b.y-c.y;
    	return p1.x*p2.y-p1.y*p2.x;
    }
    
    void slove(){
    	cnt=stop=0;
    	st[stop++]=0; st[stop++]=1;
    	for(int i=2;i<n;i++){
    		while(stop>1&&multi(p[i],p[st[stop-1]],p[st[stop-2]])>=0) stop--;
    		st[stop++]=i;
    	}
    	for(int i=0;i<stop;i++)
    	ans[cnt++]=st[i];
    	stop=0; st[stop++]=n-1; st[stop++]=n-2;
    	for(int i=n-3;i>=0;i--){
    		while(stop>1&&multi(p[i],p[st[stop-1]],p[st[stop-2]])>=0) stop--;
    		st[stop++]=i;
    	}
    	for(int i=1;i<stop-1;i++)
    	ans[cnt++]=st[i];
    /*	for(int i=0;i<cnt;i++)
    	cout<<ans[i]<<endl;
    	cout<<endl;*/
    }
    
    double Triangle(point a,point b,point c){
    	point p1; p1.x=a.x-c.x; p1.y=a.y-c.y;
    	point p2; p2.x=b.x-c.x; p2.y=b.y-c.y;
    	return fabs((p1.x*p2.y-p1.y*p2.x)*1.0)/2.0;
    }
    
    double Area(){
    	int q; int j;
    	double anst=0;
    	for(int i=0;i<cnt;i++){
    		j=(i+1)%cnt;
    		q=(j+1)%cnt;
    		while(Triangle(p[ans[i]],p[ans[j]],p[ans[q]])<=Triangle(p[ans[i]],p[ans[j]],p[ans[(q+1)%cnt]])&&q!=i)
    		q=(q+1)%cnt;    //枚举了当前最远的K点 
    		anst=max(anst,Triangle(p[ans[i]],p[ans[j]],p[ans[q]]));
    		if(q==i) continue;
    		while(j!=i&&q!=i){
    			anst=max(anst,Triangle(p[ans[i]],p[ans[j]],p[ans[q]]));
    			while(Triangle(p[ans[i]],p[ans[j]],p[ans[q]])<=Triangle(p[ans[i]],p[ans[j]],p[ans[(q+1)%cnt]])&&q!=i)
    			q=(q+1)%cnt;
    			j=(j+1)%cnt;
    		}
    		
    	}
    	return anst;
    }
    
    int main(){
    	while(scanf("%d",&n)!=EOF){
    	//	if(n==-1) break;
    		for(int i=0;i<n;i++){
    			scanf("%lf%lf",&p[i].x,&p[i].y);
    		}
    		sort(p,p+n,cmp);
    		slove();
    		double anst=0;
    		anst=max(anst,Area());
    		printf("%.2lf
    ",anst);
    	}
    	return 0;
    }
    

      

  • 相关阅读:
    Ubuntu下speedtest的安装
    Ubuntu下快速安装LAMP server
    用python3.3爬取图片
    用python3.3结合snmpget截取信息
    开源数据库Postgresql的备份和恢复
    vuedraggable 两个不同结构的数组之间拖拽
    MyBatis 使用另一个 mapper 中的 resultMap 和 sql
    怎样处理从后台传来大量数据的分页操作
    doc命令操纵数据库
    (全)Mac中mysql的root账户初始密码忘记解决办法以及更改
  • 原文地址:https://www.cnblogs.com/jie-dcai/p/3891337.html
Copyright © 2020-2023  润新知