zoukankan      html  css  js  c++  java
  • P3258 [JLOI2014]松鼠的新家(树上差分板子)

    树上差分!!!

    这是一道模板题的稍微改编

    首先,啥是树上差分

    就和名字一样,书上的差分

    倘若我们要给(u,v)的路径上加上k

    则给u+k,v+k,然后lca(u,v)-k(因为事实上lca处就加了一个,可是这么干加了俩) father(lca(---))-k(这里及更高处根本没动过,消除影响)

    所以这题就是lca+树上差分

    但是和模板题的不同是,这里,每一点(除了一开始的和重点),即使起始点也是终点,而且最后一点不用放。

    需要特别处理

          #include<iostream>
    #include<cstdio>
    using namespace std;
    int head[2*500001];
    int p;
    struct b{
    	int to;
    	int ne;
    } e[2*500001];
    int fa[2*500001][30];
    int dep[2*500001];
    int n,m,s;
    int logg[2*500001];
    int ma[300005]; 
    int cha[300005];
    int ans[300005];
    int cnt[300005];
    int x,y;
    void dfs(int now,int fat){
    	dep[now]=dep[fat]+1;
    	fa[now][0]=fat;
    	for(int i=1;(1<<i)<=dep[now];++i){
    		fa[now][i]=fa[fa[now][i-1]][i-1];
    	}
    	for(int i=head[now];i;i=e[i].ne){
    		if(e[i].to!=fat)
    		dfs(e[i].to,now);
    	}
    }
    int lca(int x,int y){
    	if(dep[x]<dep[y]){
    		swap(x,y);
    	}
    	while(dep[x]>dep[y])
    	x=fa[x][logg[dep[x]-dep[y]]-1];
    	if(x==y)
    	return x;
    	for(int k=logg[dep[x]]-1;k>=0;k--){
    		if(fa[x][k]!=fa[y][k]){
    			x=fa[x][k];
    			y=fa[y][k];
    		}
    	}
    	return fa[x][0];
    }
    void print(int x,int fa){
    //	cout<<x<<endl;
    	for(int i=head[x];i;i=e[i].ne){
    //			cout<<e[i].to;
    		if(e[i].to!=fa){
    			print(e[i].to,x);
    			ans[x]+=ans[e[i].to];
    		}
    	}
    	ans[x]+=cha[x];
    	return ;
    }
    void add(int f,int t){
    	p++;
    	e[p].to=t;
    	e[p].ne=head[f];
    	head[f]=p;
    }
    int main(){
    	scanf("%d",&n);
    	for(int i=1;i<=n;++i){
    		scanf("%d",&ma[i]);
    		if(i!=1)
    		cnt[ma[i]]++;
    		}
    	for(int i=1;i<=n-1;++i){
    		scanf("%d%d",&x,&y);
    		add(x,y);
    		add(y,x);
    	}
    	dfs(1,0);
    	for(int i=1;i<=n;++i)
    		logg[i]=logg[i-1]+(1<<logg[i-1]==i);
    	for(int i=1;i<n;++i){
    		int o=lca(ma[i],ma[i+1]);
    		int v=ma[i+1];
    		cha[ma[i]]++;
    		cha[v]++;
    		cha[o]--;
    		cha[fa[o][0]]--;
    	}
    	print(1,0);
    	for(int i=1;i<=n;++i){
    		ans[i]-=cnt[i];
    		cout<<ans[i]<<endl;
    	}
    	return 0;
    }
    
  • 相关阅读:
    用elasticsearch分析中国大学省份分布
    【翻译】Kinect v1和Kinect v2的彻底比较
    翻译 Tri-Ace:在Shader里近似渲染公式
    翻译 基于物理渲染的美术资源设计流程
    翻译 次世代基于物理渲染的反射模型
    关于Depth Bounds Test (DBT)和在CE3的运用
    使用Xcode GPU Frame Caputre教程
    如何使用Xcode分析调试在真机运行的UE4 IOS版游戏
    个人翻译的cedec2010基于物理的光照
    使用Nsight查找CE3的渲染bug
  • 原文地址:https://www.cnblogs.com/For-Miku/p/13417633.html
Copyright © 2011-2022 走看看