题目：https://www.lydsy.com/JudgeOnline/problem.php?id=4137
关于可持久化01trie树：https://www.cnblogs.com/LadyLex/p/7281110.html
看了看它的两道例题，就没写。
特殊商品可以直接用可持久化trie做。
其他部分用线段树分治。修改是单点的，询问是区间，原来想的是把询问区间定位后有 mlogn 个，在线段树的每个叶子上贡献一番；结果TLE了，因为若是在叶子处贡献，一个询问就要做 r-l+1 次。
要在线段树的每个节点上建可持久化trie，可持久化的顺序就是商店的编号；这样就能对于 mlogn 个区间，先二分找到在可持久化trie的哪一段有贡献，然后贡献给那个询问了。
做出每个节点的可持久化trie，不用trie树合并什么的，只要把当前节点表示的区间的修改排序然后一个一个加进去就行了。一共只会加 nlogn 次。
为了省时间，不用再每个节点的时候给那些修改按商店排序，而可以先按商店排好序，然后像 CDQ 分治那样稳定地把属于左边的给左边、属于右边的给右边。
在线段树的每个节点的时候，虽然可以认为同一个商店的在一棵trie树上，不同商店之间的trie树合并起来，但写的时候不能像 “if( cr==pr ) { ... sm[cr]++; ... }” 这样 “当在同一个商店的时候按普通构建trie树那样赋值” ，因为已经可持久化了， sm[cr]++ 的这个 cr 可能可以通过之前的 rt 走过来，就会出错。
注意空间不是 18*n 而是 19*n ，因为边有 18 层，节点就有 19 层。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#define ls Ls[cr]
#define rs Rs[cr]
#define pb push_back
using namespace std;
int rdn()
{
int ret=0;bool fx=1;char ch=getchar();
while(ch>‘9‘||ch<‘0‘){if(ch==‘-‘)fx=0;ch=getchar();}
while(ch>=‘0‘&&ch<=‘9‘)ret=ret*10+ch-‘0‘,ch=getchar();
return fx?ret:-ret;
}
int Mx(int a,int b){return a>b?a:b;}
const int N=1e5+5,K=17,M=19e5+5;//19 not 18 for -1
int n,m,ct,sp[N],st[N],en[N],ans[N];
int tot,Ls[N<<1],Rs[N<<1],len;
struct Node{
int x,id,t;
bool operator< (const Node &b)const
{return id<b.id;}
}vl[N],tp[N];
struct Ques{int l,r,x;}q[N];
vector<int> vt[N<<1];
namespace Tri{
int tot,c[M][2],sm[M],rt[N],bin[K+5];
void Init()
{bin[0]=1;for(int i=1;i<=K;i++)bin[i]=bin[i-1]<<1;}
void ins(int x,int cr,int pr)
{
rt[cr]=++tot; cr=rt[cr]; pr=rt[pr];
for(int t=K;t>=0;t--)
{
bool d=(x&bin[t]);
c[cr][d]=++tot; c[cr][!d]=c[pr][!d];
cr=c[cr][d]; pr=c[pr][d]; sm[cr]=sm[pr]+1;
}
}
int qry(int l,int r,int x)
{
l=rt[l-1]; r=rt[r]; int ret=0;
for(int t=K;t>=0;t--)
{
bool d=(x&bin[t]);
if(sm[c[r][!d]]-sm[c[l][!d]])
l=c[l][!d], r=c[r][!d], ret|=bin[t];
else l=c[l][d], r=c[r][d];
}
return ret;
}
}
void build(int l,int r,int cr)
{
if(l==r)return; int mid=l+r>>1;
ls=++tot; build(l,mid,ls);
rs=++tot; build(mid+1,r,rs);
}
void ins(int l,int r,int cr,int L,int R,int k)
{
if(l>=L&&r<=R){vt[cr].pb(k);return;}
int mid=l+r>>1;
if(L<=mid)ins(l,mid,ls,L,R,k);
if(mid<R)ins(mid+1,r,rs,L,R,k);
}
int fnd(int x,bool fx)
{
int l=1,r=len,ret=0;
while(l<=r)
{
int mid=l+r>>1;
if(!fx)
{if(tp[mid].id>=x)ret=mid,r=mid-1;else l=mid+1;}
else
{if(tp[mid].id<=x)ret=mid,l=mid+1;else r=mid-1;}
}
return ret;
}
void solve(int l,int r,int cr,int L,int R)
{
for(int i=L,j=1;i<=R;i++,j++)tp[j]=vl[i];//out of if()
if(vt[cr].size())
{
len=R-L+1; Tri::tot=0;//
for(int i=1;i<=len;i++)Tri::ins(tp[i].x,i,i-1);
for(int i=0,j=vt[cr].size();i<j;i++)
{
int bh=vt[cr][i]; Ques k=q[bh];
int tl=fnd(k.l,0), tr=fnd(k.r,1);//tl,tr not l,r!!!
if(!tl||!tr||tl>tr)continue;
ans[bh]=Mx(ans[bh],Tri::qry(tl,tr,k.x));
}
}
if(l==r)return; int mid=l+r>>1,p0=L-1;
int tl=R-L+1;
for(int i=1;i<=tl;i++)if(tp[i].t<=mid)vl[++p0]=tp[i];
int p1=p0;
for(int i=1;i<=tl;i++)if(tp[i].t>mid)vl[++p0]=tp[i];
solve(l,mid,ls,L,p1); solve(mid+1,r,rs,p1+1,R);
}
int main()
{
n=rdn();int T=rdn();
for(int i=1;i<=n;i++)sp[i]=rdn();
for(int i=1,op,d;i<=T;i++)
{
op=rdn();
if(!op)
{
m++; vl[m].id=rdn();vl[m].x=rdn();vl[m].t=m;
}
else
{
ct++; q[ct].l=rdn();q[ct].r=rdn();q[ct].x=rdn();
d=rdn(); st[ct]=Mx(1,m-d+1);en[ct]=m;
}
}
Tri::Init();
if(m)//
{
tot=1;build(1,m,1);
for(int i=1;i<=ct;i++)
if(st[i]<=en[i])ins(1,m,1,st[i],en[i],i);//if
sort(vl+1,vl+m+1); solve(1,m,1,1,m);
}
Tri::tot=0;
for(int i=1;i<=n;i++)Tri::ins(sp[i],i,i-1);
for(int i=1;i<=ct;i++)
{
ans[i]=Mx(ans[i],Tri::qry(q[i].l,q[i].r,q[i].x));
printf("%d\n",ans[i]);
}
return 0;
}
原文地址：https://www.cnblogs.com/Narh/p/10376459.html
时间： 2024-10-03 14:33:03