对于每个询问暴力枚举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. Spoj-FACVSPOW Factorial vs Power

    Consider two integer sequences f(n) = n! and g(n) = an, where n is a positive integer. For any integ ...

  2. 圆桌聚餐（cogs 729）

    «问题描述:假设有来自m 个不同单位的代表参加一次国际会议.每个单位的代表数分别为ri(i=1,2,3...m), .会议餐厅共有n张餐桌,每张餐桌可容纳c i(i=1,2...n) 个代表就餐.为了 ...

  3. 跳蚤 BZOJ 4310

    跳蚤 [问题描述] 很久很久以前,森林里住着一群跳蚤.一天,跳蚤国王得到了一个神秘的字符串,它想进行研究. 首先,他会把串分成不超过 k 个子串,然后对于每个子串 S,他会从S的所有子串中选择字典序最 ...

  4. Oracle11g服务详细介绍及哪些服务是必须开启的？

    中的方法成功安装Oracle 11g后,共有7个服务, 这七个服务的含义分别为: 1. Oracle ORCL VSS Writer Service: Oracle卷映射拷贝写入服务,VSS(Volu ...

  5. Python入门--2--继续学习

    继续学习小甲鱼 一.python比较操作符 == :判断左边是否等于右边 != : 判断左边是否不能右边 二. if while判断语句 栗子: temp = input ("sha shu ...

  6. SSH: Transferred 0 file(s) 解决

    Jenkins搭建过程中,使用 Publish Over SSH 插件.发生 SSH: Transferred 0 file(s). 百度.google了几个小时,终于找到答案,特此记录. 配置如下: ...

  7. grafana 安装 和 Nginx 、EL 联调

    https://blog.csdn.net/u010735147/article/details/80943593

  8. HDU 5676 ztr loves lucky numbers【DFS】

    题目链接; http://acm.hdu.edu.cn/showproblem.php?pid=5676 题意: 由4和7组成的且4和7出现次数相同的数称为幸运数字,给定n,求不大于n的最大幸运数字. ...

  9. Neo4j 第七篇：模式（Pattern）

    模式和模式匹配是Cypher的核心,使用模式来描述所需数据的形状,该模式使用属性图的结构来描述,通常使用小括号()表示节点,-->表示关系,-[]->表示关系和关系的类型,箭头表示关系的方 ...

  10. 创建ROS工作空间和包

    一.创建工作空间 mkdir -p ~/openni_ws/src cd ~/openni_ws catkin_make        //在catkin工作空间(openni_ws)下catkin_ ...