博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
SDOI2011 染色
阅读量:5298 次
发布时间:2019-06-14

本文共 4201 字,大约阅读时间需要 14 分钟。

这个题还是树剖的基本练习,不过线段树维护的东西复杂了一些。

还是树剖起手,之后我们要在线段树中维护当前区间内颜色个数,最靠左的颜色和最靠右的颜色。这样的话我们在合并的时候,如果左区间最靠右颜色和右区间最靠左颜色是相同的话,那么这个合并之后的区间的颜色个数要-1(因为原来计算的方法是左右两区间之和)。

然后在修改和查询的时候也一样,修改没啥好说的,查询的话因为我写的是三段式查询,在最后一种情况的时候,因为答案也需要合并所以需要像上面所说的一样判断一次。至于pushdown的话,这个无需担心,因为pushdown一定是完全包括的,所以直接把下一级区间的颜色个数赋为1,左右颜色为当前区间颜色即可。

最后要注意的就是在沿着链上跳的时候,两个链的相交处需要判断一下颜色是否相同,如果相同答案-1.

然后这个题我就G了,原因是我线段树传参传错了……

看一下代码。

#include
#include
#include
#include
#include
#include
#include
#define rep(i,a,n) for(int i = a;i <= n;i++)#define per(i,n,a) for(int i = n;i >= a;i--)#define enter putchar('\n')using namespace std;typedef long long ll;const int M = 200005;int read(){ int ans = 0,op = 1; char ch = getchar(); while(ch < '0' || ch > '9') { if(ch == '-') op = -1; ch = getchar(); } while(ch >= '0' && ch <= '9') { ans *= 10; ans += ch - '0'; ch = getchar(); } return ans * op;}struct seg{ int v,lazy,lcol,rcol;}t[M<<2];struct edge{ int next,to,from;}e[M<<1];int n,m,dfn[M],rk[M],a[M],size[M],dep[M],fa[M],hson[M],q,x,y,head[M],ecnt,idx,top[M],z;char s[5];void add(int x,int y){ e[++ecnt].to = y; e[ecnt].next = head[x]; head[x] = ecnt;}void dfs1(int x,int f,int depth){ size[x] = 1,fa[x] = f,dep[x] = depth; int maxson = -1; for(int i = head[x];i;i = e[i].next) { if(e[i].to == f) continue; dfs1(e[i].to,x,depth+1); size[x] += size[e[i].to]; if(size[e[i].to] > maxson) maxson = size[e[i].to],hson[x] = e[i].to; }}void dfs2(int x,int t){ dfn[x] = ++idx,rk[idx] = a[x],top[x] = t; if(!hson[x]) return; dfs2(hson[x],t); for(int i = head[x];i;i = e[i].next) { if(e[i].to == fa[x] || e[i].to == hson[x]) continue; dfs2(e[i].to,e[i].to); }}void update(int p){ t[p].v = t[p<<1].v + t[p<<1|1].v; if(t[p<<1].rcol == t[p<<1|1].lcol) t[p].v--; t[p].lcol = t[p<<1].lcol,t[p].rcol = t[p<<1|1].rcol;}void build(int p,int l,int r){ if(l == r) { t[p].v = 1; t[p].lcol = t[p].rcol = rk[l]; return; } int mid = (l+r) >> 1; build(p<<1,l,mid),build(p<<1|1,mid+1,r); update(p);}void pushdown(int p,int l,int r){ if( l == r ) while(1); t[p<<1].v = t[p<<1|1].v = 1; t[p<<1].rcol = t[p<<1].lcol = t[p<<1|1].rcol = t[p<<1|1].lcol = t[p].lazy; t[p<<1].lazy = t[p<<1|1].lazy = t[p].lazy; t[p].lazy = 0;}int acol(int p,int l,int r,int pos){ if(l == r) return t[p].lcol; int mid = (l+r) >> 1; if(t[p].lazy) pushdown(p,l,r); if(pos <= mid) return acol(p<<1,l,mid,pos); else return acol(p<<1|1,mid+1,r,pos);}int query(int p,int l,int r,int kl,int kr){ int ans = 0,pcol,qcol; if(l == kl && r == kr) return t[p].v; int mid = (l+r) >> 1; if(t[p].lazy) pushdown(p,l,r); if(kr <= mid) ans += query(p<<1,l,mid,kl,kr); else if(kl > mid) ans += query(p<<1|1,mid+1,r,kl,kr); else { ans += query(p<<1,l,mid,kl,mid),ans += query(p<<1|1,mid+1,r,mid+1,kr); pcol = acol(1,1,n,mid),qcol = acol(1,1,n,mid+1); if(pcol == qcol) ans--; } return ans;}void modify(int p,int l,int r,int kl,int kr,int col){ if(l == kl && r == kr) { t[p].v = 1; t[p].lcol = t[p].rcol = t[p].lazy = col; return; } int mid = (l+r) >> 1; if(t[p].lazy) pushdown(p,l,r); if(kr <= mid) modify(p<<1,l,mid,kl,kr,col); else if(kl > mid) modify(p<<1|1,mid+1,r,kl,kr,col); else modify(p<<1,l,mid,kl,mid,col),modify(p<<1|1,mid+1,r,mid+1,kr,col); update(p);}int qrange(int x,int y){ int ans = 0,pcol,qcol; while(top[x] != top[y]) { if(dep[top[x]] < dep[top[y]]) swap(x,y); ans += query(1,1,n,dfn[top[x]],dfn[x]); pcol = acol(1,1,n,dfn[top[x]]),qcol = acol(1,1,n,dfn[fa[top[x]]]); if(pcol == qcol) ans--; x = fa[top[x]]; } if(dep[x] > dep[y]) swap(x,y); ans += query(1,1,n,dfn[x],dfn[y]); return ans;}void mrange(int x,int y,int z){ while(top[x] != top[y]) { if(dep[top[x]] < dep[top[y]]) swap(x,y); modify(1,1,n,dfn[top[x]],dfn[x],z); x = fa[top[x]]; } if(dep[x] > dep[y]) swap(x,y); modify(1,1,n,dfn[x],dfn[y],z);}int main(){ n = read(),m = read(); rep(i,1,n) a[i] = read(); rep(i,1,n-1) x = read(),y = read(),add(x,y),add(y,x); dfs1(1,0,1),dfs2(1,1),build(1,1,n); while(m--) { scanf("%s",s); if(s[0] == 'Q') { x = read(),y = read(); printf("%d\n",qrange(x,y)); } else if(s[0] == 'C') { x = read(),y = read(),z = read(); mrange(x,y,z); } } return 0;}

 

转载于:https://www.cnblogs.com/captain1/p/9710932.html

你可能感兴趣的文章
ASP.NET 会话状态【一】
查看>>
day 38 进程 与 子进程 (开启进程 + 进程中的数据隔离 + 守护进程 + 进程中的其他方法 )...
查看>>
dubbo注册中心介绍
查看>>
在netty3.x中存在两种线程:boss线程和worker线程。
查看>>
1489 蜥蜴和地下室
查看>>
【HTML/XML 10】XML文档中的Schema文件
查看>>
C语言中整型运算取Ceiling问题
查看>>
分治——最近点对
查看>>
Spring中jdbc Template使用
查看>>
什么是你的核心竞争力?
查看>>
Network Monitoring in Software-Defined Networking :A Review(综述)
查看>>
libtiff4.04
查看>>
CentOS 7 下使用 Firewall
查看>>
软件设计原则与模式
查看>>
[微信小程序] 认识微信小程序及开发环境搭建
查看>>
json模块,pickle模块,shelve模块
查看>>
django 分页出现 UnorderedObjectListWarning 错误
查看>>
140711
查看>>
芹菜炒香菇黑木耳
查看>>
[Ynoi2016]谁的梦
查看>>