【BZOJ1018】堵塞的交通traffic（线段树，网格图，连通性）

题意：一个2行C列的矩形网格图，网格上的每个点代表一个城市，相邻的城市之间有一条道路

一开始每条道路都是堵塞的，堵塞即为不可经过。经过一些操作后，可能某些道路通畅了，也可能某些道路堵塞了

多次询问，询问两个城市是否联通

C,q<=1e5

思路：From https://www.cnblogs.com/MashiroSky/p/5973686.html

#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<map>
#include<set>
#include<queue>
#include<vector>
using namespace std;
typedef long long ll;
typedef unsigned int uint;
typedef unsigned long long ull;
typedef pair<int,int> PII;
typedef vector<int> VI;
#define fi first
#define se second 
#define MP make_pair
#define N   110000
#define M   41
#define eps 1e-8
#define pi acos(-1)

struct node
{
    int U,D,l,r,u,d,p,q;
}t[N<<2];
int n; 
char s[10];

int read()
{ 
   int v=0,f=1;
   char c=getchar();
   while(c<48||57<c) {if(c=='-') f=-1; c=getchar();}
   while(48<=c&&c<=57) v=(v<<3)+v+v+c-48,c=getchar();
   return v*f;
}

void pushup(node &p,node l,node r)
{
    p.l=l.l|(l.u&p.U&r.l&p.D&l.d);
    p.r=r.r|(r.u&p.U&l.r&p.D&r.d);
    p.u=(l.u&p.U&r.u)|(l.q&p.D&r.p);
    p.d=(l.d&p.D&r.d)|(l.p&p.U&r.q);
    p.q=(l.u&p.U&r.q)|(l.q&p.D&r.d);
    p.p=(l.d&p.D&r.p)|(l.p&p.U&r.u);
}

void build(int l,int r,int p)
{
    if(l==r)
    {
        t[p].U=t[p].D=t[p].u=t[p].d=1;
        return;
    }
    int mid=(l+r)>>1;
    build(l,mid,p<<1);
    build(mid+1,r,p<<1|1);
}

void updater(int l,int r,int x,int y,int v,int p)
{
    int mid=(l+r)>>1;
    if(x==mid)
    {
        if(y==1) t[p].U=v;
         else t[p].D=v;
         pushup(t[p],t[p<<1],t[p<<1|1]);
         return;
    }
    if(x<=mid) updater(l,mid,x,y,v,p<<1);
     else updater(mid+1,r,x,y,v,p<<1|1);
    pushup(t[p],t[p<<1],t[p<<1|1]);
}

void updatec(int l,int r,int x,int v,int p)
{
    int mid=(l+r)>>1;
    if(l==r)
    {
        t[p].l=t[p].r=t[p].p=t[p].q=v;
        return;
    }
    if(x<=mid) updatec(l,mid,x,v,p<<1);
     else updatec(mid+1,r,x,v,p<<1|1);
    pushup(t[p],t[p<<1],t[p<<1|1]);
}

node query(int l,int r,int x,int y,int p)
{
    int mid=(l+r)>>1;
    if(x<=l&&r<=y) return t[p];
    if(y<=mid) return query(l,mid,x,y,p<<1);
     else if(x>mid) return query(mid+1,r,x,y,p<<1|1);
      else
      {
          node tmp=t[p];
        pushup(tmp,query(l,mid,x,y,p<<1),query(mid+1,r,x,y,p<<1|1));
        return tmp;
      }
}

int main()
{
    freopen("bzoj1018.in","r",stdin);
    freopen("bzoj1018.out","w",stdout);
    scanf("%d",&n);
    build(1,n,1);
    int r1,r2,c1,c2; 
    while(scanf("%s",s)!=EOF)
    {
        if(s[0]=='E') break;
        scanf("%d%d%d%d",&r1,&c1,&r2,&c2);
        if(c1>c2)
        {
            swap(c1,c2);
            swap(r1,r2);
        }
        if(s[0]=='O')
        {
            if(r1==r2) updater(1,n,c1,r1,1,1);
             else updatec(1,n,c1,1,1);
        }
        if(s[0]=='C')
        {
            if(r1==r2) updater(1,n,c1,r1,0,1);
             else updatec(1,n,c1,0,1);
        }
        if(s[0]=='A')
        {
            node l=query(1,n,1,c1,1),
                 x=query(1,n,c1,c2,1),
                 r=query(1,n,c2,n,1);
            int ans;
            
            if(r1==1&&r2==1)
                 ans=x.u|(l.r&x.p)|(x.q&r.l)|(l.r&x.d&r.l);
            
            if(r1==1&&r2==2)
                ans=x.q|(l.r&x.d)|(x.u&r.l)|(l.r&x.p&r.l); 
                
            if(r1==2&&r2==1)
                ans=x.p|(l.r&x.u)|(x.d&r.l)|(l.r&x.q&r.l);
            
            if(r1==2&&r2==2)
                ans=x.d|(l.r&x.q)|(x.p&r.l)|(l.r&x.u&r.l);
            if(ans) printf("Y\n");
             else printf("N\n");
        }        
    }      
    return 0;
}