树链剖分+线段树【SDOI2011】 bzoj2243 染色

2024-04-14 23:58

本文主要是介绍树链剖分+线段树【SDOI2011】 bzoj2243 染色,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题目大意:

给一棵树,每个节点有一个颜色。写一个程序支持把两个点路径上的所有点染成一个颜色,查询两点之间的色段数量。


解题思路:

树链剖分+线段树

首先它是一颗树,而且是修改和查询两点路径上的颜色,可以想到树链剖分。

查询颜色段数可以用线段树维护区间颜色段数。

这道题涉及到区间合并,所以在线段树和lca的时候需要多记录一些东西,当前区间的最左边的颜色,最右边的颜色,已经求出的区间的左边的颜色。

我采用的做法是先求出lca,然后从左右两边各爬一边,就能有序的将所有的区间连在一起,当然大神告诉我可以左右同时爬,但是我不太会。


注意事项:

1、在树链剖分求lca的时候一定要写函数接口,不要把大段代码放在lca里,代码非常凌乱,不易差错;

2、打了标记要记得下传,时时刻刻思考要不要下传标记;

3、每次区间合并的时候要判断是否连接的部分是相同的颜色,如果相同要ans--;

4、从一侧爬的时候处理lca,另一侧就要避开lca,而且左右两条链合并的时候要判断颜色是否相同。

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<string>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 100100
using namespace std;
struct tree{int l,r,sum,ls,rs,lz;
}d[N*8];
int n,m,x,y,z,a,b;
char s[5];
int q[N];
int fir[N],nes[N*2],v[N*2];
bool pc[N*2];
int dep[N],fa[N],zon[N],sz[N],ld[N],xl[N],pos[N],top;
void dfs1(int c)
{dep[c]=dep[fa[c]]+1;sz[c]=1;for(int t=fir[c];t;t=nes[t]){if(pc[t]) continue;pc[t^1]=true;fa[v[t]]=c;dfs1(v[t]);sz[c]+=sz[v[t]];if(sz[v[t]]>sz[zon[c]]) zon[c]=v[t];}
}
void dfs2(int c)
{xl[++top]=c;pos[c]=top;ld[c]=c;if(zon[fa[c]]==c) ld[c]=ld[fa[c]];if(zon[c]) dfs2(zon[c]);for(int t=fir[c];t;t=nes[t]){if(pc[t] || v[t]==zon[c]) continue;dfs2(v[t]);}
}
void make_tree(int c,int l,int r)
{d[c].l=l;d[c].r=r;if(l==r){d[c].sum=1;d[c].ls=d[c].rs=q[xl[l]];d[c].lz=-1;return;}int mid=(l+r)/2;make_tree(c*2,l,mid);make_tree(c*2+1,mid+1,r);d[c].sum=d[c*2].sum+d[c*2+1].sum;if(d[c*2].rs==d[c*2+1].ls) d[c].sum--;d[c].ls=d[c*2].ls;d[c].rs=d[c*2+1].rs;d[c].lz=-1;
}
void push_down(int c)
{if(d[c].lz==-1) return;d[c].sum=1;d[c].ls=d[c].rs=d[c].lz;d[c*2].lz=d[c*2+1].lz=d[c].lz;d[c].lz=-1;return;
}
void change(int c)
{if(a<=d[c].l && b>=d[c].r) {d[c].lz=z;return;}push_down(c);int mid=(d[c].l+d[c].r)/2;if(a<=mid) change(c*2);if(b>mid)  change(c*2+1);push_down(c*2);push_down(c*2+1);d[c].sum=d[c*2].sum+d[c*2+1].sum;if(d[c*2].rs==d[c*2+1].ls) d[c].sum--;d[c].ls=d[c*2].ls;d[c].rs=d[c*2+1].rs;d[c].lz=-1;
}
void lca_change(int x,int y)
{while(ld[x]!=ld[y]){if(dep[ld[x]]<dep[ld[y]]) swap(x,y);a=pos[ld[x]]; b=pos[x];change(1);x=fa[ld[x]];}if(dep[x]<dep[y]) swap(x,y);a=pos[y]; b=pos[x];change(1);
}
int query(int c,int &zs,int &ys)
{push_down(c);if(a<=d[c].l && b>=d[c].r){zs=d[c].ls; ys=d[c].rs;return d[c].sum;}int mid=(d[c].l+d[c].r)/2;if(b<=mid) return query(c*2,zs,ys);if(a>mid)  return query(c*2+1,zs,ys);int zz,yy,ans=0;ans+=query(c*2,zs,yy);ans+=query(c*2+1,zz,ys);if(d[c*2].rs==d[c*2+1].ls) ans--;return ans;
}
int lca_query(int x,int y)
{int lx=x,ly=y;while(ld[lx]!=ld[ly]){if(dep[ld[lx]]<dep[ld[ly]]) swap(lx,ly);lx=fa[ld[lx]];}if(dep[lx]<dep[ly]) swap (lx,ly);int lca=ly,ans=0;int ls=-1,zs=-1,ys=-1,k1,k2;while(ld[x]!=ld[lca]){a=pos[ld[x]];b=pos[x];ans+=query(1,zs,ys);if(ls==ys) ans--;ls=zs;x=fa[ld[x]];}a=pos[lca];b=pos[x];ans+=query(1,zs,ys);if(ls==ys) ans--;k1=zs;ls=zs=ys=-1;while(ld[y]!=ld[lca]){a=pos[ld[y]];b=pos[y];ans+=query(1,zs,ys);if(ls==ys) ans--;ls=zs;y=fa[ld[y]];}if(y!=lca){a=pos[lca]+1; b=pos[y];ans+=query(1,zs,ys);if(ls==ys) ans--;}k2=zs;if(k1==k2) ans--;return ans;
}
int main()
{scanf("%d%d",&n,&m);for(int i=1;i<=n;i++) scanf("%d",&q[i]);for(int i=2;i<2*n;i+=2){scanf("%d%d",&x,&y);nes[i]=fir[x];fir[x]=i; v[i]=y;nes[i^1]=fir[y];fir[y]=i^1; v[i^1]=x;}dfs1(1);dfs2(1);make_tree(1,1,n);for(int i=1;i<=m;i++){scanf("%s",s);if(s[0]=='Q'){scanf("%d%d",&x,&y);printf("%d\n",lca_query(x,y));}else{scanf("%d%d%d",&x,&y,&z);lca_change(x,y);}}return 0;
}

这篇关于树链剖分+线段树【SDOI2011】 bzoj2243 染色的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/904388

相关文章

poj3468(线段树成段更新模板题)

题意:包括两个操作:1、将[a.b]上的数字加上v;2、查询区间[a,b]上的和 下面的介绍是下解题思路: 首先介绍  lazy-tag思想:用一个变量记录每一个线段树节点的变化值,当这部分线段的一致性被破坏我们就将这个变化值传递给子区间,大大增加了线段树的效率。 比如现在需要对[a,b]区间值进行加c操作,那么就从根节点[1,n]开始调用update函数进行操作,如果刚好执行到一个子节点,

hdu1394(线段树点更新的应用)

题意:求一个序列经过一定的操作得到的序列的最小逆序数 这题会用到逆序数的一个性质,在0到n-1这些数字组成的乱序排列,将第一个数字A移到最后一位,得到的逆序数为res-a+(n-a-1) 知道上面的知识点后,可以用暴力来解 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#in

hdu1689(线段树成段更新)

两种操作:1、set区间[a,b]上数字为v;2、查询[ 1 , n ]上的sum 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<queue>#include<set>#include<map>#include<stdio.h>#include<stdl

hdu 1754 I Hate It(线段树,单点更新,区间最值)

题意是求一个线段中的最大数。 线段树的模板题,试用了一下交大的模板。效率有点略低。 代码: #include <stdio.h>#include <string.h>#define TREE_SIZE (1 << (20))//const int TREE_SIZE = 200000 + 10;int max(int a, int b){return a > b ? a :

hdu 1166 敌兵布阵(树状数组 or 线段树)

题意是求一个线段的和,在线段上可以进行加减的修改。 树状数组的模板题。 代码: #include <stdio.h>#include <string.h>const int maxn = 50000 + 1;int c[maxn];int n;int lowbit(int x){return x & -x;}void add(int x, int num){while

poj 1127 线段相交的判定

题意: 有n根木棍,每根的端点坐标分别是 px, py, qx, qy。 判断每对木棍是否相连,当他们之间有公共点时,就认为他们相连。 并且通过相连的木棍相连的木棍也是相连的。 解析: 线段相交的判定。 首先,模板中的线段相交是不判端点的,所以要加一个端点在直线上的判定; 然后,端点在直线上的判定这个函数是不判定两个端点是同一个端点的情况的,所以要加是否端点相等的判断。 最后

HDU4737线段树

题目大意:给定一系列数,F(i,j)表示对从ai到aj连续求或运算,(i<=j)求F(i,j)<=m的总数。 const int Max_N = 100008 ;int sum[Max_N<<2] , x[Max_N] ;int n , m ;void push_up(int t){sum[t] = sum[t<<1] | sum[t<<1|1] ;}void upd

zoj 1721 判断2条线段(完全)相交

给出起点,终点,与一些障碍线段。 求起点到终点的最短路。 枚举2点的距离,然后最短路。 2点可达条件:没有线段与这2点所构成的线段(完全)相交。 const double eps = 1e-8 ;double add(double x , double y){if(fabs(x+y) < eps*(fabs(x) + fabs(y))) return 0 ;return x + y ;

圆与线段的交点

poj 3819  给出一条线段的两个端点,再给出n个圆,求出这条线段被所有圆覆盖的部分占了整条线段的百分比。 圆与线段的交点 : 向量AB 的参数方程  P = A + t * (B - A)      0<=t<=1 ; 将点带入圆的方程即可。  注意: 有交点 0 <= t <= 1 ; 此题求覆盖的部分。 则 若求得 t  满足 ; double ask(d

Codeforces 482B 线段树

求是否存在这样的n个数; m次操作,每次操作就是三个数 l ,r,val          a[l] & a[l+1] &......&a[r] = val 就是区间l---r上的与的值为val 。 也就是意味着区间[L , R] 每个数要执行 | val 操作  最后判断  a[l] & a[l+1] &......&a[r] 是否= val import ja