【题解】 BZOJ4548 小奇的糖果

本文同步在学弟ZCDHJ的个人博客发布，审核需要一段时间.

传送门

---

考虑题目中获得的糖果并不包含所有的颜色这句话,发现相当于我们可以直接选取某一个颜色强制不能选(这样子一定最优).
然后就可以考虑分开解决上面和下面.
先考虑下面:

1. 枚举颜色然后搞区间(不能包含这一种颜色)
2. 按照横坐标的顺序删点,然后再看删除的点的颜色的区间是否会被更新

以上操作的话,查询区间可以用双向链表,区间求点(和)可以使用树状数组.
差不多就这样子了.

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<iostream>
using namespace std;
inline int gi(){
    int sum=0,f=1;char ch=getchar();
    while(ch>'9' || ch<'0'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9'){sum=(sum<<3)+(sum<<1)+ch-'0';ch=getchar();}
    return f*sum;
}
int ans,n,K,c[500010],x[500010],l[500010],r[500010],last[500010],b[500010];
struct node{
    int id,x,y,z;
}p[500010];
int lowbit(int x){
    return x&-x;
}
void add(int x,int d){
    for(;x<=n+1;x+=lowbit(x))c[x]+=d;
}
int query(int x){
    int res=0;
    for(;x>0;x-=lowbit(x))res+=c[x];
    return res;
}
bool cmp1(node a,node b){
    return a.x<b.x;
}
bool cmp2(node a,node b){
    return a.y<b.y;
}
void update(int l,int r){
    if(l>r)return;
    ans=max(ans,query(r)-query(l-1));
}
void solve(){
    x[0]=0;x[n+1]=n+1;
    memset(c,0,sizeof(c));memset(last,0,sizeof(last));
    sort(p+1,p+n+1,cmp1);
    for(int i=1;i<=n;i++)add(p[i].x,1);
    for(int i=1;i<=n;i++){
        int t=p[i].id,L=last[p[i].z];
        l[t]=L;r[t]=n+1;
        if(L)r[L]=t;update(x[L]+1,x[t]-1);
        last[p[i].z]=t;
    }
    for(int i=1;i<=K;i++)update(x[last[i]]+1,n+1);
    sort(p+1,p+n+1,cmp2);
    for(int i=1,j=1;i<=n;i++){
        int t=p[i].id;
        while(j<=n && p[i].y==p[j].y)add(p[j].x,-1),j++;
        l[r[t]]=l[t];r[l[t]]=r[t];
        update(x[l[t]]+1,x[r[t]]-1);
    }
}
int main(){
#ifndef ONLINE_JUDGE
    freopen("in.in","r",stdin);
    freopen("out.out","w",stdout);
#endif
    int T=gi();
    while(T--){
        ans=0;
        n=gi();K=gi();
        for(int i=1;i<=n;i++){
            p[i].x=gi();p[i].y=gi();p[i].z=gi();
            p[i].id=i;
        }
        for(int i=1;i<=n;i++)b[i]=p[i].x;
        sort(&b[1],&b[n+1]);
        for(int i=1;i<=n;i++){
            p[i].x=lower_bound(b+1,b+n+1,p[i].x)-b;
            x[i]=p[i].x;
        }
        solve();
        for(int i=1;i<=n;i++)p[i].y=-p[i].y;
        solve();
        printf("%d\n",ans);
    }
    return 0;
}