BZOJ3658 : Jabberwocky

考虑将某线段下方的点取走：

将所有点从低到高排序

每扫描到一条水平线，对于上面每个点，找到它下面同色的前驱后继，统计中间点的个数

然后再把线上所有点插入数据结构中

最后再统计相邻的同色的点之间的点个数

用动态开点的权值线段树+树状数组维护，时间复杂度$O(nlog n)$。

考虑将某线段上方的点取走：

把扫描线的顺序反过来即可

注意特判出现颜色数没达到k的情况。

#include<cstdio>
#include<algorithm>
using namespace std;
const int N=100010,M=1800000;
int C,n,k,i,j,ans,bit[N],b[N],vis[N];
struct P{int x,y,c;}a[N];
inline bool cmpc(P a,P b){return a.c==b.c?a.x<b.x:a.c<b.c;}
inline bool cmpy(P a,P b){return a.y<b.y;}
inline void up(int x){if(ans<x)ans=x;}
inline int lower(int x){
  int l=1,r=n,t,mid;
  while(l<=r)if(b[mid=(l+r)>>1]<=x)l=(t=mid)+1;else r=mid-1;
  return t;
}
inline void add(int x){for(;x<=n;x+=x&-x)bit[x]++;}
inline int ask(int x){int t=0;for(;x;x-=x&-x)t+=bit[x];return t;}
inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';}
int T[N],vl[M],vr[M],l[M],r[M],tot,pre,nxt;
void ins(int&x,int a,int b,int c){
  if(!x)x=++tot;
  if(a==b){vl[x]=vr[x]=a;return;}
  int mid=(a+b)>>1;
  if(c<=mid)ins(l[x],a,mid,c);else ins(r[x],mid+1,b,c);
  vl[x]=vl[l[x]]?vl[l[x]]:vl[r[x]];
  vr[x]=vr[r[x]]?vr[r[x]]:vr[l[x]];
}
void getl(int x,int a,int b,int c){
  if(!x||nxt)return;
  if(c<=a){nxt=vl[x];return;}
  int mid=(a+b)>>1;
  if(c<=mid)getl(l[x],a,mid,c);
  getl(r[x],mid+1,b,c);
}
void getr(int x,int a,int b,int d){
  if(!x||pre)return;
  if(b<=d){pre=vr[x];return;}
  int mid=(a+b)>>1;
  if(d>mid)getr(r[x],mid+1,b,d);
  getr(l[x],a,mid,d);
}
int main(){
  for(read(C);C--;printf("%d
",ans)){
    read(n),read(k);
    for(i=1;i<=k;i++)vis[i]=0;
    for(ans=0,i=1;i<=n;i++)read(a[i].x),read(a[i].y),read(a[i].c),b[i]=a[i].x,vis[a[i].c]=1;
    for(sort(b+1,b+n+1),sort(a+1,a+n+1,cmpy),i=1;i<=n;i++)a[i].x=lower(a[i].x);
    for(i=1;i<=n;i++)bit[i]=0;
    for(i=1;i<=k;i++)if(!vis[i]){ans=n;break;}
    if(ans)continue;
    for(i=1;i<=tot;i++)vl[i]=vr[i]=l[i]=r[i]=0;tot=0;
    for(i=1;i<=k;i++)T[i]=0,ins(T[i],0,n+1,0),ins(T[i],0,n+1,n+1);
    for(i=1;i<=n;i=j){
      for(j=i;j<=n&&a[j].y==a[i].y;j++){
        pre=nxt=0,getl(T[a[j].c],0,n+1,a[j].x),getr(T[a[j].c],0,n+1,a[j].x);
        up(ask(nxt-1)-ask(pre));
      }
      for(j=i;j<=n&&a[j].y==a[i].y;j++)ins(T[a[j].c],0,n+1,a[j].x),add(a[j].x);
    }
    for(i=1;i<=n;i++)bit[i]=0;
    for(i=1;i<=tot;i++)vl[i]=vr[i]=l[i]=r[i]=0;tot=0;
    for(i=1;i<=k;i++)T[i]=0,ins(T[i],0,n+1,0),ins(T[i],0,n+1,n+1);
    for(i=n;i;i=j){
      for(j=i;j&&a[j].y==a[i].y;j--){
        pre=nxt=0,getl(T[a[j].c],0,n+1,a[j].x),getr(T[a[j].c],0,n+1,a[j].x);
        up(ask(nxt-1)-ask(pre));
      }
      for(j=i;j&&a[j].y==a[i].y;j--)ins(T[a[j].c],0,n+1,a[j].x),add(a[j].x);
    }
    for(sort(a+1,a+n+1,cmpc),i=1;i<=n;i++)bit[i]=0;
    for(i=1;i<=n;i++)bit[a[i].x]++;
    for(i=2;i<=n;i++)bit[i]+=bit[i-1];
    for(i=1;i<n;i++)if(a[i].c==a[i+1].c)up(bit[a[i+1].x-1]-bit[a[i].x]);
    else up(max(n-bit[a[i].x],bit[a[i+1].x-1]));
    up(max(bit[a[1].x-1],n-bit[a[n].x]));
  }
  return 0;
}