bzoj2243树链剖分+区间合并

时间:2023-12-09 15:23:37

树链上区间合并的问题比区间修改要复杂,因为每一条重链在线段树上分布一般都是不连续的,所以在进行链上操作时要手动将其合并起来,维护两个端点值

处理时的方向问题:lca->u是一个方向,lca->v是另一个方向,到最后合并这两个放向时都看左端点即可

#include<cstring>
#include<string>
#include<iostream>
#include<queue>
#include<cstdio>
#include<algorithm>
#include<map>
#include<cstdlib>
#include<cmath>
#include<vector>
//#pragma comment(linker, "/STACK:1024000000,1024000000"); using namespace std; #define INF 0x3f3f3f3f
#define maxn 200005 int v[maxn],fir[maxn],nex[maxn],e_max;
int son[maxn],fa[maxn],color[maxn],col[maxn],pos[maxn],siz[maxn],deep[maxn],top[maxn];
int tot;
int n,m; struct node
{
int l,r;
int lc,rc;
int num,tag;
node()
{
lc=rc=-;
num=;
}
} t[*maxn]; void init_()
{
memset(son,-,sizeof son);
memset(siz,,sizeof siz);
memset(fir,-,sizeof fir);
e_max=;
tot=;
} void add_edge(int s,int t)
{
int e=e_max++;
v[e]=t;
nex[e]=fir[s];
fir[s]=e;
} void dfs1(int k,int pre,int d)
{
deep[k]=d;
siz[k]++;
fa[k]=pre;
for(int i=fir[k]; ~i; i=nex[i])
{
int e=v[i];
if(e!=pre)
{
dfs1(e,k,d+);
siz[k]+=siz[e];
if(son[k]==-||siz[son[k]]<siz[e]) son[k]=e;
}
}
} void dfs2(int k,int sp)
{
top[k]=sp;
pos[k]=tot++;
col[pos[k]]=color[k];
if(son[k]==-) return ;
dfs2(son[k],sp);
for(int i=fir[k]; ~i; i=nex[i])
{
int e=v[i];
if(e!=fa[k]&&e!=son[k])
{
dfs2(e,e);
}
}
} node Merge(node A,node B)
{
if(A.num==) return B;
if(B.num==) return A;
node temp;
temp.lc=A.lc;
temp.rc=B.rc;
if(A.rc==B.lc) temp.num=A.num+B.num-;
else temp.num=A.num+B.num;
return temp;
} inline void pushdown(int k)
{
if(t[k].tag==-) return ;
t[k<<].tag=t[k<<|].tag=t[k].tag;
t[k<<].num=t[k<<|].num=;
t[k<<].lc=t[k<<].rc=t[k<<|].lc=t[k<<|].rc=t[k].tag;
t[k].tag=-;
} inline void pushup(int k)
{
t[k].lc=t[k<<].lc;
t[k].rc=t[k<<|].rc;
t[k].num=t[k<<].num+t[k<<|].num;
if(t[k<<].rc==t[k<<|].lc) t[k].num--;
} void init(int l,int r,int k)
{
t[k].l=l;
t[k].r=r;
t[k].tag=-;
if(l==r)
{
t[k].num=;
t[k].lc=col[l];
t[k].rc=col[r];
return ;
}
int mid=l+r>>;
init(l,mid,k<<);
init(mid+,r,k<<|);
pushup(k);
} void update(int d,int l,int r,int k)
{
if(t[k].l==l&&t[k].r==r)
{
t[k].num=;
t[k].tag=d;
t[k].lc=t[k].rc=d;
return ;
}
pushdown(k);
int mid=t[k].l+t[k].r>>;
if(r<=mid) update(d,l,r,k<<);
else if(l>mid) update(d,l,r,k<<|);
else
{
update(d,l,mid,k<<);
update(d,mid+,r,k<<|);
}
pushup(k);
} node query(int l,int r,int k)
{
if(t[k].l==l&&t[k].r==r)
{
return t[k];
}
pushdown(k);
int mid=t[k].l+t[k].r>>;
if(r<=mid) return query(l,r,k<<);
else if(l>mid) return query(l,r,k<<|);
else return Merge(query(l,mid,k<<),query(mid+,r,k<<|));
} void Query(int s,int t)
{
node L,R;
int f1=top[s],f2=top[t];
while(f1!=f2)
{
if(deep[f1]<deep[f2]) swap(L,R),swap(f1,f2),swap(s,t);
L=Merge(query(pos[f1],pos[s],),L);
s=fa[f1];
f1=top[s];
}
if(deep[s]<deep[t]) swap(L,R),swap(s,t);
L=Merge(query(pos[t],pos[s],),L);
if(L.lc==R.lc) L.num=L.num+R.num-;
else L.num=L.num+R.num;
printf("%d\n",L.num);
} void Change(int s,int t,int c)
{
int f1=top[s],f2=top[t];
while(f1!=f2)
{
if(deep[f1]<deep[f2]) swap(f1,f2),swap(s,t);
update(c,pos[f1],pos[s],);
s=fa[f1];
f1=top[s];
}
if(deep[s]>deep[t]) swap(s,t);
update(c,pos[s],pos[t],);
} int main()
{
while(scanf("%d%d",&n,&m)!=EOF)
{
init_();
for(int i=; i<=n; i++)
{
scanf("%d",&color[i]);
}
for(int i=; i<n; i++)
{
int a,b;
scanf("%d%d",&a,&b);
add_edge(a,b);
add_edge(b,a);
}
dfs1(,-,);
dfs2(,);
init(,tot-,);
while(m--)
{
char s[];
scanf("%s",s);
if(s[]=='Q')
{
int l,r;
scanf("%d%d",&l,&r);
Query(l,r);
}
else
{
int a,b,c;
scanf("%d%d%d",&a,&b,&c);
Change(a,b,c);
}
}
}
return ;
}