标签:Building HDU 斜率 int top 5033 atan2 &&
pro:现在在X轴上有N个摩天大楼,以及Q个人,人和大楼的坐标各不相同,保证每个人左边和右边都有楼,问每个人能看到天空的角度大小。
sol:不难想到就是维护凸包,此题就是让你模拟斜率优化,此处没有斜率来做,用几何写的。。。。
#include<bits/stdc++.h> #define rep(i,a,b) for(int i=a;i<=b;i++) using namespace std; const int maxn=200010; const double pi=acos(-1.0); struct in{ double x,h; int id; }s[maxn]; struct point{ double x,y; point(){} point(double xx,double yy):x(xx),y(yy){} }; double det(point a,point b){ return a.x*b.y-a.y*b.x;} double dot(point a,point b){ return a.x*b.x+a.y*b.y;} bool cmp(in w,in v){ return w.x<v.x;} double ans[maxn]; int q[maxn],top; void solve(int N) { sort(s+1,s+N+1,cmp); top=0; rep(i,1,N){ if(s[i].id){ while(top>1&&atan2(s[q[top]].h-s[q[top-1]].h,s[q[top]].x-s[q[top-1]].x) <atan2(-s[q[top]].h,s[i].x-s[q[top]].x)) top--; point T=point(s[q[top]].h,s[q[top]].x-s[i].x); ans[s[i].id]+=asin(fabs(det(point(0,-1),T))/sqrt(dot(T,T))); } else { while(top&&s[q[top]].h<=s[i].h) top--; while(top>1&&atan2(s[q[top]].h-s[q[top-1]].h,s[q[top]].x-s[q[top-1]].x) <atan2(s[i].h-s[q[top]].h,s[i].x-s[q[top]].x)) top--; q[++top]=i; } } } int main() { int T,N,Q,C=0; scanf("%d",&T); while(T--){ scanf("%d",&N); rep(i,1,N){ scanf("%lf%lf",&s[i].x,&s[i].h); s[i].id=0; } scanf("%d",&Q); rep(i,1,Q){ scanf("%lf",&s[N+i].x); s[N+i].id=i,ans[i]=0; s[N+i].h=0; } solve(N+Q); rep(i,1,N+Q) s[i].x=-s[i].x; solve(N+Q); printf("Case #%d:\n",++C); rep(i,1,Q) printf("%.10lf\n",180-180/pi*ans[i]); } return 0; }
标签:Building,HDU,斜率,int,top,5033,atan2,&& 来源: https://www.cnblogs.com/hua-dong/p/10686296.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。