首页 > 代码库 > UVaLive4992:Jungle Outpost
UVaLive4992:Jungle Outpost
传送门
半平面交。
首先,由显然成立法可以证明炸连续的几个总比分散火力效果更佳。
所以二分答案,转化为判定问题,即间隔$ans$个点的连线的半平面交是否为空。
半平面交判定即可。
时间复杂度:$O(Nlog^2N)$
//UVaLive4992//by Cydiater//2017.2.1#include <iostream>#include <cstdio>#include <cstring>#include <string>#include <cmath>#include <ctime>#include <cstdlib>#include <queue>#include <map>#include <algorithm>#include <iomanip>#include <bitset>#include <set>#include <vector>using namespace std;#define ll long long#define up(i,j,n) for(int i=j;i<=n;i++)#define down(i,j,n) for(int i=j;i>=n;i--)#define cmax(a,b) a=max(a,b)#define cmin(a,b) a=min(a,b)#define db double#define Vector Pointconst int MAXN=1e5+5;const int oo=0x3f3f3f3f;const db eps=1e-10;inline int read(){ char ch=getchar();int x=0,f=1; while(ch>‘9‘||ch<‘0‘){if(ch==‘-‘)f=-1;ch=getchar();} while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();} return x*f;}int dcmp(db x){if(fabs(x)<eps)return 0;else return x<0?-1:1;}struct Point{ db x,y; Point(db x=0,db y=0):x(x),y(y){}};Vector operator + (Point x,Point y){return Vector(x.x+y.x,x.y+y.y);}Vector operator - (Point x,Point y){return Vector(x.x-y.x,x.y-y.y);}Vector operator * (Vector x,db p){return Vector(x.x*p,x.y*p);}Vector operator / (Vector x,db p){return Vector(x.x/p,x.y/p);}bool operator < (const Vector &x,const Vector &y){return dcmp(x.x-y.x)==0?x.y<y.y:x.x<y.x;}bool operator == (const Vector &x,const Vector &y){return dcmp(x.x-y.x)==0&&dcmp(x.y-y.y)==0;}struct Line{ Point P; Vector v; db ang; Line(){} Line(Point P,Vector v):P(P),v(v){ang=atan2(v.y,v.x);}};bool operator < (const Line &x,const Line &y){return x.ang<y.ang;}Point V[MAXN],P[MAXN];int N,top,head,tail;Line L[MAXN],q[MAXN];namespace solution{ Point Write(){db x=read(),y=read();return Point(x,y);} db Cross(Vector x,Vector y){return x.x*y.y-x.y*y.x;} bool Onleft(Point P,Line Li){return dcmp(Cross(P-Li.P,Li.v))<0;} void Prepare(){ V[0]=Write(); down(i,N-1,1)V[i]=Write(); } void Modify(int siz){ up(i,0,N-1)L[i]=Line(V[i],V[(i+siz)%N]-V[i]); } Point LineMeet(Line x,Line y){ db t=Cross(y.v,x.P-y.P)/Cross(x.v,y.v); return x.P+x.v*t; } bool check(){ head=1;tail=0;q[++tail]=L[0]; up(i,1,N-1){ while(head<tail&&!Onleft(P[tail-1],L[i]))tail--; while(head<tail&&!Onleft(P[head],L[i]))head++; q[++tail]=L[i]; if(dcmp(q[tail-1].ang-q[tail].ang)==0){ tail--; if(Onleft(L[i].P,q[tail]))q[tail]=L[i]; } if(head<tail)P[tail-1]=LineMeet(q[tail-1],q[tail]); } while(head<tail&&!Onleft(P[tail-1],q[head]))tail--; return tail-head>1; } void Solve(){ int leftt=0,rightt=N-1,mid; while(leftt+1<rightt){ mid=(leftt+rightt)>>1; Modify(mid); if(check()) leftt=mid; else rightt=mid; } Modify(rightt); if(check()) cout<<rightt<<endl; else cout<<leftt<<endl; }}int main(){ //freopen("input.in","r",stdin); using namespace solution; while(scanf("%d",&N)!=EOF){ Prepare(); Solve(); } return 0;}
UVaLive4992:Jungle Outpost
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。