[洛谷P4585] [FJOI2015] 火星商店问题

Description

火星上的一条商业街里按照商店的编号 \(1\)，\(2\) ，…，\(n\) ，依次排列着 \(n\) 个商店。商店里出售的琳琅满目的商品中，每种商品都用一个非负整数 \(val\) 来标价。每个商店每天都有可能进一些新商品，其标价可能与已有商品相同。

火星人在这条商业街购物时，通常会逛这条商业街某一段路上的所有商店，譬如说商店编号在区间 \([L,R]\) 中的商店，从中挑选1件自己最喜欢的商品。每个火星人对商品的喜好标准各不相同。通常每个火星人都有一个自己的喜好密码 \(x\) 。对每种标价为 \(val\) 的商品，喜好密码为 \(x\) 的火星人对这种商品的喜好程度与 \(val\) 异或 \(x\) 的值成正比。也就是说，\(val\) \(xor\) \(x\) 的值越大，他就越喜欢该商品。每个火星人的购物卡在所有商店中只能购买最近 \(d\) 天内（含当天）进货的商品。另外，每个商店都有一种特殊商品不受进货日期限制，每位火星人在任何时刻都可以选择该特殊商品。每个商店中每种商品都能保证供应，不存在商品缺货的问题。

对于给定的按时间顺序排列的事件，计算每个购物的火星人的在本次购物活动中最喜欢的商品，即输出 \(val\) \(xor\) \(x\) 的最大值。这里所说的按时间顺序排列的事件是指以下 \(2\) 种事件：

事件 \(0\)， 用三个整数 \(0\) , \(s\) , \(v\) ，表示编号为 \(s\) 的商店在当日新进一种标价为 \(v\) 的商品。

事件 \(1\)，用 \(5\) 个整数 \(1\) , $ L$ , \(R\) , \(x\) , \(d\) ，表示一位火星人当日在编号为 \(L\) 到 \(R\) 的商店购买 \(d\) 天内的商品，该火星人的喜好密码为 \(x\) 。

Input

第1行中给出2个正整数 \(n\) , \(m\) ，分别表示商店总数和事件总数。

第2行中有 \(n\) 个整数，第 \(i\) 个整数表示商店 \(i\) 的特殊商品标价。

接下来的 \(m\) 行，每行表示 \(1\) 个事件。每天的事件按照先事件 \(0\)， 后事件 \(1\) 的顺序排列。

Output

将计算出的每个事件 \(1\) 的 \(val\) \(xor\) \(x\) 的最大值依次输出。

Sample Input

4 6

1 2 3 4

1 1 4 1 0

0 1 4

0 1 3

1 1 1 1 0

1 1 1 1 1

1 1 2 1 2

Sample Output

5

0

2

5

HINT

\(n, m \leq 100000\)

数据中，价格不大于 \(100000\)

---

想法

首先吐槽题面！没见过比它说得更不清楚更恶心的了。。。

“每天”的定义：出现事件0 就是新一天的开始，在第一个事件0出现前的事件可看作第0天的。

要求异或和最大有两种方法，一种是线性基，一种是可持久化 \(trie\) ，这道题中可持久化 \(trie\) 更方便。

对于不受时间限制的特殊商品，直接以商店编号为版本号建一个 可持久化 \(trie\) ，查询、更新答案。

其他操作咋一看好像很复杂。

但如果商店编号当作 \(x\) 轴，时间当作 \(y\) 轴，其实就是 添加单点+二维矩阵查询与 \(x\) 的最大异或和。

那就需要树套树（线段树+可持久化 \(trie\)）了。

以商店编号为下标建线段树，在线段树每个节点以时间为版本号建 可持久化 \(trie\) （反过来也行啦~）

但如果这样空间受不住。

于是线段树分治。

将每个每个查询拆成 \(O(logn)\) 份放到线段树节点上

遍历线段树的每个节点，现场建树、查询、更新答案。

总复杂度 \(O(nlog^2n)\)

---

代码

这代码……不太好写也不太好调。

唯一的安慰是数据弱，调过样例后就 \(1A\) 了。

（\(ps\) 指针选手似乎变为数组选手了呢……）

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<vector>
#include<cstring>

using namespace std;

int read(){
	int x=0;
	char ch=getchar();
	while(!isdigit(ch)) ch=getchar();
	while(isdigit(ch)) x=x*10+ch-'0',ch=getchar();
	return x;
}

const int N = 100005;

struct commo{ int x,d; };
struct query{ int x,l,r,id; };
int n,m;
int spi[N],ans[N];

int rt[N],cnt2,son[N*20][2],sum[N*20];
void clear(int x){
	for(int i=0;i<cnt2;i++) son[i][0]=son[i][1]=0,sum[i]=0;
	cnt2=0;
	for(int i=0;i<=x;i++) rt[i]=0;
	if(x==-1) memset(rt,0,sizeof(rt));
}
void insert(int x,int pre,int dep,int c){
	if(dep==-1) { sum[x]=sum[pre]+1; return; }
	int f=((c>>dep)&1);
	son[x][f^1]=son[pre][f^1];
	insert(son[x][f]=++cnt2,son[pre][f],dep-1,c);
	sum[x]=sum[pre]+1;
}
int ask(int x,int pre,int dep,int c){
	if(dep==-1) return 0;
	int f=((c>>dep)&1);
	if(sum[son[x][f^1]]>sum[son[pre][f^1]]) return (1<<dep)+ask(son[x][f^1],son[pre][f^1],dep-1,c);
	return ask(son[x][f],son[pre][f],dep-1,c);
}

int root,cnt;
int ch[N*2][2];
vector<commo> v[N*2];
vector<query> q[N*2];
void build(int x,int l,int r){
	if(l==r) return;
	int mid=(l+r)>>1;
	build(ch[x][0]=++cnt,l,mid);
	build(ch[x][1]=++cnt,mid+1,r);
}
void modify(int x,int l,int r,int c,commo d){
	v[x].push_back(d);
	if(l==r) return;
	int mid=(l+r)>>1;
	if(c<=mid) modify(ch[x][0],l,mid,c,d);
	else modify(ch[x][1],mid+1,r,c,d);
}
void add(int x,int l,int r,int L,int R,query d){
	if(l>=L && r<=R) { q[x].push_back(d); return; }
	int mid=(l+r)>>1;
	if(R>mid) add(ch[x][1],mid+1,r,L,R,d);
	if(L<=mid) add(ch[x][0],l,mid,L,R,d);
}

int b[N],tot;
void work(int x,int l,int r){
	if(q[x].size()){
		tot=0;
		for(int i=0;i<v[x].size();i++) b[++tot]=v[x][i].d;
		sort(b+1,b+1+tot);
		tot=unique(b+1,b+1+tot)-b-1;
		for(int i=0;i<v[x].size();i++){
			int id=lower_bound(b+1,b+1+tot,v[x][i].d)-b,pre;
			pre=(rt[id] ? rt[id] : rt[id-1]);
			insert(rt[id]=++cnt2,pre,17,v[x][i].x);
		}

		for(int i=0;i<q[x].size();i++){
			if(q[x][i].l>q[x][i].r) continue;
			int ll=lower_bound(b+1,b+1+tot,q[x][i].l)-b-1,rr=upper_bound(b+1,b+1+tot,q[x][i].r)-b-1;
			if(ll<rr) {
				ans[q[x][i].id]=max(ans[q[x][i].id],ask(rt[rr],rt[ll],17,q[x][i].x));
			}
		}

		clear(tot);
	}

	if(l==r) return;
	int mid=(l+r)>>1;
	work(ch[x][0],l,mid); work(ch[x][1],mid+1,r);
}

int main()
{
	n=read(); m=read();
	for(int i=1;i<=n;i++) spi[i]=read();

	clear(-1);
	for(int i=1;i<=n;i++)
		insert(rt[i]=++cnt2,rt[i-1],17,spi[i]);

	int opt,day=1,ques=0,l,r,d,x;
	build(root=++cnt,1,n);
	for(int i=0;i<m;i++){
		opt=read();
		if(opt==0) {
			l=read(); x=read();
			modify(root,1,n,l,(commo){x,++day});
		}
		else{
			l=read(); r=read(); x=read(); d=read();
			add(root,1,n,l,r,(query){x,max(day-d+1,1),day,ques++});
			ans[ques-1]=ask(rt[r],rt[l-1],17,x);
		}
	}
	clear(n);

	work(root,1,n);
	for(int i=0;i<ques;i++) printf("%d\n",ans[i]);

	return 0;
}