hihocoder 1236

题意:给出N个人的5个科目分数,给出q个询问~~,每次给你5个科目的具体分数。求一共有多少个人对应的5个科目都小于等于你的科目分数~~


渣渣了~~~


做法:分别单独把科目拿来看,看都有具体那些人满足条件。然后合并找出满足这五个条件的个人~~

首先将单科目的人按照从小到达的顺序排序,然后进行二分查找就可以知道一共有多少个人满足条件。

但是我们需要知道有那些人满足条件而不是数量~~。

所以我们维护一个bitset[i]:=前i个人满足条件时,都有谁~~~bitset[i][j]==1时有人~~


最后求出的五种情况&一下就可以知道都有谁剩下了~~


但是bitset<50001>b[5][50000];   ///b[i][j]:=第i维、前j个人都是谁~~~~~这样会爆内存的。实际上这里是TLE~~~

所以这里采取分块的方法进行优化~~  

即bitset<50001>b[5][250];     ///b[i][j]:=第i维、前j个块都有谁~~

于是当二分找到的位置是idx,位于块k

利用前k-1块的信息再加上当前块满足条件的几个人就可以得到相应集合了~~

#include <algorithm>
#include <iostream>
#include<string.h>
#include <fstream>
#include <math.h>
#include <vector>
#include <cstdio>
#include <string>
#include <queue>
#include <stack>
#include <map>
#include<bitset>
#include <set>
#define exp 1e-8
#define fi first
#define se second
#define ll long long
#define INF 0x3f3f3f3f
#define pb(a) push_back(a)
#define mp(a,b) make_pair(a,b)
#define all(a) a.begin(),a.end()
#define mm(a,b) memset(a,b,sizeof(a));
#define for1(a,b) for(int a=1;a<=b;a++)//1---(b)
#define rep(a,b,c) for(int a=b;a<=c;a++)//b---c
#define repp(a,b,c)for(int a=b;a>=c;a--)///
#define stl(c,itr) for(__typeof((c).begin()) itr=(c).begin();itr!=(c).end();itr++)
using namespace std;
void bug(string m="here"){cout<<m<<endl;}
template<typename __ll> inline void READ(__ll &m){__ll 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();}m=x*f;}
template<typename __ll>inline void read(__ll &m){READ(m);}
template<typename __ll>inline void read(__ll &m,__ll &a){READ(m);READ(a);}
template<typename __ll>inline void read(__ll &m,__ll &a,__ll &b){READ(m);READ(a);READ(b);}
template<typename __ll>inline void read(__ll &m,__ll &a,__ll &b,__ll &c){READ(m);READ(a);READ(b);READ(c);}
template<typename __ll>inline void read(__ll &m,__ll &a,__ll &b,__ll &c,__ll &d){READ(m);READ(a);READ(b);READ(c);read(d);}
template < class T > T gcd(T a, T b) { return b ? gcd(b, a % b) : a; }
template < class T > T lcm(T a, T b) { return a / gcd(a, b) * b; }
template < class T > inline void rmin(T &a, const T &b) { if(a > b) a = b; }
template < class T > inline void rmax(T &a, const T &b) { if(a < b) a = b; }
template < class T > T pow(T a, T b) { T r = 1; while(b > 0) { if(b & 1) r = r * a; a = a * a; b /= 2; } return r; }
template < class T > T pow(T a, T b, T mod) { T r = 1; while(b > 0) { if(b & 1) r = r * a % mod; a = a * a % mod; b /= 2; } return r; }

const int maxn=50001;
bitset<maxn> b[6][300];
bitset<maxn>ans[6];
int now[6];
struct DAT
{
    int val,idx;
    DAT(){}
    DAT(int k){val=k,idx=0;}
    bool operator < (const DAT &rhs)const{
        return val<rhs.val;
    }
}dat[6][maxn];
int n,m,q;
int len,num;
int l[300],r[300];
int main()
{
    int cas;
    read(cas);
    while(cas--)
    {
        read(n,m);
        for1(j,n)for1(i,5)
        {
            read(dat[i][j].val);
            dat[i][j].idx=j;
        }
        for1(i,5)sort(dat[i]+1,dat[i]+n+1);
        for(int i=1;i<=5;i++)
            for(int j=1;j<250;j++)
                b[i][j].reset();

        len=sqrt(n);
        num=n/len;
        if(n%len)num++;
        for(int i=1;i<=num;i++)
            l[i]=(i-1)*len+1,r[i]=i*len;
        r[num]=n;

        for(int i=1;i<=5;i++)
        {
            for(int j=1;j<=num;j++)
            {
                b[i][j]|=b[i][j-1];
                for(int k=l[j];k<=r[j];k++)
                    b[i][j][dat[i][k].idx]=1;

            }
        }
        read(q);
        int lastans=0;
        while(q--)
        {
            for1(i,5)read(now[i]);
            for1(i,5)now[i]^=lastans;
            for1(i,5)ans[i].reset();
            for(int i=1;i<=5;i++)
            {
                int idx=upper_bound(dat[i]+1,dat[i]+1+n,DAT(now[i]))-dat[i]-1;
                if(idx==0)continue;
                int belong=(idx-1)/len+1;
                ans[i]|=b[i][belong-1];
                for(int k=l[belong];k<=idx;k++)
                    ans[i][dat[i][k].idx]=1;
            }
            for(int i=2;i<=5;i++)
                ans[1]&=ans[i];
            lastans=ans[1].count();
            printf("%d\n",lastans);
        }
    }
}



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值