bzoj4184-shallot

题目

共\(n\)次操作，每次向数集中插入或删除一个数，并询问此时的最大异或和。

分析

看到最大异或和就想到线性基，但是线性基是显然不支持删除的。
题目所需要的就是按时间顺序维护一个线性基。对于这一类插入，删除，整体查询的离线问题，我们可以用线段树（同理cdq分治）解决。
通过离线询问，我们可以记录每个数的插入和删除时间，用线段树维护时间线。对于一个数，我们在线段树上他插入时间和删除时间的区间中插入这个数，然后dfs线段树，在路上累计节点的信息，到达每个叶子节点的时候就计算并输出答案。这里的时间复杂度为\(O(f(n) \* n \* logn)\)。
所以在这道题中，我们在线段树上在每个数出现和完全消失（所有相同的数都被删除）的时间段插入这个数，一次dfs，动态插入数到线性基中，到达叶子节点计算答案，这样我们以\(O(logn)\)的时间代价去除了删除操作。

代码

#include<cstdio>
#include<cctype>
#include<algorithm>
#include<map>
#include<list>
#include<cstring>
using namespace std;
int read() {
    int x=0,f=1;
    char c=getchar();
    for (;!isdigit(c);c=getchar()) if (c=='-') f=-1;
    for (;isdigit(c);c=getchar()) x=x*10+c-'0';
    return x*f;
}
const int maxn=5e5+10;
const int maxj=31;
int t[maxn*3];
struct lb {
    int p[maxj];
    void insert(int x) {
        for (register int j=maxj-1;~j;--j) if (x&(1<<j)) {
            if (!p[j]) {
                for (register int i=j-1;i>=0;--i) if (x&(1<<i)) x^=p[i];
                for (register int i=j+1;i<maxj;++i) if (p[i]&(1<<j)) p[i]^=x;
                p[j]=x;
                break;
            } else x^=p[j];
        }   
    }
    int get() {
        int ret=0;
        for (register int j=0;j<maxj;++j) ret^=p[j];
        return ret;
    }
};
void print(int x) {
    if (!x) {
        puts("0");
        return;
    }
    static char s[12];
    int tot=0;
    while (x) s[++tot]=x%10+'0',x/=10;
    while (tot--) putchar(s[tot+1]);
    puts("");
}
struct edge {
    int v,nxt;
} e[maxn<<4];
int tot=0;
void dfs(int now,int l,int r,lb b) {
    for (register int j=t[now];j;j=e[j].nxt) b.insert(e[j].v);
    if (l==r) {
        print(b.get());
        return;
    }
    int mid=(l+r)>>1;
    dfs(now<<1,l,mid,b);
    dfs(now<<1|1,mid+1,r,b);
}
void insert(int x,int L,int R,int l,int r,int d) {
    if (L==l && R==r) {
        e[++tot]=(edge){d,t[x]};
        t[x]=tot;
        return;
    }
    int mid=(L+R)>>1;
    if (r<=mid) insert(x<<1,L,mid,l,r,d); else 
    if (l>mid) insert(x<<1|1,mid+1,R,l,r,d); else {
        insert(x<<1,L,mid,l,mid,d);
        insert(x<<1|1,mid+1,R,mid+1,r,d);
    }
}
int a[maxn],bas[maxn],b[maxn],st[maxn];
int main() {
    #ifndef ONLINE_JUDGE
        freopen("test.in","r",stdin);
        freopen("my.out","w",stdout);
    #endif
    int n=read(),tk=0;
    for (register int i=1;i<=n;++i) {
        a[i]=read();
        if (a[i]>0) b[++tk]=a[i];
    }
    sort(b+1,b+tk+1);
    tk=unique(b+1,b+tk+1)-b-1;
    for (register int i=1;i<=n;++i) {
        if (a[i]>0) a[i]=lower_bound(b+1,b+tk+1,a[i])-b; else a[i]=-(lower_bound(b+1,b+tk+1,-a[i])-b);
        if (a[i]>0) {
            if (!bas[a[i]]) st[a[i]]=i;
            bas[a[i]]++;
        } else {
            a[i]=-a[i];
            bas[a[i]]--;
            if (!bas[a[i]]) insert(1,1,n,st[a[i]],i-1,b[a[i]]),st[a[i]]=0;
        }
    }
    for (register int i=1;i<=tk;++i) if (st[i]) insert(1,1,n,st[i],n,b[i]);
    dfs(1,1,n,lb());
}

转载于:https://www.cnblogs.com/owenyu/p/6724742.html