bzoj4103: [Thu Summer Camp 2015]异或运算

对于每个询问暴力枚举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]异或运算的更多相关文章

[BZOJ4103][Thu Summer Camp 2015]异或运算可持久化Trie树
4103: [Thu Summer Camp 2015]异或运算 Time Limit: 20 Sec Memory Limit: 512 MB Description 给定长度为n的数列X={x1 ...
BZOJ4103 [Thu Summer Camp 2015]异或运算【可持久化trie树】
题目链接 BZOJ4103 题解一眼看过去是二维结构,实则未然需要树套树之类的数据结构区域异或和,就一定是可持久化$trie$树观察数据,$m$非常大,而$n$和$p$比较小,甚 ...
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 ...
【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 ...
【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 ...
[BZOJ 4103] [Thu Summer Camp 2015] 异或运算【可持久化Trie】
题目链接:BZOJ - 4103 题目分析 THUSC滚粗之后一直没有写这道题,从来没写过可持久化Trie,发现其实和可持久化线段树都是一样的.嗯,有些东西就是明白得太晚. 首先Orz ZYF-ZYF ...
【BZOJ 4103】 [Thu Summer Camp 2015]异或运算可持久化01Trie
我们观察数据:树套树 PASS 主席树 PASS 一层一个Trie PASS 再看,异或!我们就把目光暂时定在01Tire然后我们发现,我们可以带着一堆点在01Trie上行走,因为O(n*q* ...
BZOJ 4103: [Thu Summer Camp 2015]异或运算可持久化trie
开始想了一个二分+可持久化trie验证,比正解多一个 log 仔细思考,你发现你可以直接按位枚举,然后在可持久化 trie 上二分就好了. code: #include <bits/stdc++ ...
【BZOJ 4104】【Thu Summer Camp 2015】解密运算
http://www.lydsy.com/JudgeOnline/problem.php?id=4104 网上题解满天飞,我也懒得写了 #include<cstdio> #include& ...

随机推荐

回顾基础知识，类，fbv,cbv
一类中绑定方法的传参,不需要self class Foo(object): def __init__(self,name): self.name = name def foo(self,x): se ...
Linux（10）：期中架构（2）--- NFS存储服务 & 实时同步
1. 共享存储服务概念: # NFS是Network File System的缩写,中文意思是网络文件系统, # 它的主要功能是通过网络(一般是局域网)让不同的主机系统之间可以共享文件或目录. 2. ...
ElasticSearch 多索引
1.用逗号将索引隔开,如: $ curl -XPOST http://localhost:9200/aaa,website/_search/ { "took": 1, " ...
如何让Gridview在没有数据的时候显示表头[没有使用SqlDataSource控件时]
原文发布时间为:2008-08-03 -- 来源于本人的百度文章 [由搬家工具导入] 要看全文请点击http://blog.csdn.net/windok2004/archive/2007/10/28 ...
Mac快速查看隐藏文件
使用终端显示隐藏文件的最简单方法是使用终端.只要打开终端(位于应用程序--实用工具),将以下代码复制进去然后回车 defaults write com.apple.finder AppleShowA ...
POJ 1509 Glass Beads【字符串最小表示法】
题目链接: http://poj.org/problem?id=1509 题意: 求循环字符串的最小表示. 分析: 浅析"最小表示法"思想在字符串循环同构问题中的应用判断两字符串 ...
Cookie 和 Session 有什么区别呢？
Cookie 和 Session 有什么区别呢?大部分的面试者应该都可以说上一两句,比如:什么是 Cookie?什么是 Session?两者的区别等但如果再往深入探讨的话,就慢慢有一些朋友不太了解了 ...
LightOJ1234 Harmonic Number 调和级数求和
[题目] [预备知识] ,其中r是欧拉常数,const double r= 0.57721566490153286060651209; 这个等式在n很大的时候比较精确. [解法]可以在 n较小的时 ...
使用NSOperationQueue简化多线程开发
多线程开发是一件需要特别精心的事情,即使是对有多年开发经验的工程师来说. 为了能让初级开发工程师也能使用多线程,同时还要简化复杂性.各种编程工具提供了各自的办法.对于iOS来说,建议在尽可能的情况下避 ...
android 计时器
<?xml version="1.0" encoding="utf-8"?> <LinearLayout xmlns:android=&quo ...

bzoj4103: [Thu Summer Camp 2015]异或运算

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

随机推荐

热门专题