hdu 5033 单调栈的维护

题意：在一条直线上，给定一系列摩天大楼的位置和高度，每次询问自己所在位置所能看到的视野。

分析：关键在于楼与楼之间，到当前将要加进去的楼或人的位置，若之前的楼比将要进来的楼矮，可以忽略，因为对后面有影响的，只有比将要加进去的楼高的楼，最终构成递减序列，这是最基本的要求。对于从左边开始连续的三幢楼a,b,c,其中以c为基点，若a与c构成的仰角>=b与c构成的仰角，则b对于后面更远的楼层或人的位置将不起作用，故b可以忽略此时。

#include<iostream>#include<cstdio>#include<algorithm>#include<cmath>#include<cstring>using namespace std;struct node{    double x,h;    int id;    bool operator<(const node &p)const    {        return x<p.x;    }};const int MAX=100005;const double PI=acos(-1.0),EPS=1e-8;node p[MAX<<1],st[MAX<<1];double ans[MAX];int n,q;bool is_out(node &a,node &b,node &c){    return (a.h-c.h)*(c.x-b.x)>=(b.h-c.h)*(c.x-a.x);//以c为基点，比较a,b两幢楼的倾斜角的大小，若b较小或等于，说明b被覆盖，出栈！}double getAngle(node & a,node &b){    double an=90.0-180.0*atan((double)a.h/(double)(b.x-a.x))/PI;//返回人所在位置为b时，与左边楼a的连线与竖直方向构成的倾斜角大小！    //printf("%d %d %lf\n",a.h,b.x-a.x,an);    return an;}/*实时维护单调递减的栈，即楼高要递减，且一定包含到目前为止将要入栈的楼，不满足的出栈，同时若栈顶所在楼被覆盖，则不断出栈！直到栈顶所在楼是目前为止视线最窄的。*/void getAns(){    int top=0;    for(int i=0; i<n+q; i++)    {        //printf("current %d %d\n",p[i].x,p[i].h);        while(top>=1&&st[top-1].h<=p[i].h)        {            top--;            //printf("%d %d 高度低\n",st[top-1].x,st[top-1].h);        }        while(top>=2&&is_out(st[top-2],st[top-1],p[i]))        {            top--;            //printf("%d %d 倾斜角低\n",st[top-1].x,st[top-1].h);        }        if(p[i].h<EPS)            ans[p[i].id]+=getAngle(st[top-1],p[i]);        else            st[top++]=p[i];        //printf("end\n");    }}int main(){    int t,i,cas=1;    //freopen("in.txt","r",stdin);    scanf("%d",&t);    while(t--)    {        scanf("%d",&n);        for(i=0; i<n; i++)            scanf("%lf%lf",&p[i].x,&p[i].h);        scanf("%d",&q);        for(i=0; i<q; i++)        {            scanf("%lf",&p[i+n].x);            p[i+n].h=0;            p[i+n].id=i;//id保存的是人的位置。        }        sort(p,p+n+q);        memset(ans,0,sizeof(ans));        getAns();//位置x递增求左边视线角。        reverse(p,p+n+q);        for(i=0; i<n+q; i++)            p[i].x=10000005.0-p[i].x;//原来位置右边x最大，现在变最小，求右边视线角！        getAns();        printf("Case #%d:\n",cas++);        for(i=0; i<q; i++)            printf("%.10lf\n",ans[i]);    }    return 0;}