[BZOJ3189][Coci2011][扫描线][线段树]Slika

首先发现LOAD和与其对应的SAVE之间的操作是没用的，所以可以O(m)扫两遍出去SAVE和LOAD操作。

对于这种矩形覆盖的问题，可以用扫面线加线段树处理。
开两棵线段树，用来记录(x+y)为奇数或偶数的时候被覆盖的情况，我们不妨扫描线扫x坐标，那么随着x+1，只要交换一下两颗线段树就可以了。

对于线段树中的一个结点，开一个set记录一下覆盖了这个结点表示的区间的颜色以及时间，以时间为关键字（我直接记录时间，在外面开个数组记录每个时间对应的颜色），这样就相当于区间覆盖，单点查询。

本地跑得贼慢，交到OJ上跑得比被抓的余畅还快

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <set>
#include <vector>
#define N 1010
#define M 100010

using namespace std;

int n,m,k,t;
int S[M],cl[M];
char c[N];
struct stp{
  int op,x1,y1,x2,y2,c,t;
}A[M];
vector<stp> b[N],e[N];
struct seg{
  int l,r;
  set<int,greater<int> > S;
};

struct tree{
  int L,R,Maxx;
  seg T[2][N<<3];
  tree(){ L=0; R=1; }
  void Swap(){ swap(L,R); }
  void insert(int g,int l,int r,int t,int p){
    if(T[p][g].l==l&&T[p][g].r==r) return (void)T[p][g].S.insert(t);
    int mid=T[p][g].l+T[p][g].r>>1;
    if(r<=mid) insert(g<<1,l,r,t,p);
    else if(l>mid) insert(g<<1|1,l,r,t,p);
    else insert(g<<1,l,mid,t,p),insert(g<<1|1,mid+1,r,t,p);
  }
  void erase(int g,int l,int r,int x,int p){
    if(T[p][g].l==l&&T[p][g].r==r) return (void)T[p][g].S.erase(x); 
    int mid=T[p][g].l+T[p][g].r>>1;
    if(r<=mid) erase(g<<1,l,r,x,p);
    else if(l>mid) erase(g<<1|1,l,r,x,p);
    else erase(g<<1,l,mid,x,p),erase(g<<1|1,mid+1,r,x,p);
  }
  void erase(stp x){
    if(((x.x2-x.x1+1)&1)^(x.y1&1))
      erase(1,x.y1,x.y2,x.t,L);
    else
      erase(1,x.y1,x.y2,x.t,R);
  }
  void insert(stp x){
    if(x.y1&1) insert(1,x.y1,x.y2,x.t,L);
    else insert(1,x.y1,x.y2,x.t,R);
  }
  void query(int g,int x,int p){
    if(!T[p][g].S.empty()) Maxx=max(Maxx,*T[p][g].S.begin());
    if(T[p][g].l==T[p][g].r) return ;
    int mid=T[p][g].l+T[p][g].r>>1;
    if(x<=mid) query(g<<1,x,p);
    else query(g<<1|1,x,p);
  }
  int query(int x){
    Maxx=0; query(1,x,(x&1)?L:R);
    return Maxx;
  }
  void Build(int g,int l,int r,int p){
    T[p][g].l=l; T[p][g].r=r; T[p][g].S.clear();
    if(l==r) return ;
    int mid=l+r>>1;
    Build(g<<1,l,mid,p); Build(g<<1|1,mid+1,r,p);
  }
  void build(){
    Build(1,1,n,0); Build(1,1,n,1);
  }
}T;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}

inline void rea(int &x){
  char c=nc(); x=0;
  for(;c>'9'||c<'0';c=nc());for(;c>='0'&&c<='9';x=x*10+c-'0',c=nc());
}

inline void rea(char *x){
  char c=nc(); int len=0;
  for(;c>'Z'||c<'A';c=nc());for(;c>='A'&&c<='Z';x[len++]=c,c=nc());
}

int main(){
  rea(n); rea(k); rea(m); cl[0]=1;
  for(int i=1;i<=m;i++){
    rea(c); A[i].t=i;
    if(c[0]=='S') A[i].op=1;
    else if(c[0]=='L') A[i].op=2,rea(A[i].c);
    else{
      rea(A[i].c),rea(A[i].x1),rea(A[i].y1),rea(A[i].x2),rea(A[i].y2);
      A[i].x1++; A[i].x2++; A[i].y1++; A[i].y2++;
      cl[i]=A[i].c; 
    }
  }
  for(int i=1;i<=m;i++) if(A[i].op==1) S[++t]=i;
  /*for(int i=m;i>0;i--){
    if(A[i].op==0) printf("%d %d\n",A[i].x1,A[i].x2);
    else if(A[i].op==2) i=S[A[i].c];
  }
  return 0;*/
  for(int i=m;i>0;i--){
    if(A[i].op==0){
      b[A[i].x1].push_back(A[i]);
      e[A[i].x2+1].push_back(A[i]);
    }
    else if(A[i].op==2) i=S[A[i].c];
  }
  T.build();
  for(int i=1;i<=n;i++,putchar('\n')){
    T.Swap();
    for(int j=0;j<e[i].size();j++)
      T.erase(e[i][j]);
    for(int j=0;j<b[i].size();j++)
      T.insert(b[i][j]);
    for(int j=1;j<=n;j++)
      printf("%d ",cl[T.query(j)]);
  }
  return 0;
}