BZOJ 2957: 楼房重建

来源：互联网发布：微信秒抢红包软件苹果编辑：程序博客网时间：2024/06/08 06:58

这道题还不错….我太蠢了想不来QAQ
膜了一下线段树做法之后发现还有简单粗暴的分块算法好劲啊
线段树做法比较巧妙推荐去看看黄学长的有个细节要注意取右边的贡献时要用区间贡献减去左边的贡献因为右边的贡献在算的时候并没有被左边约束
ORZ HZWER
然后是分块块内维护一个单调上升序列跟普通的不一样意会一下就懂了
每次就直接块内重新构造序列每个块依次二分求值 Po姐那里会告诉你怎样分块更快。。
ORZ PoPoQQQ
这道题线段树的速度快一点。。
线段树

#include<bits/stdc++.h>using namespace std;const int N=400005;inline int read(){    char ch=getchar(); int x=0,f=1;    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}    while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}    return x*f;}int c[N];double mx[N];int query(double k,int x,int l,int r){    int mid=(l+r)>>1,lc=x<<1,rc=lc|1;    if(l==r)return k<mx[x];    if(mx[lc]<=k)return query(k,rc,mid+1,r);    return c[x]-c[lc]+query(k,lc,l,mid);}void change(double k,int u,int x,int l,int r){    int mid=(l+r)>>1,lc=x<<1,rc=lc|1;    if(l==r){c[x]=1,mx[x]=k; return;}    if(u<=mid)change(k,u,lc,l,mid);    else change(k,u,rc,mid+1,r);    mx[x]=max(mx[lc],mx[rc]);    c[x]=c[lc]+query(mx[lc],rc,mid+1,r);}int main(){    int n=read(),m=read(); double x,y;    while(m--)    {        x=read(),y=read();        change(y/x,(int)x,1,1,n);        printf("%d\n",c[1]);    }    return 0;}

分块

#include<bits/stdc++.h>using namespace std;const int N=100010;const double eps=1e-10;inline int read(){    char ch=getchar(); int x=0,f=1;    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}    while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}    return x*f;}int bl[N]; double a[N];struct node{    double s[1010];    int l,r,tot;    void bulid()    {        int i; double u=0;        tot=0;        for(i=l;i<=r;i++)            if(a[i]>u+eps)                u=a[i],s[++tot]=a[i];    }    int query(double &x)    {        int r=tot+1-(upper_bound(s+1,s+1+tot,x+eps)-s);        if(r)x=s[tot];        return r;    }}e[120];int main(){    int n=read(),m=read(),b,i,j,x,y;    b=sqrt(n*log(n)/log(2)/2.0);    for(i=1;(i-1)*b+1<=n;i++)        e[i].l=(i-1)*b+1,e[i].r=min(i*b,n);    for(i=1;i<=n;i++)        bl[i]=(i-1)/b+1;    while(m--)    {        x=read(),y=read();        a[x]=(double)y/x; e[bl[x]].bulid();        double u=0; int ans=0;        for(j=1;e[j].l;j++)ans+=e[j].query(u);        printf("%d\n",ans);    }    return 0;}

0 0