羅区[P4688] [Ynoi2016]ウサギの穴(ビットセット、MOチーム)ダウン

トピックリンク
まずトラックYnoi
もちろん答えは、* 3数3の間隔を引くことによって、各問い合わせや公人の長さで
、その後3つの区間ビット単位、数は公人であれば、あなたは、各区間のMoの+ビットセットメモリの状態を持つチームのようなものができ、ライン上。
しかし、今の数で、それぞれの番号を保存するだけでビットセットが表示されますが、また、出現回数を保存します。
その後、我々は、各番号の数が表示され、ことがわかった(= N \)\
ので、グリッド内の各離散的なデジタルアカウントのビットセット後に期待しています。
覚えておいてください\(SA \)で基数ソートそれを?これはセクション可能\(N- \)固定された取付位置を有するデジタル時間間隔を加えています。
だから、チームはミズーリことができるようになります
しかし、データの範囲、もし劣らず、オープンを見て\(1E5 \)の長さ\(1E5 \)ビットセットああ。
これは、答えは三つのグループ、リスト上のグループへの1つに分割した、重要ではありません。

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <bitset>
#include <cmath>
#include <iostream>
using namespace std;
const int MAXN = 100010;
const int MAXM = 34010;
bitset <MAXN> ans[MAXM], now;
int n, m, a[MAXN];
inline int read(){
    int s = 0;
    char ch = getchar();
    while(ch < '0' || ch > '9') ch = getchar();
    while(ch >= '0' && ch <= '9'){ s = s * 10 + ch - '0'; ch = getchar(); }
    return s;
}
int Q;
struct lsh{
    int val, id;
    int operator < (const lsh A) const{
        return val < A.val;
    }
}p[MAXN];
struct ask{
    int l, r, id;
    int operator < (const ask A) const{
        return l / Q == A.l / Q ? r < A.r : l < A.l;
    }
}q[MAXM * 3];
int val[MAXN], cnt, sum[MAXN], v[MAXN], all, Ans[MAXM];
void work(){
    sort(q + 1, q + all + 1);
    now.reset();
    memset(v, 0, sizeof v);
    int l = 1, r = 1; now.set(val[1]); v[val[1]] = 1;
    for(int i = 1; i <= all; ++i){
        while(r < q[i].r){ ++r; now.set(val[r] - v[val[r]]); ++v[val[r]]; }
        while(l > q[i].l){ --l; now.set(val[l] - v[val[l]]); ++v[val[l]]; }
        while(r > q[i].r){ --v[val[r]]; now.set(val[r] - v[val[r]], 0); --r; }
        while(l < q[i].l){ --v[val[l]]; now.set(val[l] - v[val[l]], 0); ++l; }
        ans[q[i].id] &= now;
    }
}
int main(){
    n = read(); m = read(); Q = sqrt(n);
    for(int i = 1; i <= n; ++i)
        p[i].val = read(), p[i].id = i;
    for(int i = 1; i <= 33337; ++i)
       ans[i].set();
    sort(p + 1, p + n + 1);
    for(int i = 1; i <= n; ++i)
        if(p[i].val != p[i - 1].val)
            val[p[i].id] = ++cnt;
        else val[p[i].id] = cnt;
    for(int i = 1; i <= n; ++i)
        ++sum[val[i]];
    for(int i = 1; i <= cnt; ++i)
        sum[i] += sum[i - 1];
    for(int i = 1; i <= n; ++i)
        val[i] = sum[val[i]];
    int o = m / 3;
    if(o){
        for(int qqc = 1; qqc <= 2; ++qqc){
            cnt = 0;
            for(int i = 1; i <= o; ++i){
                q[++cnt].id = i; q[cnt].l = read(); q[cnt].r = read(); Ans[i] += q[cnt].r - q[cnt].l + 1;
                q[++cnt].id = i; q[cnt].l = read(); q[cnt].r = read(); Ans[i] += q[cnt].r - q[cnt].l + 1;
                q[++cnt].id = i; q[cnt].l = read(); q[cnt].r = read(); Ans[i] += q[cnt].r - q[cnt].l + 1;
            }
            all = o * 3; work();
            for(int i = 1; i <= o; ++i){
                printf("%d\n", Ans[i] - int(ans[i].count()) * 3);
                Ans[i] = 0; ans[i].set();
            }
        }
    }
    m -= o * 2; cnt = 0;
    for(int i = 1; i <= m; ++i){
        q[++cnt].id = i; q[cnt].l = read(); q[cnt].r = read(); Ans[i] += q[cnt].r - q[cnt].l + 1;
        q[++cnt].id = i; q[cnt].l = read(); q[cnt].r = read(); Ans[i] += q[cnt].r - q[cnt].l + 1;
        q[++cnt].id = i; q[cnt].l = read(); q[cnt].r = read(); Ans[i] += q[cnt].r - q[cnt].l + 1;
    }
    all = m * 3; work();
    for(int i = 1; i <= m; ++i)
        printf("%d\n", Ans[i] - int(ans[i].count()) * 3);
    return 0;
}

おすすめ

転載: www.cnblogs.com/Qihoo360/p/11026356.html