HDU 4946 Area of Mushroom 凸包 第八次多校

题目链接:hdu 4946

题意:一大神有N个学生,各个都是小神,大神有个二次元空间,每一个小神都有一个初始坐标,如今大神把这些空间分给徒弟们,规则是假设这个地方有一个人比谁都先到这,那么这个地方就是他的,问谁的地盘面积是无穷大的.

思路:由于空间是无限大的,所以仅仅要速度最大就有可能有无限的地盘,重合的点不能严格的比对方快,也不符合规定,然后求速度最大的点围城的凸包,凡是在凸包上的点且不重合的,地盘都是无穷大.

比赛的时候,不会凸包的,现找的模板,不了解性能,不知道计算凸包不能传重合的点进去结果WA出翔

#include <cstdio>
#include <cstring>
#include <iostream>
#include <cmath>
#include <queue>
#include <vector>
#include <algorithm> #define INF 1e9
#define maxn 1000
#define EPS 1e-6 #define N 1000
using namespace std;
struct point
{
int x,y; //横纵坐标 : x,y
double len,theta; //与參考点的距离 len 与參考点构成的向量与 (1,0)向量构成的夹角的余弦值 theta
}g[N]; //定义了一个全局变量,记录凸包中的点
/*--------按余弦值,从大到小高速排序--------*/
void qsort(int st,int en)
{
int i=st,j=en;
g[0]=g[i];
while(i<j)
{
while(i<j && g[0].theta>=g[j].theta) j--;
if(i<j) { g[i]=g[j]; i++; }
while(i<j && g[0].theta<=g[i].theta) i++;
if(i<j) { g[j]=g[i]; j--; }
}
g[i]=g[0];
if(st<i-1) qsort(st,i-1);
if(i+1<en) qsort(i+1,en);
} /*-----------Graham 扫描法-------------*/
void graham(int *n)
{
/*第一步,寻找y坐标最小,然后x坐标最小的点*/
int p=1;
for(int i=2;i<=*n;i++)
if((g[i].y<g[p].y)||(g[i].y==g[p].y && g[i].x<g[p].x)) p=i;
g[0]=g[p]; g[p]=g[1]; g[1]=g[0];
/*找到该点,并把它存放在 g 中的第一个元素的位子上*/ /*第二步,计算全部的点距离參考点的距离(len) 还有夹角的余弦值 (theta)*/
for(int i=2;i<=*n;i++)
{
g[i].len=sqrt((g[i].x-g[1].x)*(g[i].x-g[1].x)+(g[i].y-g[1].y)*(g[i].y-g[1].y));
g[i].theta=100*(g[i].x-g[1].x)/g[i].len;
}
qsort(2,*n);//先依据夹角的余弦值从大到小排序 /*第三步,将全部theta值相等的点,仅仅保存len值最大的,存放在数组map中*/
point map[N];
int tot=0; p=1;
while(p<=*n)
{
int k=p;
while(fabs(g[p].theta-g[p+1].theta)<=EPS)
{
if(g[p+1].len>g[k].len) k=p+1;
p++;
}
map[++tot]=g[k];
p++;
} /*第四步,对map中的元素扫描一遍,确定凸包的元素,放在数组g中*/
*n=tot; tot=3; //先做了一个小小的处理,使得自己更好理解
memset(g,0,sizeof(g));
g[1]=map[1]; g[2]=map[2]; g[3]=map[3]; //先将前三个点入栈 g
for(int i=4;i<=*n;i++) //依次用map中的每一个点对g中的点进行一次推断,看是否是属于凸包
{
double chaji=(g[tot].x-g[tot-1].x)*(map[i].y-g[tot].y)-(map[i].x-g[tot].x)*(g[tot].y-g[tot-1].y);
for(;chaji<=0 && tot>=1;) //假设旋转的方向不同,g[tot]这个点就不是,删除,并继续推断 g 中下一个点是不是
{
tot--;
chaji=(g[tot].x-g[tot-1].x)*(map[i].y-g[tot].y)-(map[i].x-g[tot].x)*(g[tot].y-g[tot-1].y);
}
g[++tot]=map[i]; //将map[i]这个点入栈,至于是否是属于凸包中的点,等待以后的点来推断
}
*n=tot;//凸包处理完,总共同拥有tot个凸包上的点
} /*---------------------------------------------------------------------------------------------*/ struct Node{
int x,y,v,rank,bl;
}tmp;
int maxe;
vector<Node> stu;
int cas=0;
bool cmp(Node a,Node b)
{
return a.rank<b.rank;
}
void init(int n)
{
stu.clear();
maxe=0;
for(int i=0;i<n;i++)
{
scanf("%d %d %d",&tmp.x,&tmp.y,&tmp.v);
tmp.rank=i;
tmp.bl=0;
stu.push_back(tmp);
maxe=max(tmp.v,maxe);
}
}
bool check(int a,int b,int tt)
{
int c;
if(b==tt)
{
c=1;
}else
{
c=b+1;
}
double chaji=(g[c].x-g[b].x)*(stu[a].y-g[c].y)-(stu[a].x-g[c].x)*(g[c].y-g[b].y);
if(chaji<=EPS) return 1;
else return 0;
} int main()
{
//freopen("data.in","r",stdin);
int n;
while (~scanf("%d",&n))
{
if(n==0) break;
printf("Case #%d: ",++cas);
init(n); if(maxe==0)
{
for(int i=0;i<n;i++)
{
printf("0");
}
}
else
{
for(int i=0;i<n;i++)
{
if(stu[i].v==maxe)
{
stu[i].bl=1;
}
}
for(int i=0;i<n-1;i++)
{
if(stu[i].bl==1)//错误写法: if(stu[i].v=maxe)
{
for(int j=i+1;j<n;j++)
{
if(stu[i].x==stu[j].x&&stu[i].y==stu[j].y&&stu[i].v==stu[j].v&&i!=j)
{
stu[j].bl=0;
stu[i].bl=-1;
}
}
}
}
int tt=0;
for(int i=0;i<n;i++)
{
if(stu[i].v==maxe&&stu[i].bl!=0)
{ g[++tt].x=stu[i].x;
g[tt].y=stu[i].y;
}
}
//printf("%d\n",tt);
if(tt>=3)
{
graham(&tt);
for(int i=0;i<n;i++){
if(stu[i].bl==1){
stu[i].bl=0;
for(int j=1;j<=tt;j++)
{
if((stu[i].x==g[j].x&&stu[i].y==g[j].y)||check(i,j,tt))
{
stu[i].bl=1;
break;
}
}
}
}
}
sort(stu.begin(),stu.end(),cmp);
for(int i=0;i<n;i++)
{
if(stu[i].bl==-1)
stu[i].bl=0;
printf("%d",stu[i].bl);
}
} puts("");
}
return 0;
}
上一篇:Mybatis学习(8)逆向工程


下一篇:bzoj1071[SCOI2007]组队