BZOJ2957: 楼房重建

来源:互联网 发布:《张国荣的时光》知乎 编辑:程序博客网 时间:2024/06/08 17:13

楼房i被看见的条件是他和起点连线的斜率是1~i里最大的

可以用线段树做
维护每个区间的ans,max
合并区间的时候,
左区间的ans是可以直接加的,
对于右区间,当左区间的max>右区间max时,整个区间被挡住,没有贡献
当左区间的max<右区间max时,需要递归下右区间找
递归过程中,带一个这个区间前方的斜率max,
当max<左子区间的max时,这个带下来的max对右子区间的贡献无影响,通过整个区间的ans-左子区间的ans算粗右子区间的答案直接加上就行,对于左子区间再递归下去
当max>左子区间的max时,左子区间无贡献,递归下右子区间
这个做法..我总觉得是log2

然后….说了这么多…其实我没打线段树…..(憋打我)

我写了分块
分块的话,每个块内维护一个斜率单调增的子序列,修改一个值时,这个块内暴力重构,因为修改一个值只会对右方的贡献产生影响,向右顺着块扫一遍,记录之前的最大值,在每个块内二分找有多少个值比他大
code:

#include<set>#include<map>#include<deque>#include<queue>#include<stack>#include<cmath>#include<ctime>#include<bitset>#include<string>#include<vector>#include<cstdio>#include<cstdlib>#include<cstring>#include<climits>#include<complex>#include<iostream>#include<algorithm>#define ll long longusing namespace std;const int maxn = 110000;const int maxk = 110;const double eps=1e-10;int id[maxn],L[maxk];int len[maxk],N;double a[maxn],q[maxk][1200];int find_(int k,double x){    int l=1,r=len[k];    while(l<=r)    {        int mid=(l+r)>>1;        if(q[k][mid]-x>eps) r=mid-1;        else l=mid+1;    }    return r+1;}void rebuild_(int k){    int llen=0;    for(int i=L[k];i<L[k+1];i++)if(a[i]-q[k][llen]>eps)             q[k][++llen]=a[i];    len[k]=llen;}int n,m;int main(){    scanf("%d%d",&n,&m); N=(3.0*sqrt(n));    for(int i=1;i<=n;i++) id[i]=(i-1)/N+1;    for(int i=1;i<=id[n];i++) L[i]=(i-1)*N+1;    L[id[n]+1]=n+1;    while(m--)    {        int x,y; scanf("%d%d",&x,&y);        a[x]=(double)y/x;        rebuild_(id[x]);        double las=0.0; int ret=0;        for(int i=1;i<=id[n];i++)if(len[i])        {            int k=find_(i,las);            if(k<=len[i]) ret+=len[i]-k+1,las=q[i][len[i]];        }        printf("%d\n",ret);    }    return 0;}
0 0