[hihoCoder1236 Scores 2015BeijingOnline]简单粗暴的分块+简单粗暴的bitset

题意：50000个5维向量，50000次询问每一维都不大于某一向量的向量个数，强制在线。

思路：做完这题才知道bitset效率这么高，自己本地测试了下1s可以操作10¹⁰个bit，orz简单粗暴

• 令S(i)表示第i维比当前向量的i维小的向量集，则答案为count(∩S(i))，0≤i‹5
• 每个向量都可以和一个id绑定（取下标就行了），绑定后就可以考虑bitset了。分别按每一维排序，每隔√N个位置预处理下当前位置之前的id的bitset
• 查询时，对每一维二分得到最大位置，然后用1个预处理的结果+最多√N次暴力单bit插入就能得到S(i)。然后就是and和count了。

 

#include <bits/stdc++.h>
using namespace std;
#ifndef ONLINE_JUDGE
    #include "local.h"
#endif
#define X first
#define Y second
#define pb(x) push_back(x)
#define mp(x, y) make_pair(x, y)
#define all(a) (a).begin(), (a).end()
#define mset(a, x) memset(a, x, sizeof(a))
#define mcpy(a, b) memcpy(a, b, sizeof(a))
typedef long long ll;
template<typename T>bool umax(T&a, const T&b){return b<=a?false:(a=b,true);}
template<typename T>bool umin(T&a, const T&b){return b>=a?false:(a=b,true);}

int k;
struct Node {
    int a[5];
    int id;
    void read() {
        for (int i = 0; i < 5; i ++) {
            scanf("%d", a + i);
        }
    }
    bool operator < (const Node &that) const {
        return a[k] < that.a[k];
    }
};
Node stu[1 << 16];
int r[5][1 << 16], id[5][1 << 16];
bitset<1 << 16> bs[5][1 << 8];

int main() {
#ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
    //freopen("out.txt", "w", stdout);
#endif // ONLINE_JUDGE
    int T, n, q, m;
    cin >> T;
    while (T --) {
        cin >> n >> m;
        for (int i = 0; i < n; i ++) {
            stu[i].read();
            stu[i].id = i;
        }
        int L = sqrt(n + 0.5);
        for (int i = 0; i < 5; i ++) {
            k = i;
            sort(stu, stu + n);
            bitset<1 << 16> bbs;
            for (int j = 0; j < n; j ++) {
                r[i][j] = stu[j].a[i];
                id[i][j] = stu[j].id;
                bbs[stu[j].id] = 1;
                if (j % L == L - 1 || j == n - 1)  {
                    bs[i][j / L] = bbs;
                }
            }
        }
        int lastans = 0;
        cin >> q;
        while (q --) {
            Node qry;
            qry.read();
            bitset<1 << 16> ans;
            for (int i = 0; i < 5; i ++) {
                qry.a[i] ^= lastans;
                bitset<1 << 16> buf;
                int pos = upper_bound(r[i], r[i] + n, qry.a[i]) - r[i];
                if (pos / L) buf = bs[i][pos / L - 1];
                while (pos % L) buf[id[i][-- pos]] = 1;
                if (i) ans &= buf;
                else ans = buf;
            }
            printf("%d\n", lastans = ans.count());
        }
    }
    return 0;
}