zoukankan      html  css  js  c++  java
  • poj3237 Tree

    Description

    You are given a tree with N nodes. The tree’s nodes are numbered 1 through N and its edges are numbered 1 through N − 1. Each edge is associated with a weight. Then you are to execute a series of instructions on the tree. The instructions can be one of the following forms:

    CHANGE i v Change the weight of the ith edge to v
    NEGATE a b Negate the weight of every edge on the path from a to b
    QUERY a b Find the maximum weight of edges on the path from a to b

    Input

    The input contains multiple test cases. The first line of input contains an integer t (t ≤ 20), the number of test cases. Then follow the test cases.

    Each test case is preceded by an empty line. The first nonempty line of its contains N (N ≤ 10,000). The next N − 1 lines each contains three integers ab and c, describing an edge connecting nodes a and b with weight c. The edges are numbered in the order they appear in the input. Below them are the instructions, each sticking to the specification above. A lines with the word “DONE” ends the test case.

    Output

    For each “QUERY” instruction, output the result on a separate line.

    Sample Input

    1
    
    3
    1 2 1
    2 3 2
    QUERY 1 2
    CHANGE 1 3
    QUERY 1 2
    DONE

    Sample Output

    1
    3

    好惨啊……有一个加1的没写结果最裸的树链剖分调了一整天
    线段树维护最大值和最小值,对于一条链上权值取反的操作只要mx=-mn,mn=-mx即可
    #include<cstdio>
    #include<iostream>
    #include<cstring>
    #define LL long long
    #define inf 1000000000
    #define N 10010
    using namespace std;
    int bin[15]={1,2,4,8,16,32,64,128,256,512,1024,2048,4096,8192,16384};
    inline LL read()
    {
        LL 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-'0';ch=getchar();}
        return x*f;
    }
    struct segtree{int l,r,mx,mn,tag;}tree[4*N];
    struct edge{int to,next,v;}e[2*N];
    int T,n,cnt,tt,mx;
    char ch[10];
    int head[N];
    int fa[N][15],son[N],depth[N],v[N];
    bool mrk[N];
    int place[N],pplace[N],belong[N];
    int query[N];
    inline void revswap(int &a,int &b){int t=a;a=-b;b=-t;}
    inline void ins(int u,int v,int w)
    {
    	e[++cnt].to=v;
    	e[cnt].v=w;
    	e[cnt].next=head[u];
    	head[u]=cnt;
    }
    inline void insert(int u,int v,int w)
    {
    	ins(u,v,w);
    	ins(v,u,w);
    }
    inline void dfs(int x,int dep)
    {
    	if (mrk[x])return;mrk[x]=1;
    	depth[x]=dep;son[x]=1;
    	for (int i=1;i<15;i++)
    		if (bin[i]<=depth[x])
    		  	fa[x][i]=fa[fa[x][i-1]][i-1];
    		else break;
    	for (int i=head[x];i;i=e[i].next)
    		if (!mrk[e[i].to])
    		{
    			fa[e[i].to][0]=x;
    			dfs(e[i].to,dep+1);
    			v[e[i].to]=e[i].v;
    			son[x]+=son[e[i].to];
    			query[(i+1)>>1]=e[i].to;
    		} 
    }
    inline void dfs2(int x,int chain)
    {
    	place[x]=++tt;pplace[tt]=x;
    	belong[x]=chain;
    	int res=0,mx=-inf;
    	for (int i=head[x];i;i=e[i].next)
    		if (fa[x][0]!=e[i].to)
    		{
    			if (son[e[i].to]>mx)
    			{
    				mx=son[e[i].to];
    				res=e[i].to;
    			}
    		}
    	if (!res)return;
    	dfs2(res,chain);
    	for (int i=head[x];i;i=e[i].next)
    		if (fa[x][0]!=e[i].to&&res!=e[i].to)
    			dfs2(e[i].to,e[i].to);
    }
    inline int LCA(int a,int b)
    {
    	if (depth[a]<depth[b])swap(a,b);
    	int res=depth[a]-depth[b];
    	for (int i=0;i<15;i++)
    		if (res & bin[i]) a=fa[a][i];
    	for (int i=14;i>=0;i--)
    		if (fa[a][i]!=fa[b][i])
    		{
    			a=fa[a][i];
    			b=fa[b][i];
    		}
    	if (a==b)return a;
    	return fa[a][0];
    }
    inline void pushdown(int k)
    {
    	tree[k].tag=0;
    	tree[k<<1].tag^=1;tree[k<<1|1].tag^=1;
    	revswap(tree[k<<1].mx,tree[k<<1].mn);
    	revswap(tree[k<<1|1].mx,tree[k<<1|1].mn);
    }
    inline void update(int k)
    {
    	tree[k].mx=max(tree[k<<1].mx,tree[k<<1|1].mx);
    	tree[k].mn=min(tree[k<<1].mn,tree[k<<1|1].mn);
    }
    inline void buildtree(int now,int l,int r)
    {
    	tree[now].l=l;tree[now].r=r;
    	if (l==r)
    	{
    		tree[now].mx=tree[now].mn=v[pplace[l]];
    		return;
    	}
    	int mid=(l+r)>>1;
    	buildtree(now<<1,l,mid);
    	buildtree(now<<1|1,mid+1,r);
    	update(now);
    }
    inline void change_in_tree(int now,int x,int y)
    {	
    	int l=tree[now].l,r=tree[now].r;
    	if (l!=r&&tree[now].tag)pushdown(now);
    	if (l==r)
    	{
    		tree[now].mx=tree[now].mn=y;
    		return;
    	}
    	int mid=(l+r)>>1;
    	if (x<=mid)change_in_tree(now<<1,x,y);
    	else change_in_tree(now<<1|1,x,y);
    	update(now);
    }
    inline void negate_in_tree(int now,int x,int y)
    {
    	int l=tree[now].l,r=tree[now].r;
    	if (l!=r&&tree[now].tag)pushdown(now);
    	if (l==x&&r==y)
    	{
    		revswap(tree[now].mx,tree[now].mn);
    		tree[now].tag=1;
    		return;
    	}
    	int mid=(l+r)>>1;
    	if (y<=mid)negate_in_tree(now<<1,x,y);
    	else if (x>mid)negate_in_tree(now<<1|1,x,y);
    	else
    	{
    		negate_in_tree(now<<1,x,mid);
    		negate_in_tree(now<<1|1,mid+1,y);
    	}
    	update(now);
    }
    inline void ask_in_tree(int now,int x,int y)
    {
    	int l=tree[now].l,r=tree[now].r;
    	if (l!=r&&tree[now].tag)pushdown(now);
    	if (l==x&&r==y)
    	{
    		mx=max(mx,tree[now].mx);
    		return;
    	}
    	int mid=(l+r)>>1;
    	if (y<=mid)ask_in_tree(now<<1,x,y);
    	else if (x>mid)ask_in_tree(now<<1|1,x,y);
    	else 
    	{
    		ask_in_tree(now<<1,x,mid);
    		ask_in_tree(now<<1|1,mid+1,y);
    	}
    }
    inline void reverse(int from,int to)
    {
    	if (from==to)return;
    	int l,r;
    	while (belong[from]!=belong[to])
    	{
    		l=place[belong[from]];r=place[from];
    		negate_in_tree(1,l,r);
    		from=fa[belong[from]][0];
    	}
    	if (place[to]+1<=place[from])
    	negate_in_tree(1,place[to]+1,place[from]);
    }
    inline int ask(int from,int to)
    {
    	int l,r;mx=-inf;
    	if (from==to)return mx;
    	while (belong[from]!=belong[to])
    	{
    		l=place[belong[from]];r=place[from];
    		ask_in_tree(1,l,r);
    		from=fa[belong[from]][0];
    	}
    	if (place[to]+1<=place[from])
    	ask_in_tree(1,place[to]+1,place[from]);
    	return mx;
    }
    inline void work()
    {
    	memset(tree,0,sizeof(tree));
    	memset(head,0,sizeof(head));
    	memset(e,0,sizeof(e));
    	memset(mrk,0,sizeof(mrk));
    	memset(query,0,sizeof(query));
    	memset(place,0,sizeof(place));
    	memset(pplace,0,sizeof(pplace));
    	memset(belong,0,sizeof(belong));
    	memset(fa,0,sizeof(fa));
    	memset(son,0,sizeof(son));
    	memset(depth,0,sizeof(depth));
    	memset(v,0,sizeof(v));
    	tt=cnt=0;
    	n=read();
    	for (int i=1;i<n;i++)
    	{
    		int x=read(),y=read(),z=read();
    		insert(x,y,z);
    	}
    	dfs(1,0);
    	dfs2(1,1);
    	buildtree(1,1,n);
    	while (scanf("%s",ch+1)&&ch[1]!='D')
    	{
    		int a=read(),b=read();
    		if (ch[1]=='C')change_in_tree(1,place[query[a]],b);
    		if (ch[1]=='N')
    		{
    			int c=LCA(a,b);
    			reverse(a,c);
    			reverse(b,c);
    		}
    		if (ch[1]=='Q')
    		{
    			int c=LCA(a,b);
    			printf("%d
    ",max( ask(a,c),ask(b,c) ));
    		}
    	}
    }
    int main()
    {
    	T=read();
    	while (T--)work();
    }
    
    ——by zhber,转载请注明来源
  • 相关阅读:
    如何解决秒杀的性能问题和超卖的讨论
    Redis作者谈Redis应用场景
    Net开源工作流Roadflow的使用与集成
    VisualStudioOnline协同工作流程
    中文乱码?不,是 HTML 实体编码!(转)
    Nginx搭建反向代理服务器过程详解(转)
    2016年1月15日面试某互联网公司总结(转)
    浏览器缓存知识小结及应用(转)
    Gradle多项目配置的一个demo
    Android getResources的作用和须要注意点
  • 原文地址:https://www.cnblogs.com/zhber/p/4177023.html
Copyright © 2011-2022 走看看