Cod sursa(job #2887373)

Utilizator Vlad_AnicaAnica-Popa Vlad-Ioan Vlad_Anica Data 9 aprilie 2022 14:32:28
Problema A+B Scor 0
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 5.77 kb
#include <bits/stdc++.h>
#define ld long double
#pragma GCC optimize("Ofast")
#define lsb(x) ((x) & (-(x)))

using namespace std;

const int K = 700;

struct Query {
    int l, r, ind;

    bool operator < (Query &other) const {
        if(l / K != other.l / K)
            return (l < other.l);

        return r > other.r;
    }
} v[500005];

mt19937 gen(time(0));
uniform_int_distribution <int> rng;

int n, x;
int l, r;
int nr, op;

int s[1000005];
int tmp[1500005], vals[500005];
unordered_map <int, int> cod;
int ans[500005];
int nxt[1500005], prv[1500005];
int nxt1[1000005], prv1[1000005];
int lst[500005];
int aib[500005];

void update(int ind, int val) {
    if(!ind)
        return;
    for(; ind <= n; ind += lsb(ind))
        aib[ind] += val;
}

int query(int ind) {
    int ans = 0;
    for(; ind; ind -= lsb(ind))
        ans += aib[ind];

    return ans;
}

void print() {
    cout << "Print\n";
    for(int i = 1; i <= nr + n; i++)
        cout << i - nr << " " << prv[i] - nr << " " << nxt[i] - nr << "\n";
}

void del(int poz, int st) {
    poz += nr;

    if(st) {
        if(nxt[poz] <= nr + n) {
            //w.erase(w.find(nxt[poz] - poz));
            update(nxt[poz] - poz, -1);

            if(prv[poz] >= nr) {
                //w.insert(nxt[poz] - prv[poz]);
                update(nxt[poz] - prv[poz], 1);
            }
        }
        prv[nxt[poz]] = prv[poz];
        nxt[prv[poz]] = nxt[poz];
    }
    else {
        if(prv[poz] >= nr) {
            //w.erase(w.find(poz - prv[poz]));
            update(poz - prv[poz], -1);

            if(nxt[poz] <= nr + n) {
                //w.insert(nxt[poz] - prv[poz]);
                update(nxt[poz] - prv[poz], 1);
            }
        }

        //cout << prv[poz] << " " << nxt[poz] << " " << poz << "\n";
        nxt[prv[poz]] = nxt[poz];
        prv[nxt[poz]] = prv[poz];
    }

    prv[poz] = nxt[poz] = -1;
}

void add(int poz, int st) {
    //cout << poz << " add\n";
    poz += nr;

    if(st) {
        int pz = prv1[poz];
        //cout << poz << " " << pz << "\n";
        if(pz >= nr && nxt[pz] <= nr + n) {
            //w.erase(w.find(nxt[pz] - pz));
            update(nxt[pz] - pz, -1);
        }
        prv[poz] = pz;
        nxt[poz] = nxt[pz];
        prv[nxt[pz]] = poz;
        if(nxt[pz] <= nr + n) {
            //w.insert(nxt[pz] - poz);
            update(nxt[pz] - poz, 1);
        }
        nxt[pz] = poz;
    }
    else {
        int pz = nxt1[poz];
        if(pz <= nr + n && prv[pz] >= nr) {
           // w.erase(w.find(pz - prv[pz]));
            update(pz - prv[pz], -1);
        }
        nxt[prv[pz]] = poz;
        if(prv[pz] >= nr) {
            //w.insert(poz - prv[pz]);
            update(poz - prv[pz], 1);
        }
        nxt[poz] = pz;
        prv[poz] = prv[pz];
        prv[pz] = poz;
    }

    //print();
}

bool kek = 0;

void solve() {
    if(kek) {
        cin >> n;
    } else {
        n = 500000;
    }

    ld t1 = clock();
    cod.clear();
    vals[0] = 0;
    for(int i = 1; i <= n; i++) {
        if(kek) {
            cin >> x;
        } else {
            x = abs(rng(gen));
        }
        s[i] = s[i - 1] ^ x;
        vals[i] = s[i];
    }

    sort(vals + 1, vals + n + 1);

    nr = 0;

    cod[0] = ++nr;

    for(int i = 1; i <= n; i++) {
        if(vals[i] != vals[i - 1])
            cod[vals[i]] = ++nr;
    }

    nr++;

    for(int i = nr; i <= nr + n; i++) {
        tmp[i] = cod[s[i - nr]];
    }

    for(int i = nr; i <= nr + n; i++)
        s[i] = tmp[i];

    int q;
    if(kek) {
        cin >> q;
    }
    else {
        q = 500000;
    }

    for(int i = 1; i <= q; i++) {
        if(kek) {
            cin >> v[i].l >> v[i].r;
        } else {
            v[i].l = abs(rng(gen)) % n;
            v[i].r = abs(rng(gen)) % (n - v[i].l) + v[i].l;
        }
        v[i].ind = i, v[i].l++;
    }

    sort(v + 1, v + q + 1);

    for(int i = 0; i <= nr; i++)
        lst[i] = i;

    for(int i = 0; i <= n + 2 * nr; i++)
        prv[i] = nxt[i] = -1;

    for(int i = nr; i <= nr + n; i++) {

        //cout << nr << " " << i << " " << s[i] << " " << lst[s[i]] << "\n";

        nxt[lst[s[i]]] = i;
        prv[i] = lst[s[i]];
        lst[s[i]] = i;
    }

    for(int i = n + nr + 1; i <= n + 2 * nr - 1; i++) {
        prv[i] = lst[i - n - nr];
        nxt[lst[i - n - nr]] = i;
    }

    for(int i = 0; i <= n + 2 * nr; i++)
        nxt1[i] = nxt[i], prv1[i] = prv[i];

    //print();

    //print();

    int l = 0, r = -1;
    for(int i = 1; i <= q; i++) {
        //cout << v[i].l << " " << v[i].r << " xdddddddddd\n";
        while(l > v[i].l) {
            l--;
            del(l, 1);
        }
        while(v[i].r > r) {
            r++;
            del(r, 0);
        }
        while(l < v[i].l) {
            add(l, 1);
            l++;
        }
        while(v[i].r < r) {
            add(r, 0);
            r--;
        }

        //print();

        int st = 1, dr = n, mid, ansu = 0;
        int lg = log2(n);

        for(int i = lg; i >= 0; i--) {
            if(!aib[ansu | (1 << i)])
                ansu |= (1 << i);
        }

        ans[v[i].ind] = (ansu > n ? -1 : ansu + 1);
        //print();
    }

    if(kek) {
        for(int i = 1; i <= q; i++)
            cout << ans[i] << "\n";
    }
    else {
        ld t2 = clock();
        cout << (t2 - t1) / CLOCKS_PER_SEC << "\n";
        cout << op << "\n";
    }
}

int main() {
    ios_base :: sync_with_stdio(false); cin.tie(0); cout.tie(0);
    srand(time(0));

    int T = 100;

    kek = 1;

    if(kek)
        T = 1;

    while(T--) {
        solve();
    }

    return 0;
}