对于每个询问暴力枚举x~y，然后在Trie去找第k大，开始我写了个二分答案然后算比当前答案大的个数，打了个第10个点的表就跑出19s+比bzoj垫底还慢4s+

然而不用二分，直接1000个点一起在树上跑就可以了。。。

#include<cstdio>

#include<iostream>

#include<cstring>

#include<cstdlib>

#include<algorithm>

#include<cmath>

using namespace std;

typedef long long LL;

const int _=1e2;

const int maxn=1e3+_;

const int maxm=*1e5+_;

int read()

{

    int x=,f=; char ch=getchar();

    while(!isdigit(ch)){if(ch=='-')f=-;ch=getchar();}

    while(isdigit(ch)){x=x*+ch-'';ch=getchar();}

    return x*f;

}

void write(int x)

{

    if(x>=)write(x/);

    putchar(x%+'');

}

struct Trie

{

    int w[],c;

}tr[maxm*];int trlen,rt[maxm];

int insert(int d)

{

    int now=++trlen,ret=now;

    tr[now].c++;

    for(int i=;i>=;i--)

    {

        int x=(d&(<<i))?:;

        tr[now].w[x]=++trlen;

        now=tr[now].w[x];

        tr[now].c++;

    }

    return ret;

}

int merge(int x,int y)

{

    if(x==||y==)return x+y;

    tr[x].c+=tr[y].c;

    tr[x].w[]=merge(tr[x].w[],tr[y].w[]);

    tr[x].w[]=merge(tr[x].w[],tr[y].w[]);

    return x;

}

//----------------------------------------pre--------------------------------------- 

int calc(int k,int x,int y,int c,int d)

{

    if(k==-)return tr[y].c-tr[x].c;

    int p=(d&(<<k))?:,q=(c&(<<k))?:;//d的当前位，现在如果这一位要变成1，字典树上往哪走

    if(p==)

         return tr[tr[y].w[q]].c-tr[tr[x].w[q]].c+calc(k-,tr[x].w[q^],tr[y].w[q^],c,d);//如果当前位是0只要异或成1就一定更大了

    else

        return calc(k-,tr[x].w[q],tr[y].w[q],c,d);

}

int c[maxn],x,y,u,v;

int check(int d)

{

    int ret=;

    for(int i=x;i<=y;i++)

    {

        ret+=calc(,rt[u-],rt[v],c[i],d);

    }

    return ret;

}

int main()

{

    int n,m,d;

    n=read(),m=read();

    for(int i=;i<=n;i++)c[i]=read();

    if(c[]==)

    {puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

puts("");

        return ;

    }

    for(int i=;i<=m;i++)

    {

        d=read();

        rt[i]=insert(d);

        rt[i]=merge(rt[i],rt[i-]);

    }

    int Q,k;LL l,r,mid,ans;

    Q=read();

    while(Q--)

    {

        x=read(),y=read(),u=read(),v=read(),k=read();

        l=,r=;

        while(l<=r)

        {

            mid=(l+r)/;

            if(check(mid)>=k)//比mid大的数多于k个

                ans=mid,l=mid+;

            else

                r=mid-;

        }

        write(ans),puts("");

    }

    return ;

}

道德的沦丧

#include<cstdio>

#include<iostream>

#include<cstring>

#include<cstdlib>

#include<algorithm>

#include<cmath>

using namespace std;

typedef long long LL;

const int _=1e2;

const int maxn=1e3+_;

const int maxm=*1e5+_;

int read()

{

    int x=,f=; char ch=getchar();

    while(!isdigit(ch)){if(ch=='-')f=-;ch=getchar();}

    while(isdigit(ch)){x=x*+ch-'';ch=getchar();}

    return x*f;

}

void write(int x)

{

    if(x>=)write(x/);

    putchar(x%+'');

}

struct Trie

{

    int w[],c;

}tr[maxm*];int trlen,rt[maxm];

int insert(int d)

{

    int now=++trlen,ret=now;

    tr[now].c++;

    for(int i=;i>=;i--)

    {

        int x=(d&(<<i))?:;

        tr[now].w[x]=++trlen;

        now=tr[now].w[x];

        tr[now].c++;

    }

    return ret;

}

int merge(int x,int y)

{

    if(x==||y==)return x+y;

    tr[x].c+=tr[y].c;

    tr[x].w[]=merge(tr[x].w[],tr[y].w[]);

    tr[x].w[]=merge(tr[x].w[],tr[y].w[]);

    return x;

}

int c[maxn],nx[maxn],ny[maxn],g[maxn];

int main()

{

    freopen("a.in","r",stdin);

    freopen("a.out","w",stdout);

    int n,m,d;

    n=read(),m=read();

    for(int i=;i<=n;i++)c[i]=read();

    for(int i=;i<=m;i++)

    {

        d=read();

        rt[i]=insert(d);

        rt[i]=merge(rt[i],rt[i-]);

    }

    int Q,k,x,y,u,v;

    Q=read();

    while(Q--)

    {

        x=read(),y=read(),u=read(),v=read(),k=read();

        int ans=;

        for(int i=x;i<=y;i++)nx[i]=rt[u-],ny[i]=rt[v];

        for(int i=;i>=;i--)

        {

            int num=;

            for(int p=x;p<=y;p++)//看看如果这一位变成1的方案数是否已经超过k了

            {

                g[p]=(c[p]&(<<i))?:;//走那边变成1

                num+=tr[tr[ny[p]].w[g[p]]].c-tr[tr[nx[p]].w[g[p]]].c;

            }

            if(num>=k)//超过了，动用1就行了

            {

                ans|=(<<i);

                for(int p=x;p<=y;p++)

                    nx[p]=tr[nx[p]].w[g[p]],ny[p]=tr[ny[p]].w[g[p]];

            }

            else //不够，要用0补

            {

                k-=num;

                for(int p=x;p<=y;p++)

                    nx[p]=tr[nx[p]].w[g[p]^],ny[p]=tr[ny[p]].w[g[p]^];

            }

        }

        write(ans),puts("");

    }

    return ;

}

bzoj4103: [Thu Summer Camp 2015]异或运算的更多相关文章

  1. [BZOJ4103][Thu Summer Camp 2015]异或运算 可持久化Trie树

    4103: [Thu Summer Camp 2015]异或运算 Time Limit: 20 Sec  Memory Limit: 512 MB Description 给定长度为n的数列X={x1 ...

  2. BZOJ4103 [Thu Summer Camp 2015]异或运算 【可持久化trie树】

    题目链接 BZOJ4103 题解 一眼看过去是二维结构,实则未然需要树套树之类的数据结构 区域异或和,就一定是可持久化\(trie\)树 观察数据,\(m\)非常大,而\(n\)和\(p\)比较小,甚 ...

  3. bzoj4103 [Thu Summer Camp 2015]异或运算（可持久化trie）

    内存限制:512 MiB 时间限制:1000 ms 题目描述 给定长度为n的数列X={x1,x2,...,xn}和长度为m的数列Y={y1,y2,...,ym},令矩阵A中第i行第j列的值Aij=xi ...

  4. 【BZOJ 4103】 4103: [Thu Summer Camp 2015]异或运算 （可持久化Trie）

    4103: [Thu Summer Camp 2015]异或运算 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 474  Solved: 258 De ...

  5. 【bzoj4103】[Thu Summer Camp 2015]异或运算 可持久化trie树

    Description 给定长度为n的数列X={x1,x2,...,xn}和长度为m的数列Y={y1,y2,...,ym},令矩阵A中第i行第j列的值Aij=xi xor yj,每次询问给定矩形区域i ...

  6. [BZOJ 4103] [Thu Summer Camp 2015] 异或运算 【可持久化Trie】

    题目链接:BZOJ - 4103 题目分析 THUSC滚粗之后一直没有写这道题,从来没写过可持久化Trie,发现其实和可持久化线段树都是一样的.嗯,有些东西就是明白得太晚. 首先Orz ZYF-ZYF ...

  7. 【BZOJ 4103】 [Thu Summer Camp 2015]异或运算 可持久化01Trie

    我们观察数据:树套树 PASS    主席树 PASS  一层一个Trie PASS 再看,异或!我们就把目光暂时定在01Tire然后我们发现,我们可以带着一堆点在01Trie上行走,因为O(n*q* ...

  8. BZOJ 4103: [Thu Summer Camp 2015]异或运算 可持久化trie

    开始想了一个二分+可持久化trie验证,比正解多一个 log 仔细思考,你发现你可以直接按位枚举,然后在可持久化 trie 上二分就好了. code: #include <bits/stdc++ ...

  9. 【BZOJ 4104】【Thu Summer Camp 2015】解密运算

    http://www.lydsy.com/JudgeOnline/problem.php?id=4104 网上题解满天飞,我也懒得写了 #include<cstdio> #include& ...

随机推荐

  1. Split The Tree

    Split The Tree 时间限制: 1 Sec  内存限制: 128 MB 题目描述 You are given a tree with n vertices, numbered from 1 ...

  2. java 汉字保存到mysql 乱码

    保存之前正常,插入数据乱码 确认jsp mysql编码都确定为utf8 在连接数据库是加上编码 jdbc:mysql://localhost:3306/test?useUnicode=true& ...

  3. P3102 [USACO14FEB]秘密代码Secret Code

    题目描述 Farmer John has secret message that he wants to hide from his cows; the message is a string of ...

  4. Bootstrap开启模态框后对数据处理(标记模态框的开启与关闭状态)

    JS用全局变量标记状态,方法中动态修改全局变量以标记状态是一个重要思想. 需求:组合条件查询数据,查询完之后填充到模态框中,开启模态框,模态框中有组合条件查询,此时查询只需要更新模态框表格数据不需要开 ...

  5. JS设置页面中方法执行一次的思想

    思想:在JS中定义一全局变量,在方法执行的时候根据全局变量的值判断是否需要执行,在方法中修改全局变量的值,可以使得方法只执行一次.: 例如: 定义全局变量: var isLoad = false;// ...

  6. SQL行转列 （及EAV模型获取数据）

    参考文章: http://www.williamsang.com/archives/1508.html 情景简介 学校里面记录成绩,每个人的选课不一样,而且以后会添加课程,所以不需要把所有课程当作列. ...

  7. nginx原配置

    #原配置 server { listen ; server_name localhost; #charset koi8-r; #access_log logs/host.access.log main ...

  8. MD5进行文件完整性校验的操作方法

    我组产品包含大量音频和图片资源,MD5主要就用来检测这些资源文件的完整性.主要思路是:先计算出所有资源文件的MD5值,存到一个xml文件中,作为标准的MD5值.然后把这个xml文件放到我们的产品中,每 ...

  9. 【转载】5种网络IO模型

    同步(synchronous) IO和异步(asynchronous) IO,阻塞(blocking) IO和非阻塞(non-blocking)IO分别是什么,到底有什么区别?这个问题其实不同的人给出 ...

  10. iOS开发核心语言Objective C —— 全部知识点总结

    本分享是面向有意向从事iOS开发的伙伴及苹果产品的发烧友,亦或是已经从事了iOS的开发人员,想进一步提升者.假设您对iOS开发有极高的兴趣,能够与我一起探讨iOS开发.一起学习,共同进步.假设您是零基 ...