A - Circuits(线段树维护区间最大值)

Problem

求两条平行线穿过的最大矩形。

Solution

• 贪心走不通之后场上想到了枚举，考虑到了定一条之后枚举剩下的，然后剪纸。实际上从剩下的里面找最大就是个区间查询最值问题，线段树可维护。
• 不过还有一些细节要处理：首先区间很长，要对矩形进行离散化；每次确定一条的时候不可能把穿过的所有矩形的对应区间进行修改（这样做可能比暴力还差），实际上这样产生的来回修改毫无意义，顺着这条优化思路，能不能想个办法保证每个矩形只对线段树做两次修改操作（离散化产生两个点）？选择一个矩形下边界的时候将其在区间删掉，在选择到其上边界之前矩形应当一直是删掉的状态（第一条平行线穿过）。具体来说，先将优先队列中上边界小于当前位置的矩形恢复，单点查询第一条线穿过的矩形数，以当前位置为下边界的矩形将其在区间删掉，并将其加进优先队列，维护优先队列的元素数。当前答案为恢复后优先队列元素数+单点查询+删除后的区间最值。
• 其它：从线段树原理以及在此情形下的应用出发，lazy_tag不能省；注意数组越界！

Code

#include<cstdio>
#include<queue>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<vector>
using namespace std;
const int maxn=1e5+5;
int n,m,tree[maxn*8],mx,cnt,tag[maxn*8];
struct posi{
    int x,id;
}pos[maxn*2];
struct rectangle{
    int yl,yr;
    bool operator<(const rectangle &a)const{
        return yr>a.yr;
    }
}rec[maxn];
bool cmp(const posi &a,const posi &b)
{
    return a.x<b.x;
}
vector<int>points[maxn*2];
void pushup(int x,int l,int r)
{
    tree[x]=max(tree[x<<1],tree[x<<1|1]);
}
void pushdown(int x)
{
	int ll=x<<1,rr=x<<1|1;

	tag[ll]+=tag[x];
	tree[ll]+=tag[x];
	tag[rr]+=tag[x];
	tree[rr]+=tag[x];
	tag[x]=0;

}
void update(int x,int l,int r,int ll,int rr,int k)
{
    if(l>rr||r<ll) return ;
    if(l>=ll&&r<=rr)
    {
        tree[x]+=k;
        tag[x]+=k;
        return ;
    }
    pushdown(x);
    int mid=(l+r)>>1;
    update(x<<1,l,mid,ll,rr,k);
    update(x<<1|1,mid+1,r,ll,rr,k);
    pushup(x,l,r);
}
int query(int x,int l,int r,int ll,int rr)
{
    if(l>rr||r<ll) return 0;
    if(l>=ll&&r<=rr)
    return tree[x];
    pushdown(x);
    int mid=(l+r)>>1;
    return max(query(x<<1,l,mid,ll,rr),query(x<<1|1,mid+1,r,ll,rr));
}
int main()
{
    pos[0].x=-1e8-5;
    scanf("%d",&n);
    for(int i=1,x1,x2,y1,y2;i<=n;i++)
    {
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        rec[i].yl=-1e8-5;
        pos[(i-1)*2+1].x=y1;
        pos[i*2].x=y2;
        pos[(i-1)*2+1].id=pos[i*2].id=i;
    }
    sort(pos+1,pos+1+2*n,cmp);
    for(int i=1;i<=2*n;i++)
    {
        int num=pos[i].id;
        if(pos[i].x!=pos[i-1].x)
            m++;
        points[m].push_back(pos[i].id);
        if(rec[num].yl==-1e8-5)
            rec[num].yl=m;
        else
            rec[num].yr=m;
    }
    for(int i=1;i<=n;i++)
        update(1,1,m,rec[i].yl,rec[i].yr,1);
    priority_queue<rectangle>q;
    cnt=0;
    for(int i=1;i<=m;i++)
    {
        while(!q.empty())
        {

            rectangle rk=q.top();
            if(rk.yr>=i)
                break;
            cnt--;
            q.pop();
            update(1,1,m,rk.yl,rk.yr,1);
        }
        int ans=query(1,1,m,i,i)+cnt;
        for(int j=0;j<points[i].size();j++)
        {
            if(i!=rec[points[i][j]].yl)
                continue;
            update(1,1,m,rec[points[i][j]].yl,rec[points[i][j]].yr,-1);
            cnt++;
            q.push({rec[points[i][j]].yl,rec[points[i][j]].yr});
        }
        ans+=query(1,1,m,1,m);
        mx=max(mx,ans);
    }
    cout<<mx<<endl;
    return 0;
}