HDU5033 building 单调栈+计算几何

正解:单调栈

解题报告:

哇生气辽QAQ本来打了半天feel good都快调出来了然后说换题了QAQ(所以可能那题的代码会过一阵子再放上来了QAQ

不过还是大爆手速打了一通拿到首杀了嘻嘻

美滋滋辽

然后放个传送门先qwq

有一点点解题思想和oil是有点儿像的,就是先把询问和修改放一个结构体里排序,这样当到达某个询问的时候就可以直接解决不用进行其他修改了qwq

然后就是做两次单调栈一次算出往左能看到的角度一次算出往右能看到的角度

然后就欧克辽?

ummm好像是没什么好说的了只要推出了角度计算公式然后按斜率维护栈,然后就差不多了?

顺便是不是可以总结下?就是这种计算几何询问+修改的经常会有放在一个数据结构里然后排序的操作?这样通常都会方便一些呢qwq

over?

代码放下我就做后面的题目辣qwq

啊还有就是这题我不是港说计算几何嘛?其实并没有用到计算几何的方法的说感觉,只是在这个范畴内所以就随手cue了一下qwq

#include <bits/stdc++.h>
using namespace std;
#define rp(i,x,y) for(register int i=x;i<=y;++i)

const double PI=acos(-1.0) ;
+,inf=;
struct node
{
    int x,h;
    bool operator<(const node &R)const{return x<R.x;}
}bd[N<<],stk[N<<];
double ans[N];
int n,q;

inline int read()
{
    ;;
    '))ch=getchar();
    ;
    )+(x<<)+(ch^'),ch=getchar();
    return y?x:-x;
}
inline bool judge(node a,node b,node c)
{
    )c.h=;
    return (long long)(a.h-c.h)*(b.x-c.x)<=(long long)(b.h-c.h)*(a.x-c.x);
}
inline double xl(node a,node b){return atan(1.0*(b.x-a.x)/a.h);}
inline void solve()
{
    ;
    rp(i,,n+q)
    {
        )
        {
             && judge(stk[top-],stk[top-],bd[i]))--top;
            ans[-bd[i].h]+=xl(stk[top-],bd[i]);
            continue;
        }
        ].h<=bd[i].h) || (top>= && judge(stk[top-],stk[top-],bd[i])))--top;
        stk[top++]=bd[i];
    }
}

int main()
{
    int T=read();
    rp(j,,T)
    {
        n=read();rp(i,,n)bd[i].x=read(),bd[i].h=read();
        q=read();rp(i,,q)bd[i+n].x=read(),bd[i+n].h=-i;
        memset(ans,,,bd+n+q+);solve();
        reverse(bd+,bd+n+q+);rp(i,,n+q)bd[i].x=inf-bd[i].x;solve() ;
        printf(,q)printf(/PI) ;
    }
    ;
}
//话说插入代码中编程语言c#和c++看起来有什么区别嘛qwq我有时候选c#有时候选c++感觉界面都一样鸭qwq

这儿!是!代码!QwQ