博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
Bzoj 3631: [JLOI2014]松鼠的新家(树链剖分+线段树)
阅读量:4935 次
发布时间:2019-06-11

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

3631: [JLOI2014]松鼠的新家

Time Limit: 10 Sec Memory Limit: 128 MB
Description
松鼠的新家是一棵树,前几天刚刚装修了新家,新家有n个房间,并且有n-1根树枝连接,每个房间都可以相互到达,且俩个房间之间的路线都是唯一的。天哪,他居然真的住在“树”上。松鼠想邀请小熊维尼前来参观,并且还指定一份参观指南,他希望维尼能够按照他的指南顺序,先去a1,再去a2,……,最后到an,去参观新家。
可是这样会导致维尼重复走很多房间,懒惰的维尼不听地推辞。可是松鼠告诉他,每走到一个房间,他就可以从房间拿一块糖果吃。维尼是个馋家伙,立马就答应了。
现在松鼠希望知道为了保证维尼有糖果吃,他需要在每一个房间各放至少多少个糖果。因为松鼠参观指南上的最后一个房间an是餐厅,餐厅里他准备了丰盛的大餐,所以当维尼在参观的最后到达餐厅时就不需要再拿糖果吃了。
Input
第一行一个整数n,表示房间个数
第二行n个整数,依次描述a1-an
接下来n-1行,每行两个整数x,y,表示标号x和y的两个房间之间有树枝相连。
Output
一共n行,第i行输出标号为i的房间至少需要放多少个糖果,才能让维尼有糖果吃。
Sample Input
5
1 4 5 3 2
1 2
2 4
2 3
4 5
Sample Output
1
2
1
2
1
HINT
2<= n <=300000

/*裸树剖.*/#include
#include
#define MAXN 300001using namespace std;int n,m,tot,cut,maxsize,head[MAXN],a[MAXN],s[MAXN];int size[MAXN],pos[MAXN],top[MAXN],fa[MAXN],deep[MAXN];struct data{
int l,r,lc,rc,sum,bj,size;}tree[MAXN<<2];struct edge{
int v,next;}e[MAXN<<1];int read(){ int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){
if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar(); return x*f;}void add(int u,int v){ e[++cut].v=v;e[cut].next=head[u];head[u]=cut;}void build(int l,int r){ int k=++tot; tree[k].l=l,tree[k].r=r,tree[k].size=r-l+1; if(l==r) return ; int mid=(l+r)>>1; tree[k].lc=tot+1;build(l,mid); tree[k].rc=tot+1;build(mid+1,r); return ;}void push(int k){ tree[tree[k].lc].bj+=tree[k].bj; tree[tree[k].rc].bj+=tree[k].bj; tree[tree[k].lc].sum+=tree[tree[k].lc].size*tree[k].bj; tree[tree[k].rc].sum+=tree[tree[k].rc].size*tree[k].bj; tree[k].bj=0;return ;}void change(int k,int l,int r,int z){ if(l<=tree[k].l&&tree[k].r<=r) { tree[k].bj+=z;tree[k].sum+=z*tree[k].size; return ; } if(tree[k].bj) push(k); int mid=(tree[k].l+tree[k].r)>>1; if(l<=mid) change(tree[k].lc,l,r,z); if(r>mid) change(tree[k].rc,l,r,z); tree[k].sum=tree[tree[k].lc].sum+tree[tree[k].rc].sum; return ;}int query(int k,int l,int r){ if(l<=tree[k].l&&tree[k].r<=r) return tree[k].sum; if(tree[k].bj) push(k); int mid=(tree[k].l+tree[k].r)>>1,total=0; if(l<=mid) total+=query(tree[k].lc,l,r); if(r>mid) total+=query(tree[k].rc,l,r); return total;}void slove(int k,int l,int r){ if(l==r){
s[l]=tree[k].sum;return ;} if(tree[k].bj) push(k); int mid=(l+r)>>1; slove(tree[k].lc,l,mid),slove(tree[k].rc,mid+1,r); return;}void dfs1(int u){ size[u]=1; for(int i=head[u];i;i=e[i].next) { int v=e[i].v; if(!fa[v]) fa[v]=u,deep[v]=deep[u]+1,dfs1(v),size[v]+=size[u]; } return ;}void dfs2(int u,int top1){ int k=0;pos[u]=++maxsize;top[u]=top1; for(int i=head[u];i;i=e[i].next) { int v=e[i].v; if(fa[v]==u&&size[v]>size[k]) k=v; } if(!k) return ; dfs2(k,top1); for(int i=head[u];i;i=e[i].next) { int v=e[i].v; if(fa[v]==u&&v!=k) dfs2(v,v); } return ;}void slovechange(int x,int y){ change(1,pos[y],pos[y],-1); while(top[x]!=top[y]) { if(deep[top[x]]
deep[y]) swap(x,y); change(1,pos[x],pos[y],1);return ;}int main(){ int x,y; n=read(); for(int i=1;i<=n;i++) a[i]=read(); for(int i=1;i<=n-1;i++) x=read(),y=read(),add(x,y),add(y,x); fa[1]=1;dfs1(1),dfs2(1,1),build(1,n); for(int i=1;i<=n-1;i++) slovechange(a[i],a[i+1]); slove(1,1,n); //for(int i=1;i<=n;i++) printf("%d\n",query(1,pos[i],pos[i])); for(int i=1;i<=n;i++) printf("%d\n",s[pos[i]]); return 0;}

转载于:https://www.cnblogs.com/nancheng58/p/10068079.html

你可能感兴趣的文章
ionic 入门学习
查看>>
[python]pickle和cPickle
查看>>
末日了,天是灰色的。
查看>>
Vuejs vm对象详解
查看>>
自定义RatingBar的一个问题(只显示显示一个星星)
查看>>
剑指Offer--二叉树的镜像
查看>>
PAT-BASIC-1031-查验身份证
查看>>
Python笔记5----集合set
查看>>
连连看小游戏
查看>>
js二级联动
查看>>
谜题32:循环者的诅咒
查看>>
RMI
查看>>
动态切换多数据源的配置
查看>>
win7电脑调整分区后分区不见的文件寻回法子
查看>>
《第一行代码》学习笔记2-Android开发特色
查看>>
bzoj3396 [Usaco2009 Jan]Total flow 水流
查看>>
20165231 2017-2018-2 《Java程序设计》第3周学习总结
查看>>
(180905)如何通过梯度下降法降低损失----Google机器学习速成课程笔记
查看>>
(响应式PC端媒体查询)电脑屏幕分辨率尺寸大全
查看>>
Crystal Reports拉报表报错:Error detected by database DLL
查看>>