BZOJ4422[Cerc2015]Cow Confinement（扫描线+线段树）

很容易发现一个O(n²)DP，f[i][j]=f[i][j+1]+f[i+1][j]-f[i+1][j+1]。然后由于有栅栏，一些位置没办法走，然后就可以用类似差分的方法，f[i]表示当前行f[i+1]无法到达的花朵，然后对于每个点找到其下方第一个栅栏。分情况讨论，需要支持单点修改（出现花），区间标记覆盖（出现一个栅栏），以及区间归零（栅栏走了），当然还要区间查询（牛出现了），然后还要支持第一个障碍物的查询，不过这些用一个操作较多的线段树就可以跑过了，复杂度一个log可以“轻松”跑过。

#include<cstdio>
#include<algorithm>
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
using namespace std;
const int N=2e6+,M=1e6;
struct fence{int xl,xr,y,i;bool tp;}a[N<<];
bool operator<(fence a,fence b){return a.y!=b.y?a.y>b.y:a.xl<b.xl;}
struct flower{int x,y;}b[N];
bool operator<(flower a,flower b){return a.y>b.y;}
struct cow{int x,y,i;}c[N];
bool operator<(cow a,cow b){return a.y>b.y;}
struct seg{int num;bool cov,cut;}tr[N<<];
int n,m,val[N],ans[N];
void pushup(int rt)
{
    tr[rt].num=tr[rt<<].num+tr[rt<<|].num;
    tr[rt].cut=tr[rt<<].cut|tr[rt<<|].cut;
}
void modify(int rt){tr[rt].cov=,tr[rt].num=;}
void pushdown(int rt){if(tr[rt].cov)modify(rt<<),modify(rt<<|),tr[rt].cov=;}
void add(int k,int v,int l,int r,int rt)
{
    tr[rt].num+=v;
    if(l==r)return;
    pushdown(rt);
    int mid=l+r>>;
    if(k<=mid)add(k,v,lson);else add(k,v,rson);
    pushup(rt);
}
void cover(int L,int R,int l,int r,int rt)
{
    if(L<=l&&r<=R){modify(rt);return;}
    pushdown(rt);
    int mid=l+r>>;
    if(L<=mid)cover(L,R,lson);
    if(R>mid)cover(L,R,rson);
    pushup(rt);
}
int query(int L,int R,int l,int r,int rt)
{
    if(L<=l&&r<=R)return tr[rt].num;
    pushdown(rt);
    int mid=l+r>>,ret=;
    if(L<=mid)ret+=query(L,R,lson);
    if(R>mid)ret+=query(L,R,rson);
    return ret;
}
void update(int k,int l,int r,int rt)
{
    if(l==r){tr[rt].cut^=;return;}
    pushdown(rt);
    int mid=l+r>>;
    if(k<=mid)update(k,lson);else update(k,rson);
    pushup(rt);
}
int getnxt(int L,int l,int r,int rt)
{
    if(l>=L)
    {
        if(tr[rt].cut)
        {
            while(l!=r)
            if(tr[rt<<].cut)rt<<=,r=l+r>>;
            else rt=rt<<|,l=(l+r>>)+;
            return l;
        }
        return ;
    }
    int tmp,mid=l+r>>;
    pushdown(rt);
    if(L<=mid&&(tmp=getnxt(L,lson)))return tmp;
    return getnxt(L,rson);
}
int main()
{
    int f;scanf("%d",&f);
    for(int i=,x1,y1,x2,y2;i<f;i++)
    {
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        a[i<<]=(fence){x1,x2,y1-,i,};
        a[i<<|]=(fence){x1,x2,y2,i,};
    }
    sort(a,a+(f<<));
    scanf("%d",&m);
    for(int i=;i<=m;i++)scanf("%d%d",&b[i].x,&b[i].y);
    sort(b+,b+m+);
    scanf("%d",&n);
    for(int i=;i<=n;i++)scanf("%d%d",&c[i].x,&c[i].y),c[i].i=i;
    sort(c+,c+n+);
    f=;
    update(M,,M,);
    for(int i=M,p=,q=;i;i--)
    {
        int sum,cut;
        while(a[f].y==i)
        {
            if(!a[f].tp)
            {
                cover(a[f].xl,a[f].xr,,M,);
                if(a[f].xl!=)add(a[f].xl-,-val[a[f].i],,M,);
                if(a[f].xl!=)update(a[f].xl-,,M,);
                if(a[f].xr!=M)update(a[f].xr,,M,);
            }
            else{
                cut=getnxt(a[f].xr,,M,),sum=query(a[f].xl,a[f].xr,,M,);
                val[a[f].i]=query(a[f].xr+,cut,,M,);
                cover(a[f].xl,a[f].xr,,M,);
                if(a[f].xl>)add(a[f].xl-,sum+val[a[f].i],,M,);
                if(a[f].xl!=)update(a[f].xl-,,M,);
                if(a[f].xr!=M)update(a[f].xr,,M,);
            }
            f++;
        }
        while(b[p].y==i)add(b[p].x,,,M,),++p;
        while(c[q].y==i)cut=getnxt(c[q].x,,M,),ans[c[q].i]=query(c[q].x,cut,,M,),++q;
    }
    for(int i=;i<=n;i++)printf("%d\n",ans[i]);
}