P5445 [APIO2019]路灯

传送门·

对于询问 $(a,b)$ ，感觉一维很不好维护，考虑把询问看成平面上的一个点，坐标为 $(a,b)$

每个坐标 $(x,y)$ 的值表示到当前 $x$ 和 $y$ 联通的时间和

考虑一个修改的贡献，它其实就是把左边一段区间 $[l,x]$ 和右边一段区间 $[x+1,r]$ 联通或断开

放到平面上发现其实就是横坐标在 $[l,x]$ ，纵坐标在 $[x+1,r]$ 的矩形里修改，那么矩形左下角为 $[l,x+1]$，右上角为 $[x,r]$

如果每个时间点都把相应矩形 $+1$ 的话显然是不可行的，考虑起点和终点的时间差

如果当前操作是联通，设当前时间为 $t$，则把相应矩形 $-t$，如果是断开则把矩形加 $t$

这样我们询问时直接矩形单点查值即可，但是要注意，如果查询时当前区间仍然联通，那么对应矩形还没加 $t$ ，我们查完值以后答案还要再加 $t$

矩形加具体就是差分成四个前缀修改，查询的时候直接查左下所有位置的和即可

如果把左下角为 $(xa,ya)$，右上角 $(xb,yb)$ 的闭区间都加一个 $v$，那么其实就是把 $(0,0),(xa,ya)$ 加一个 $v$，$(0,0),(xa,yb+1)$ 减 $v$，$(0,0),(xb+1,ya)$ 减 $v$，$(0,0),(xb+1,yb+1)$ 加 $v$

二维平面加一维时间轴， $CDQ$ 分治即可

然后就是具体的维护了，我维护亮灯区间用的是闭区间，好像比较多细节...

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<set>
using namespace std;
typedef long long ll;
inline int read()
{
    int x=,f=; char ch=getchar();
    while(ch<''||ch>'') { if(ch=='-') f=-; ch=getchar(); }
    while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
    return x*f;
}
const int N=6e5+;
char s[N];
int n,Q,tot;
struct dat{//区间
    int l,r;
    dat (int a=,int b=) { l=a,r=b; }
    inline bool operator < (const dat &tmp) const {
        return l<tmp.l;
    }
};
set <dat> S;//维护开灯的区间
set <dat>::iterator ita,itb;
struct Poi{//CDQ分治的点
    int x,y,id,v;//id=0表示是修改，v是修改的值
    Poi (int a=,int b=,int c=,int d=) { x=a,y=b,id=c,v=d; }
}d[N<<],tmp[N<<];
int t[N];
inline void add(int x,int v) { while(x<=n+) t[x]+=v,x+=x&-x; }
inline int ask(int x) { int res=; while(x) res+=t[x],x-=x&-x; return res; }
inline bool fc(Poi &a,Poi &b) { if(a.x!=b.x) return a.x<b.x; return a.y!=b.y ? a.y<b.y : a.id<b.id; }
int ans[N];
void solve(int l,int r)
{
    if(l==r) return; int mid=l+r>>;
    solve(l,mid); solve(mid+,r);
    int i=l,j=mid+,p=l-;
    while(i<=mid&&j<=r)
    {
        if(fc(d[i],d[j]))
        {
            if(!d[i].id) add(d[i].y,d[i].v);
            tmp[++p]=d[i++]; continue;
        }
        if(d[j].id) ans[d[j].id]+=ask(d[j].y);
        tmp[++p]=d[j++];
    }
    while(i<=mid) { if(!d[i].id) add(d[i].y,d[i].v); tmp[++p]=d[i++];  }
    while(j<=r) { if(d[j].id) ans[d[j].id]+=ask(d[j].y); tmp[++p]=d[j++]; }
    for(int k=l;k<=mid;k++) if(!d[k].id) add(d[k].y,-d[k].v);//撤销
    for(int k=l;k<=r;k++) d[k]=tmp[k];
}
int main()
{
    n=read(),Q=read();
    scanf("%s",s+);
    S.insert(dat(,)); S.insert(dat(n+,n+));//初始虚区间防越界
    for(int i=;i<=n;i++)
    {
        if(s[i]=='') continue;
        int R=i; while(s[R+]=='') R++;
        S.insert(dat(i,R)); i=R;//初始情况
    }
    char ch[]; int a,b,cnt=;
    for(int i=;i<=Q;i++)
    {
        scanf("%s",ch); a=read();
        if(ch[]=='q')
        {
            b=read();
            d[++tot]=Poi(a,b,++cnt,);//询问点
            ita=S.lower_bound(dat(a+,)); ita--;
            itb=S.lower_bound(dat(b,)); itb--;
            if( ita==itb && (*ita).r>=b- ) ans[cnt]+=i;//判断区间是否还联通
            continue;
        }
        if(s[a]=='')//联通
        {
            ita=S.lower_bound(dat(a+,)); ita--; itb=ita; ita++;//找到左右区间
            dat A=(*ita),B=(*itb); int l=a,r=a;
            if(B.r==l-&&B.l>=) l=B.l,S.erase(B);//注意判断,注意边界
            if(A.l==r+&&A.r<=n) r=A.r,S.erase(A);
            S.insert(dat(l,r));
            d[++tot]=Poi(l,a+,,-i); d[++tot]=Poi(a+,r+,,-i);
            d[++tot]=Poi(l,r+,,i); d[++tot]=Poi(a+,a,,i);
            //联通时矩形-i
            s[a]=''; continue;
        }
        //断开
        ita=S.lower_bound(dat(a+,)); ita--; dat A=(*ita);
        S.erase(A); if(A.l<a) S.insert(dat(A.l,a-)); if(A.r>a) S.insert(dat(a+,A.r));//注意边界
        d[++tot]=Poi(A.l,a+,,i); d[++tot]=Poi(a+,A.r+,,i);
        d[++tot]=Poi(A.l,A.r+,,-i); d[++tot]=Poi(a+,a,,-i);
        //断开时矩形+i
        s[a]='';
    }
    solve(,tot);
    for(int i=;i<=cnt;i++) printf("%d\n",ans[i]);
    return ;
}