【BZOJ】2741: 【FOTILE模拟赛】L
题意:给定一个长度为n的序列,m次询问,每次询问一个区间[l, r],求max(Ai xor Ai+1 xor Ai+2 ... xor Aj),其中l<=i<=j<=r。(n<=12000, m<=6000, Ai在signed longint范围内)
#include <bits/stdc++.h>
using namespace std; const int nT=5000005, N=12005;
struct node *null;
struct node {
node *c[2];
int s;
}pool[nT], *it=pool;
node *newT() { node *x; x=it++; x->c[0]=x->c[1]=null; x->s=0; return x; }
node *add(node *p, int y, int dep=30) {
node *x=newT();
if(dep==-1) return x;
*x=*p;
int f=(y>>dep)&1;
x->c[f]=add(p->c[f], y, dep-1); ++x->c[f]->s;
return x;
}
int query(node *a, node *b, int k, int dep=30) {
if(dep==-1) return 0;
int f=(k>>dep)&1, ff=f^1;
if(b->c[ff]->s>a->c[ff]->s) return (1<<dep)|query(a->c[ff], b->c[ff], k, dep-1);
return query(a->c[f], b->c[f], k, dep-1);
}
node *Root[N], **root=Root+1;
int a[N], n, m, size, num, pos[N], d[120][N], last; void init() {
null=new node; null->c[0]=null->c[1]=null; null->s=0;
root[0]=null;
root[-1]=null; scanf("%d%d", &n, &m); //++n;
for(int i=1; i<=n; ++i) scanf("%d", &a[i]), a[i]^=a[i-1];
for(int i=1; i<=n; ++i) root[i]=add(root[i-1], a[i]);
size=sqrt(n+0.5);
num=(n-1)/size+1;
for(int i=1; i<=n; ++i) pos[i]=(i-1)/size+1;
for(int i=1; i<=num; ++i) {
int pi=(i-1)*size+1;
for(int j=pi; j<=n; ++j) d[i][j]=max(d[i][j-1], query(root[pi-1], root[j], a[j]));
}
}
int getans(int l, int r) {
--l;
int bl=pos[l], ret=0;
if(pos[r]>bl) ret=d[bl+1][r];
for(int i=l, end=min(r, bl*size); i<=end; ++i) ret=max(ret, query(root[l-1], root[r], a[i]));
return ret;
}
void solve() {
int x, y, l, r;
while(m--) {
scanf("%d%d", &x, &y);
x=(((long long)last+x)%n)+1;
y=(((long long)last+y)%n)+1;
l=min(x, y); r=max(x, y); //printf("%d %d\n", l, r);
printf("%d\n", last=getans(l, r));
}
}
int main() {
init();
solve();
return 0;
}
(自己的傻×做法,线段树套可持久化trie,mle成翔)
#include <bits/stdc++.h>
using namespace std; typedef long long ll;
const int M=8000005;
struct node *null;
struct node {
node *l, *r;
int cnt;
}pool[M], *bin[M], *it=pool, *nod;
int top;
node *newnode() {
node *x;
if(!top) { if(pool+M==it) exit(0); x=it++; }
else x=bin[--top];
x->l=x->r=null; x->cnt=1;
return x;
}
void clean(node *&x) {
if(x!=null && x!=nod && !(--x->cnt)) { bin[top++]=x; clean(x->l); clean(x->r); }
x=null;
}
node *merge(node *a, node *b) {
if(a==null) { ++b->cnt; return b; }
if(b==null) { ++a->cnt; return a; }
node *x=newnode();
x->l=merge(a->l, b->l);
x->r=merge(a->r, b->r);
return x;
}
node *add(node *p, ll y, int dep) {
if(dep==-1) return nod;
node *x=newnode();
x->l=p->l;
x->r=p->r;
if((y>>dep)&1) ++x->l->cnt, x->r=add(p->r, y, dep-1);
else ++x->r->cnt, x->l=add(p->l, y, dep-1);
return x;
}
node *add(node *p, ll y) { node *x=add(p, y, 63); clean(p); return x; }
node *build(node *a, node *b) {
if(a==null || b==null) return null;
node *x=newnode();
if(b->l!=null) {
x->l=build(a->l, b->l);
x->r=build(a->r, b->l);
}
if(b->r!=null) {
node *t;
x->l=merge(t=x->l, build(a->r, b->r)); clean(t);
x->r=merge(t=x->r, build(a->l, b->r)); clean(t);
}
return x;
}
void P(node *x, ll now, int dep=63) {
if(x==null) return;
if(dep==-1) printf("now:%lld\n", now);
P(x->l, now, dep-1);
P(x->r, now|(1<<dep), dep-1);
}
void D(node *x) {
P(x, 0); puts("");
} struct T {
node *l, *r, *all;
ll sum;
void pushup(T &lc, T &rc) {
sum=lc.sum^rc.sum;
static node *t;
l=merge(lc.l, build(rc.l, add(null, lc.sum)));
r=merge(rc.r, build(lc.r, add(null, rc.sum)));
all=merge(lc.all, rc.all); // printf("top:%d\n", top);
all=merge(t=all, build(lc.r, rc.l)); clean(t); // printf("sum:%lld , have:\n", sum); D(all);
}
ll get(node *x, int dep, bool flag) {
if(dep==-1) return 0;
ll ret=0;
if(x->l!=null && x->r!=null)
if(flag) ret=get(x->l, dep-1, 0);
else ret=(1ll<<dep), ret|=get(x->r, dep-1, 0);
else if(x->l==null && x->r!=null)
ret=(1ll<<dep), ret|=get(x->r, dep-1, 0);
else if(x->r==null && x->l!=null)
ret=get(x->l, dep-1, 0);
else puts("error");
return ret;
}
ll get() { return get(all, 63, 1); }
void clr() { clean(l); clean(r); clean(all); }
T& operator=(const T &a) { l=a.l; r=a.r; all=a.r; sum=a.sum; ++l->cnt; ++r->cnt; ++all->cnt; return *this; }
}t[12005<<2]; int n, q;
ll a[12005];
void build(int l, int r, int x) {
if(l==r) {
t[x].l=t[x].r=t[x].all=add(null, a[l]); // printf("%lld\n", a[l]); D(t[x].all);
t[x].sum=a[l];
return;
}
int mid=(l+r)>>1, lc=x<<1, rc=lc|1;
build(l, mid, lc);
build(mid+1, r, rc);
t[x].pushup(t[lc], t[rc]);
}
void query(int l, int r, int x, int L, int R, T &ret) {
if(L<=l && r<=R) { ret=t[x]; return; }
int mid=(l+r)>>1;
if(R<=mid) { query(l, mid, x<<1, L, R, ret); return; }
else if(mid<L) { query(mid+1, r, x<<1|1, L, R, ret); return; }
T lc, rc;
query(l, mid, x<<1, L, R, lc);
query(mid+1, r, x<<1|1, L, R, rc);
ret.pushup(lc, rc); lc.clr(); rc.clr();
}
void init() {
null=new node; null->l=null->r=null; null->cnt=1;
nod=newnode();
}
int main() {
scanf("%d%d", &n, &q);
init();
for(int i=1; i<=n; ++i) scanf("%lld", &a[i]);
build(1, n, 1);
int last=0, x, y, l, r;
T t;
while(q--) {
scanf("%d%d", &x, &y);
x=((x+last)%n)+1;
y=((y+last)%n)+1;
l=min(x, y);
r=max(x, y); // printf("l:%d, r:%d\n", l, r);
query(1, n, 1, l, r, t);
printf("%d\n", last=t.get());
t.clr();
}
return 0;
}
题解:
我看到大家的tag是可持久化trie后..我就往这个方向思考了下...就yy出了第一种sb做法.......................
即:发现我们只需要维护二进制位...然后查询就是在对应区间一直向右走即可(特判符号位...),用线段树维护区间........可持久化合并trie.......可是你会发现.........有O(nlogn)次合并..每次合并O(size(trie)).................然后可能又tle又mle。反正我开了引用计数的垃圾回收也跪了...............玛雅,rewrite的节奏啊...3k啊.............然后查题解..................发现是分块= =....什么鬼................
吐槽:最近到底怎么了.....................越是放假就越颓废??我感觉一放假,什么事情都会发生= =...这几天天天家里来客人妈呀...........能不能愉快的做题了........而且由于脑袋不知怎么得..转的特别慢...脑洞到是还好..可是一些简单的东西我要想好久啊啊>_<。。。还有本题竟然没有负数.....而且在int范围内............
ydc狂D蒟蒻我TAT....我已经被D成狗了.........ydc(先发了我这篇还没写完的博文....又发了一个自己讲莫队分块的博文的链接...下边是有我20天前的脑残留言的):“我怎么感到是在羞辱我”...............................(妈呀...还是在半群里...
真正的题解:
首先转化问题(妈呀完全没想到啊...),由于xor操作有特殊的性质,即$a \otimes b \otimes b = a$所以我们要求一段区间[a,b]的xor和,可以由[0, b]和[0, a-1]的xor值得来,设$s[i]=a[1] \otimes a[2] \otimes \cdots \otimes a[i]$,那么[l, r]区间的xor和就是 $s[r] \otimes s[l-1]$
果然我是脑残吗...
然后考虑区间最大的....即我们要求区间$[l, r]$最大的xor值,就是要求
$$max\{ s[j-1] \otimes s[i], l<=j<=i<=r \}$$
转化一下,其实就是要求$[l-1, r]$这个区间取两个数$i, j$求最大的$s[j] \otimes s[i]$(可以直接忽略i=j的情况,对结果无影响...)
容易得到dp方程,设d[i, j]表示区间$[i, j]$取两个数最大的xor值
$$d[i,j]=max\{d[i,j-1], f(s[j], \{s[l-1 \cdots r] \} ) \}, f表示将s[j]放到那个集合中找到一个最大的xor和\}$$
而这个方程外层都已经是$O(n^2)$的...
先考虑决策....容易想到这是个类似trie的东西....然后自己就能yy出来可持久化trie了...trie的边就代表一个二进制位..然后从上往下走即可....请自行yy...(....可持久化用在两个地方...1、省内存...2、区间信息...
然后就是分块大法!
引用ydc:“一般莫队能做的,分块+可持久化线段树/可持久化块状链表都是能做的”,orz
我们分块处理上边的方程..先分成$\sqrt n$块..然后方程$d[i,j]$表示第$i$块的起点到第$j$个元素的最大异或值...方程转移一样...这样预处理就是$O(31*n^{1.5})$了...(31是位数别说你看不懂..即log(2^31)
然后对于连续的我们直接得到答案,最后在处理一下非连续的即可...
总复杂度
$O(31*n^{1.5})$预处理,$O(31*mn^{0.5})$查询...可以水过..
【BZOJ】2741: 【FOTILE模拟赛】L的更多相关文章
- BZOJ.2741.[FOTILE模拟赛]L(分块 可持久化Trie)
题目链接 首先记\(sum\)为前缀异或和,那么区间\(s[l,r]=sum[l-1]^{\wedge}sum[r]\).即一个区间异或和可以转为求两个数的异或和. 那么对\([l,r]\)的询问即求 ...
- bzoj 2741 [FOTILE模拟赛] L
Description 多个询问l,r,求所有子区间异或和中最大是多少 强制在线 Solution 分块+可持久化trie 1.对于每块的左端点L,预处理出L到任意一个i,[L,j] 间所有子区间异或 ...
- 【BZOJ2741】【块状链表+可持久化trie】FOTILE模拟赛L
Description FOTILE得到了一个长为N的序列A,为了拯救地球,他希望知道某些区间内的最大的连续XOR和. 即对于一个询问,你需要求出max(Ai xor Ai+1 xor Ai+2 .. ...
- BZOJ2741 FOTILE模拟赛L(分块+可持久化trie)
显然做个前缀和之后变成询问区间内两个数异或最大值. 一种暴力做法是建好可持久化trie后直接枚举其中一个数查询,复杂度O(nmlogv). 观察到数据范围很微妙.考虑瞎分块. 设f[i][j]为第i个 ...
- 【bzoj2741】[FOTILE模拟赛]L 可持久化Trie树+分块
题目描述 FOTILE得到了一个长为N的序列A,为了拯救地球,他希望知道某些区间内的最大的连续XOR和. 即对于一个询问,你需要求出max(Ai xor Ai+1 xor Ai+2 ... xor A ...
- BZOJ2741:[FOTILE模拟赛]L
Description FOTILE得到了一个长为N的序列A,为了拯救地球,他希望知道某些区间内的最大的连续XOR和. 即对于一个询问,你需要求出max(Ai xor Ai+1 xor Ai+2 .. ...
- 【bzoj2741】[FOTILE模拟赛] L
Portal --> bzoj2741 Solution 突然沉迷分块不能自拔 考虑用分块+可持久化trie来解决这个问题 对于每一块的块头\(L\),预处理\([L,i]\)区间内的所有子区间 ...
- 【BZOJ】【2741】【FOTILE模拟赛】L
可持久化Trie+分块 神题……Orz zyf & lyd 首先我们先将整个序列搞个前缀异或和,那么某一段的异或和,就变成了两个数的异或和,所以我们就将询问[某个区间中最大的区间异或和]改变成 ...
- bzoj 2741: 【FOTILE模拟赛】L 分塊+可持久化trie
2741: [FOTILE模拟赛]L Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 1116 Solved: 292[Submit][Status] ...
- BZOJ2741: 【FOTILE模拟赛】L
2741: [FOTILE模拟赛]L Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 1170 Solved: 303[Submit][Status] ...
随机推荐
- linux中创建gpio节点
转自:http://blog.chinaunix.net/uid-29165999-id-4296162.html #define GPIO_MAJOR 230 // major device NO. ...
- C# 面试宝典
1.简述 private. protected. public. internal 修饰符的访问权限. private 私有成员 只有类成员才能访问 protected 保护成员 只有该类及该类的 ...
- MVC4 使用概述
1.Bundle使用: http://www.cnblogs.com/inline/p/3897256.html 2.MVC总结: http://www.cnblogs.com/xlhblogs/ar ...
- MVC4 WEBAPI(一)使用概述
所谓概述,也就是总结一些WEB API常用的使用用法.MVC APIWEB是一个轻量级的服务接口,完全符合RestFul框架设计,每个URL代表一种资源,使用方便,没有WCF那么庞大,但是麻雀虽小五脏 ...
- 库函数系统调用文件方式,王明学learn
库函数系统调用文件方式 基于C函数库的文件编程是独立于具体的操作系统平台的,不管是在Windows.Linux还是其他的操作系统中,都是使用这些函数.使用库函数进行程序设计可提高程序的可移植性. 对于 ...
- C字符数组赋值(转)
举例如下: char a[10];1.定义的时候直接用字符串赋值char a[10]="hello";注意:不能先定义再给它赋值,如 char a[10]; a[10]=" ...
- ios 多文件上传
/** * 上传多个文件 * * @param url 请求接口地址 * @param filedata 文件名称和数据(key:value) * @param btnName 上 ...
- Arduino101学习笔记(十四)—— Flash库
一.一些API 1.打开文件 SerialFlashFile file; file = SerialFlash.open("filename.bin"); if (file) { ...
- 软引用SoftReference异步加载图片
HashMap<String, SoftReference<Drawable>> imageCache 关于SoftReference这个类多少知道些机制,会用就ok了. 机制 ...
- SQL Server2008函数大全(完整版)
SQL2008 表达式:是常量.变量.列或函数等与运算符的任意组合. 1. 字符串函数 函数 名称 参数 示例 说明 ascii(字符串表达式) select ascii('abc') 返回 97 返 ...