CIA11 毒瘤题(七维偏序,分块+bitset)

首先bitset很强呀x,我们针对每一维,对于每一个点预处理出比他小的集合。然后他的答案就是这k个集合的交。但是O(n^2k)的空间很gg。我们可以分块,对于每一维,只预处理出若干块集合,然后查询的时候找到离他最近的块,用这个块的答案暴力得出正确的集合。应该是 O(nn+n27/32)

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <bitset>
using namespace std;
#define ll long long
#define N 50010
#define nn 233
#define inf 0x3f3f3f3f
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,K,L[8][240],R[8][240],D;
bitset<N>b[8][240],ans[8];
struct point{
    int d[8];
    int& operator[](int x){return d[x];}
    friend bool operator<(point a,point b){return a[D]<b[D];}
}P[N];
struct data{
    int x,id;
    friend bool operator<(data a,data b){return a.x<b.x;}
}a[8][N];
inline void solve(){
    sort(a[D]+1,a[D]+n+1);int now=1;L[D][now]=1;
    for(int i=1;i<=n;++i){
        b[D][now][a[D][i].id]=1;
        while(a[D][i+1].x==a[D][i].x) ++i,b[D][now][a[D][i].id]=1;
        if(i-L[D][now]+1>=nn) R[D][now]=i,L[D][++now]=i+1,b[D][now]=b[D][now-1];
    }R[D][now]=n;L[D][0]=now;
}
inline void get(int x){
    int now=1;ans[D].reset();
    while(now<=L[D][0])
        if(a[D][R[D][now]].x>=x) break;else now++;
    ans[D]=b[D][now-1];
    for(int i=L[D][now];i<=R[D][now];++i){
        if(a[D][i].x<=x) ans[D][a[D][i].id]=1;
        else break;
    }
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();K=read();
    for(int i=1;i<=n;++i){
        for(int k=1;k<=K;++k) a[k][i].x=P[i][k]=read(),a[k][i].id=i;
    }for(int i=1;i<=K;++i) D=i,solve();
    for(int i=1;i<=n;++i){
        for(int k=1;k<=K;++k) D=k,get(P[i][k]);
        for(int k=2;k<=K;++k) ans[1]&=ans[k];
        printf("%d\n",ans[1].count()-1);
    }return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值