BZOJ 3261 最大异或和 可持久化Trie

前两天被带修改的可持久化线段树给搞毛了,弃疗了QaQ,数据结构搞累了准备换一换其它的。然而,还是遇到了数据结构。0.0
可持久化Trie——如果理解了主席树这个也就很好理解了,建树的方式是相同的,对于这道题来说也是logn的空间来新建一棵树。然而,为了减少不必要的空间浪费,选择了指针——很少用真正的指针,导致无限的调试。
感觉自己的代码还是臃肿,比起piupiuqqq神犇的代码。。简直没得比。

对于题目要求一个后缀异或和suf[p] ^ x (l<=p<=r)的最大值。然而后缀异或和不好维护,因为有插入操作,所以我们可以维护前缀异或和。每次求的可以是suf[p] ^ x = pre[n] ^ pre[p2] ^ x (l-1 <= p2 <= r-1)。x^pre[n]可以直接知道,那么我们的任务就十分明确了:要找到一个pre[p2]使得pre[p2]^x2最大(x2=x^pre[n])。使用可持久化Trie,维护前缀和,在以root[i]为根的Trie中储存pre[1]~pre[i]。范围[l-2,r-1]等同于存在于root[r-1]中而不在root[l-2]中。需要注意边界l-2<0的情况。

QaQ 臃肿的100+的代码。

// 2^24 > 10^7

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
using namespace std;

int get ()
{
    char c = getchar(); int x = 0;
    while (c < '0' || c > '9') c = getchar();
    while (c <= '9' && c >= '0') x = x*10+c-48, c = getchar();
    return x;
}
int getcmd ()
{
    char c = getchar();
    while (c != 'A' && c != 'Q') c = getchar();
    return c == 'A';
}

struct node 
{
    int tot; node *lc, *rc;
    node () {tot = 0, lc = rc = NULL;}
};

struct PersistableTrie
{
    int cnt; node *root[600006];

    PersistableTrie () {root[0]=NULL;}

    void insert (node *rt, node *nrt, int x, int t)
    {
        node *rtlc = rt==NULL ? NULL : rt->lc;
        node *rtrc = rt==NULL ? NULL : rt->rc;
        nrt->tot = 1+(rt==NULL?0:rt->tot);
        if (!t) return ;
        if (x&t) 
        {
            nrt->rc = rtrc;
            insert (rtlc, nrt->lc=new node, x, t>>1);
        }
        else
        {
            nrt->lc = rtlc; 
            insert (rtrc, nrt->rc=new node, x, t>>1);
        }
    }
    void insert (int x)
    {
        root[++cnt] = new node;
        insert (root[cnt-1], root[cnt], x, 1<<24);
    }

    int  query (node *rt, node *nrt, int x, int t)
    {
        if (!t) return 0;
        if (nrt == NULL) return x;
        node *rtlc = rt==NULL ? NULL : rt->lc;
        node *rtrc = rt==NULL ? NULL : rt->rc;
        int ltot  =  rtlc==NULL ? 0 :  rtlc->tot;
        int rtot  =  rtrc==NULL ? 0 :  rtrc->tot;
        int nltot = nrt->lc==NULL ? 0 : nrt->lc->tot;
        int nrtot = nrt->rc==NULL ? 0 : nrt->rc->tot;
        if (x&t)
        {
            if (nrtot-rtot) 
                return t+query(rtrc,nrt->rc,x,t>>1);
            return query(rtlc,nrt->lc,x,t>>1);
        }
        else
        {
            if (nltot-ltot)
                return t+query(rtlc,nrt->lc,x,t>>1);
            return query(rtrc,nrt->rc,x,t>>1);
        }
    }
    void query (int l, int r, int x)
    {
        if (l==1) l = 2;
        int res = query (root[l-2], root[r-1], x, 1<<24);
        char s[15]; int len = 0;
        if (!res) putchar ('0');
        while (res) s[len++] = (res%10)+48, res /= 10;
        while (len) putchar (s[--len]);
        putchar ('\n');
    }
}Solve;

int n, m, sumxor;

int main ()
{
    n = get(), m = get();       
    for (int i = 1; i <= n; i++) Solve.insert(sumxor=sumxor^get());
    for (int i = 1; i <= m; i++)
    {
        int cmd = getcmd();
        if (cmd) Solve.insert(sumxor=sumxor^get());
        else
        {
            int l = get(), r = get(), x = get();
            Solve.query(l, r, x^sumxor);
        }
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值