POJ1264 SCUD Busters 凸包

POJ1264

有m个国家(m<=20)对每个国家给定n个城镇 这个国家的围墙是保证围住n个城镇的周长最短的多边形 必然是凸包

进行若干次导弹发射 落到一个国家内则国家被破坏

最后回答总共有多少面积被破坏

首先求凸包

然后判断点是否在凸包内 要用O(logn)的判断方法 不然会超时

这道题常数卡的有点紧 TLE三次才过

蒟蒻没救了 2017年做1991年的题还被卡常数

#include<iostream>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;

const double eps=1e-9;

int cmp(double x)
{
 if(fabs(x)<eps)return 0;
 if(x>0)return 1;
 	else return -1;
}

const double pi=acos(-1.0);

inline double sqr(double x)
{
 return x*x;
}

struct point
{
 double x,y;
 point (){}
 point (double a,double b):x(a),y(b){}
 void input()
 	{
 	 scanf("%lf%lf",&x,&y);
	}
 friend point operator +(const point &a,const point &b)
 	{
 	 return point(a.x+b.x,a.y+b.y);
	}
 friend point operator -(const point &a,const point &b)
 	{
 	 return point(a.x-b.x,a.y-b.y);
	}
 friend bool operator ==(const point &a,const point &b)
 	{
 	 return cmp(a.x-b.x)==0&&cmp(a.y-b.y)==0;
	}
 friend point operator *(const point &a,const double &b)
 	{
 	 return point(a.x*b,a.y*b);
	}
 friend point operator*(const double &a,const point &b)
 	{
 	 return point(a*b.x,a*b.y);
	}
 friend point operator /(const point &a,const double &b)
 	{
 	 return point(a.x/b,a.y/b);
	}
 double norm()
 	{
 	 return sqrt(sqr(x)+sqr(y));
	}
};

struct line
{
 point a,b;
 line(){};
 line(point x,point y):a(x),b(y)
 {

 }
};
double det(const point &a,const point &b)
{
 return a.x*b.y-a.y*b.x;
}

double dot(const point &a,const point &b)
{
 return a.x*b.x+a.y*b.y;
}

double dist(const point &a,const point &b)
{
 return (a-b).norm();
}

point rotate_point(const point &p,double A)
{
 double tx=p.x,ty=p.y;
 return point(tx*cos(A)-ty*sin(A),tx*sin(A)+ty*cos(A));
}

bool parallel(line a,line b)
{
 return !cmp(det(a.a-a.b,b.a-b.b));
}

bool line_joined(line a,line b,point &res)
{
 if(parallel(a,b))return false;
 double s1=det(a.a-b.a,b.b-b.a);
 double s2=det(a.b-b.a,b.b-b.a);
 res=(s1*a.b-s2*a.a)/(s1-s2);
 return true;
}

bool pointonSegment(point p,point s,point t)
{
 return cmp(det(p-s,t-s))==0&&cmp(dot(p-s,p-t))<=0;
}

void PointProjLine(const point p,const point s,const point t,point &cp)
{
 double r=dot((t-s),(p-s))/dot(t-s,t-s);
 cp=s+r*(t-s);
}

struct polygon_convex
{
 vector<point>P;
 polygon_convex(int Size=0)
 	{
 	 P.resize(Size);
	}
};

bool comp_less(const point &a,const point &b)
{
 return cmp(a.x-b.x)<0||cmp(a.x-b.x)==0&&cmp(a.y-b.y)<0;

}

polygon_convex convex_hull(vector<point> a)
{
 polygon_convex res(2*a.size()+5);
 sort(a.begin(),a.end(),comp_less);
 a.erase(unique(a.begin(),a.end()),a.end());//删去重复点
 int m=0;
 for(int i=0;i<a.size();i++)
 	{
 	 while(m>1&&cmp(det(res.P[m-1]-res.P[m-2],a[i]-res.P[m-2]))<=0)--m;
 	 res.P[m++]=a[i];
	}
 int k=m;
 for(int i=int(a.size())-2;i>=0;--i)
 	{
 	 while(m>k&&cmp(det(res.P[m-1]-res.P[m-2],a[i]-res.P[m-2]))<=0)--m;
 	 res.P[m++]=a[i];
	}
 res.P.resize(m);
 if(a.size()>1)res.P.resize(m-1);
 return res;
}

bool is_convex(vector<point> &a)
{
 for(int i=0;i<a.size();i++)
 	{
 	 int i1=(i+1)%int(a.size());
 	 int i2=(i+2)%int(a.size());
 	 int i3=(i+3)%int(a.size());
 	 if((cmp(det(a[i1]-a[i],a[i2]-a[i1]))*cmp(det(a[i2]-a[i1],a[i3]-a[i2])))<0)
	  	return false;
	}
 return true;
}
int containO(const polygon_convex &a,const point &b)
{
 int n=a.P.size();
 point g=(a.P[0]+a.P[n/3]+a.P[2*n/3])/3.0;
 int l=0,r=n;
 while(l+1<r)
 	{
 	 int mid=(l+r)/2;
 	 if(cmp(det(a.P[l]-g,a.P[mid]-g))>0)
 	 	{
 	 	 if(cmp(det(a.P[l]-g,b-g))>=0&&cmp(det(a.P[mid]-g,b-g))<0)r=mid;
 	 	 	else l=mid;
		}else
			{
			 if(cmp(det(a.P[l]-g,b-g))<0&&cmp(det(a.P[mid]-g,b-g))>=0)l=mid;
 	 	 		else r=mid;
			}
	}
 r%=n;
 int z=cmp(det(a.P[r]-b,a.P[l]-b))-1;
 if(z==-2)return 1;
 return z;
}

polygon_convex pc[30];
double area(int n)
{
 double ans=0;
 vector<point>a;
 a.clear();
 for(int i=0;i<pc[n].P.size();i++)
    a.push_back(pc[n].P[i]);
 a.push_back(a[0]);
 for(int i=0;i<a.size()-1;i++)ans+=det(a[i+1],a[i]);
 //cout<<ans/2<<endl;
 return ans/2.0;
}

int tot;
double are[30];

vector<point> pp;
int shoot[600][600];
bool damed[30];
int main()
{//freopen("t.txt","r",stdin);
 int n;
 tot=1;
 while(scanf("%d",&n)&&n!=-1)
 	{
	 pp.clear();
 	 for(int i=0;i<n;i++)
 	 	{
 	 	 point p;
 	 	 p.input();
 	 	 pp.push_back(p);
		}
	 pc[tot]=convex_hull(pp);
	 //cout<<pc[tot].P.size()<<endl;
	 are[tot]=area(tot);
	 //cout<<are[tot]<<endl;
	 tot++;
	}
 int x,y;
 memset(damed,0,sizeof(damed));
 int sum=0;
 memset(shoot,0,sizeof(shoot));
 while(scanf("%d%d",&x,&y)!=EOF)
 	{
 	 if(sum==tot-1)break;
 	 if(shoot[x][y]==0)
 	 	{
 	 	 for(int i=1;i<tot;i++)
		   	if(containO(pc[i],point(x,y)))
			   	{
			   	 shoot[x][y]=i;
			   	 break;
				}
		}
 	 if(!damed[shoot[x][y]]){sum++;damed[shoot[x][y]]=1;}
	}
 double tarea=0;
 for(int i=1;i<tot;i++)
 	{
 	 //cout<<are[i]<<endl;
 	 if(damed[i])tarea+=are[i];
	}
 printf("%.2lf\n",-tarea+0.0005);
 return 0;
}

  

时间: 2024-10-20 10:08:22

POJ1264 SCUD Busters 凸包的相关文章

POJ3528 HDU3662 三维凸包模板

POJ3528 HDU3662 第一道题 给定若干点 求凸包的表面积,第二题 给定若干点就凸包的面数. 简单说一下三维凸包的求法,首先对于4个点假设不共面,确定了唯一四面体,对于一个新的点,若它不在四面体内,为了让它进入凸包, 则对于所有凸包上的边,若边的一面是该点可以看到的而另一面看不到,则该点与该边构成的面要加入凸包. 模板代码非常清晰, #include<stdio.h> #include<algorithm> #include<string.h> #includ

关于2016.12.12——T1的反思:凸包的意义与应用

2016.12.12 T1 给n个圆,保证圆圆相离,求将圆围起来的最小周长.n<=100 就像上图.考场上,我就想用切线的角度来做凸包.以圆心x,y排序,像点凸包一样,不过用两圆之间的下切线角度来判断. 这就是下切线(我自己瞎编的名字): 好像是对的啊: 然后我就保证必AC的希望,用这种写法交了,然后就只得了N=2的暴力分... 自以为是正解,却落得如此下场... 为什么?这样不对吗?借用学长的力量,果然被Hack掉了: 这种情况,圆心排序后,检测的顺序并不是圆上的切点的顺序,自然就会挂. 蓝瘦

poj 1113 凸包周长

Wall Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 33888   Accepted: 11544 Description Once upon a time there was a greedy King who ordered his chief Architect to build a wall around the King's castle. The King was so greedy, that he w

UVA 11800 Determine the Shape --凸包第一题

题意: 给四个点,判断四边形的形状.可能是正方形,矩形,菱形,平行四边形,梯形或普通四边形. 解法: 开始还在纠结怎么将四个点按序排好,如果直接处理的话,有点麻烦,原来凸包就可搞,直接求个凸包,然后点就自动按逆时针排好了,然后就判断就可以了,判断依据题目下面有,主要是用到点积和叉积,判断垂直用点积,判断平行用叉积. 代码: #include <iostream> #include <cstdio> #include <cstring> #include <cstd

POJ 3348 最直接的凸包问题

题目大意: 给定一堆树的点,找到能组合成的最大面积,一个物体占50面积,求最多放多少物体 1 #include <cstdio> 2 #include <cstring> 3 #include <iostream> 4 #include <algorithm> 5 using namespace std; 6 const double eps = 1e-10; 7 const int N = 10005; 8 double myabs(double x) 9

hdu 1348 Wall(凸包模板题)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1348 Wall Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Submission(s): 3386    Accepted Submission(s): 968 Problem Description Once upon a time there was a gre

POJ1113 Wall【凸包】

Wall Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 24604   Accepted: 8183 Description King想在自己的n个城堡外建Wall,使Wall与任一城堡距离至少为L且能围住它的城堡. 求Wall最短长度. Input 第一行为 N (3 <= N <= 1000) 和 L(1 <= L <= 1000). 接下来 N 行,每行为城堡坐标Xi,Yi (-10000 <=

[凸包]Triangles

https://nanti.jisuanke.com/t/15429 题目大意:给出平面内$n$个整数坐标点,保证无三点共线.可以进行若干次连线,每次选择一个点对连接线段,但是任意两条线段都不得在给定的$n$个点之外有交点.问连线完成后,最多能构造出多少个三角形. 解题关键: 小于三个点的情况答案为零.考虑三个点的情况,由于三点不共线,必然构成一个三角形.现加入第四个点,若其在原三角形外部,则称其为外点,可以新构造$1$个三角形:若其在原三角形内部,则称其为内点,可以新构造$3$个三角形.故要尽

ZOJ - 3537 Cake (凸包+区间DP+最优三角剖分)

Description You want to hold a party. Here's a polygon-shaped cake on the table. You'd like to cut the cake into several triangle-shaped parts for the invited comers. You have a knife to cut. The trace of each cut is a line segment, whose two endpoin