cf613 E Delete a Segment

题意

有$(n<=1e5)$个区间，要求删掉一个区间，使剩下的区间不覆盖的个数。

思路

2中做法

线段树

首先给出的区间是乱序的，对其排个序，然后我们可以一次枚举删掉这个区间，然后求剩下的区间所剩的个数。可以用线段树来进行更新查询，可以类似求区间连续01个数那种方式，但我们由于这里是区间，所以把点化为区间，把区间扩大2倍$[2\ast l,2\ast r]$，然后多加入这个$[2\ast l+1,2\ast r+1]$来进行离散化，最后就是区间线段树的操作，只不过pushup的时候我们有细节，如果这个区间全部被覆盖，那么就不用从下面得到答案，直接更新为1就ok了

code

#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;
const int man = 4e5+10;
#define IOS ios::sync_with_stdio(0)
template <typename T>
inline T read(){T sum=0,fl=1;int ch=getchar();
for(;!isdigit(ch);ch=getchar())if(ch=='-')fl=-1;
for(;isdigit(ch);ch=getchar())sum=sum*10+ch-'0';
return sum*fl;}
template <typename T>
inline void write(T x) {static int sta[35];int top=0;
do{sta[top++]= x % 10, x /= 10;}while(x);
while (top) putchar(sta[--top] + 48);}
template<typename T>T gcd(T a, T b) {return b==0?a:gcd(b, a%b);}
template<typename T>T exgcd(T a,T b,T &g,T &x,T &y){if(!b){g = a,x = 1,y = 0;}
else {exgcd(b,a%b,g,y,x);y -= x*(a/b);}}
#ifndef ONLINE_JUDGE
#define debug(fmt, ...) {printf("debug ");printf(fmt,##__VA_ARGS__);puts("");}
#else
#define debug(fmt, ...)
#endif
typedef long long ll;
const ll mod = 1e9+7;
int a[man*2];

struct node{
    int cnt;
    bool lc,rc;//这段区间的左右区间是否被覆盖
    node(){
        cnt = 0;
        lc = rc = 0;
    }
}tr[man*30];
int flag[man*30];

void build(int l,int r,int rt){
    flag[rt] = 0;
    tr[rt].cnt = 0;
    tr[rt].lc = tr[rt].rc = 0;
    if(l==r){
        //把叶子的下面给清掉
        tr[rt<<1].cnt = tr[rt<<1|1].cnt = 0;
        tr[rt<<1].lc = tr[rt<<1].rc = 0;
        tr[rt<<1|1].lc = tr[rt<<1|1].rc = 0;
        return;
    }
    int m = l + r >> 1;
    build(l,m,rt<<1);
    build(m+1,r,rt<<1|1);
}

void pushup(int rt){
    if(flag[rt]){//如果这整个区间都覆盖，则就为1了，
        tr[rt].cnt = 1;
        tr[rt].lc = tr[rt].rc = 1;
    }else{
        tr[rt].cnt = tr[rt<<1].cnt  + tr[rt<<1|1].cnt;
        if(tr[rt<<1].rc&&tr[rt<<1|1].lc)tr[rt].cnt--;
        tr[rt].lc = tr[rt<<1].lc;
        tr[rt].rc = tr[rt<<1|1].rc;
    }
}

void update(int l,int r,int L,int R,int rt,int v){
    if(L<=l&&R>=r){
        flag[rt]+= v;//打标机并更新。
        pushup(rt);
        return;
    }
    int m = l + r >>1;  
    if(L<=m)update(l,m,L,R,rt<<1,v);
    if(R>m)update(m+1,r,L,R,rt<<1|1,v);
    pushup(rt);
}

pair<int,int>b[man];
int main() {
    #ifndef ONLINE_JUDGE
        //freopen("in.txt", "r", stdin);
        //freopen("out.txt","w",stdout);
    #endif
    int t;
    cin >> t;
    while(t--){
        int n;
        cin >> n;
        int cnt = 0;
        for(int i = 1;i <= n;i++){
            int l,r;
            cin >> l >> r;
            l = l*2;r = 2*r;
            b[i].first = l;
            b[i].second = r;
            a[cnt++] = l;
            a[cnt++] = l+1;
            a[cnt++] = r;
            a[cnt++] = r+1;
        }
        sort(a,a+cnt);
        cnt = unique(a,a+cnt) - a;
        build(1,cnt,1);
        int ans = -1;
        for(int i = 1;i <= n;i++){
            int l = lower_bound(a,a+cnt,b[i].first) - a + 1;
            int r = lower_bound(a,a+cnt,b[i].second) - a + 1;
            update(1,cnt,l,r,1,1);
        }
        for(int i = 1;i <= n;i++){
            int l = lower_bound(a,a+cnt,b[i].first) - a + 1;
            int r = lower_bound(a,a+cnt,b[i].second) - a + 1;
            update(1,cnt,l,r,1,-1);
            ans = max(ans,tr[1].cnt);
            update(1,cnt,l,r,1,1);
        }
        cout << ans << endl;
    }   
    return 0;
}

单调栈

首先考虑我们删掉中间的一个 线段，左边的个数我们可以通过前缀num得到，num[i]就是前i个区间有多少个不相邻的区间，然后可以用max_r来得到,如果当前的L大于max_r[i-1]，那么num就会增加一个，依次累加就ok，那么右边的区间又如何得到，因为我们的l从后往前，依次递减，然后我们只需要最小的L，与最小的L的区间相互重叠的区间我们是用不到的，因为可以看成一个区间，所以只有这个最小的L有用，然后右边的一部分我们可以用单调栈维护L（后往前）。
整个过程就是从后往前得到答案，然后维护单调栈里的元素，每次的答案就是num[i-1] + cnt
cnt就是单调栈里大于max_r[i-1]的个数，可以二分得到

code

#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;
const int man = 2e5+10;
#define IOS ios::sync_with_stdio(0)
template <typename T>
inline T read(){T sum=0,fl=1;int ch=getchar();
for(;!isdigit(ch);ch=getchar())if(ch=='-')fl=-1;
for(;isdigit(ch);ch=getchar())sum=sum*10+ch-'0';
return sum*fl;}
template <typename T>
inline void write(T x) {static int sta[35];int top=0;
do{sta[top++]= x % 10, x /= 10;}while(x);
while (top) putchar(sta[--top] + 48);}
template<typename T>T gcd(T a, T b) {return b==0?a:gcd(b, a%b);}
template<typename T>T exgcd(T a,T b,T &g,T &x,T &y){if(!b){g = a,x = 1,y = 0;}
else {exgcd(b,a%b,g,y,x);y -= x*(a/b);}}
#ifndef ONLINE_JUDGE
#define debug(fmt, ...) {printf("debug ");printf(fmt,##__VA_ARGS__);puts("");}
#else
#define debug(fmt, ...)
#endif
typedef long long ll;
const ll mod = 1e9+7;
pair<int,int>sp[man];
int num[man],sta[man],max_r[man];
int top;

int bin(int x){
    int l = 1,r = top;
    int ans = 0;
    while(l<=r){
        int mid = l + r >> 1;
        if(sta[mid]>=x){
            l = mid + 1;
            ans = mid;
        }else r = mid - 1;
    }
    if(sta[ans]==x)ans--;
    return ans;
}


int main() {
    #ifndef ONLINE_JUDGE
        //freopen("in.txt", "r", stdin);
        //freopen("out.txt","w",stdout);
    #endif
    int t;
    scanf("%d",&t);
    while(t--){
        int n;
        scanf("%d",&n);
        max_r[0] = -INT_MAX;
        num[0] = 0;top= 0;
        for(int i = 1;i <= n;++i){
            scanf("%d%d",&sp[i].first,&sp[i].second);
        }
        sort(sp+1,sp+1+n);
        for(int i = 1;i <= n;++i){
            max_r[i] = max(sp[i].second,max_r[i-1]);
            if(sp[i].first>max_r[i-1])num[i] = 1;
            num[i] += num[i-1];
        }
        int ans = 0;
        for(int i = n;i >= 1;i--){
            ans = max(ans,num[i-1]+bin(max_r[i-1]));
            while(top&&sp[i].second>=sta[top])top--;
            sta[++top] = sp[i].first;
            num[i] = 0;
        }
        printf("%d\n",ans);
    }
    return 0;
}