BSOJ1495(BZOJ2157) 旅遊 樹鏈剖分

1495 -- 【BZOJ2157】旅遊
Description
  Ray 樂忠於旅遊,這次他來到了T 城。T 城是一個水上城市,一共有 N 個景點,有些景點之間會用一座橋連接。爲了方便遊客到達每個景點但又爲了節約成本,T 城的任意兩個景點之間有且只有一條路徑。換句話說, T 城中只有N ? 1 座橋。Ray 發現,有些橋上可以看到美麗的景色,讓人心情愉悅,但有些橋狹窄泥濘,令人煩躁。於是,他給每座橋定義一個愉悅度w,也就是說,Ray 經過這座橋會增加w 的愉悅度,這或許是正的也可能是負的。有時,Ray 看待同一座橋的心情也會發生改變。現在,Ray 想讓你幫他計算從u 景點到v 景點能獲得的總愉悅度。有時,他還想知道某段路上最美麗的橋所提供的最大愉悅度,或是某段路上最糟糕的一座橋提供的最低愉悅度。
Input
輸入的第一行包含一個整數N,表示T 城中的景點個數。景點編號爲 0...N ? 1。
接下來N ? 1 行,每行三個整數u、v 和w,表示有一條u 到v,使 Ray 愉悅度增加w 的橋。橋的編號爲1...N ? 1。|w| <= 1000。
輸入的第N + 1 行包含一個整數M,表示Ray 的操作數目。
接下來有M 行,每行描述了一個操作,操作有如下五種形式: C i w,表示Ray 對於經過第i 座橋的愉悅度變成了w。 N u v,表示Ray 對於經過景點u 到v 的路徑上的每一座橋的愉悅度都變成原來的相反數。 SUM u v,表示詢問從景點u 到v 所獲得的總愉悅度。 MAX u v,表示詢問從景點u 到v 的路徑上的所有橋中某一座橋所提供的最大愉悅度。 MIN u v,表示詢問從景點u 到v 的路徑上的所有橋中某一座橋所提供的最小愉悅度。測試數據保證,任意時刻,Ray 對於經過每一座橋的愉悅度的絕對值小於等於1000。
Output
對於每一個詢問(操作S、MAX 和MIN),輸出答案。
Sample Input
30
1 1
1 2 2
8
SUM 0 2
MAX 0 2
N 0 1
SUM 0 2
MIN 0 2
C 1 3
SUM 0 2
MAX 0 2
Sample Output
3
2
1
-1
5
3
Hint

一共有10 個數據,對於第i (1 <= i <= 10) 個數據, N = M = i * 2000。


樹鏈剖分模板題,都是很基礎的操作..

好久沒寫過樹剖了,調了半天,發現一開始size沒設爲1,還少寫了一句DFS2,有負邊權maxx初值竟然設爲0...

一臉尷尬,本來準備快速AC的,慘遭打臉。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
//#include<cmath>
#define L(x) (x<<1)
#define R(x) (x<<1|1)
#define INF 0x3fffffff
using namespace std;
struct seg
{
	int l,r,maxx,minn,sum,lazy;
}tree[100005];
struct TREE
{
	int x,y,val;
}T[50005];
struct node
{
	int to,next,val;
}w[1000005];
int h[1000005]={0},cnt=0,dep[200005]={0},fa[1000005]={0},size[1000005]={0},top[1000005]={0},tid[1000005]={0},tim=0,n,son[100005]={0};
void Addarc(int x,int y,int val)
{
	cnt++;w[cnt].to=y;w[cnt].next=h[x];h[x]=cnt;
}
void DFS1(int x,int f,int d)
{
	dep[x]=d;
	fa[x]=f;
	size[x]=1;
	for(int i=h[x];i;i=w[i].next)
	{
		int to=w[i].to;
		if(to==f)continue;
		DFS1(to,x,d+1);
		size[x]+=size[to];
		if(son[x]==0||(size[son[x]]<size[to]))son[x]=to;
	}
}
void DFS2(int x,int tp)
{
	top[x]=tp;
	tid[x]=++tim;
	if(son[x]==0)return;
	DFS2(son[x],tp);
	for(int i=h[x];i;i=w[i].next)
	{
		int to=w[i].to;
		if(to==fa[x]||to==son[x])continue;
		DFS2(to,to);
	}
}
void Pushup(int root)
{
	tree[root].sum=tree[L(root)].sum+tree[R(root)].sum;
	tree[root].maxx=max(tree[L(root)].maxx,tree[R(root)].maxx);
	tree[root].minn=min(tree[L(root)].minn,tree[R(root)].minn);
}
void Pushdown(int root)
{
	if(tree[root].lazy)
	{
		tree[root].lazy=0;
		tree[L(root)].lazy^=1;
		tree[R(root)].lazy^=1;
		tree[L(root)].sum*=-1;
		tree[R(root)].sum*=-1;
		tree[L(root)].maxx*=-1;tree[L(root)].minn*=-1;swap(tree[L(root)].maxx,tree[L(root)].minn);
		tree[R(root)].maxx*=-1;tree[R(root)].minn*=-1;swap(tree[R(root)].maxx,tree[R(root)].minn);
	}
}
void Build(int l,int r,int root)
{
	tree[root].l=l;tree[root].r=r;
        tree[root].maxx=-INF;tree[root].minn=INF;tree[root].lazy=0;tree[root].sum=0;
        if(l==r)return;
        int mid=(l+r)>>1;
        Pushdown(root);
	Build(l,mid,L(root));
	Build(mid+1,r,R(root));
	Pushup(root);
}
void Change(int x,int data,int root)
{
	if(tree[root].l>x||tree[root].r<x)return;
	if(tree[root].l==tree[root].r&&tree[root].l==x)
	{
		tree[root].maxx=data;
		tree[root].minn=data;
		tree[root].sum=data;
		return;
	}
	Pushdown(root);
	Change(x,data,L(root));
	Change(x,data,R(root));
	Pushup(root);
}
void Neg(int l,int r,int root)
{
	if(tree[root].l>r||tree[root].r<l)return;
	if(l<=tree[root].l&&tree[root].r<=r)
	{
		tree[root].lazy^=1;
		tree[root].maxx*=-1;
		tree[root].minn*=-1;
		tree[root].sum*=-1;
		swap(tree[root].minn,tree[root].maxx);
		return;
	}
	Pushdown(root);
	Neg(l,r,L(root));
	Neg(l,r,R(root));
	Pushup(root);
}
int Ask_sum(int l,int r,int root)
{
	if(tree[root].l>r||tree[root].r<l)return 0;
	if(l<=tree[root].l&&tree[root].r<=r)return tree[root].sum;
	int sum=0;
	Pushdown(root);
	sum+=Ask_sum(l,r,L(root));
	sum+=Ask_sum(l,r,R(root));
	Pushup(root);
	return sum;
}
int Ask_max(int l,int r,int root)
{
	if(tree[root].l>r||tree[root].r<l)return -INF;
	if(l<=tree[root].l&&tree[root].r<=r)return tree[root].maxx;
	int maxx=-INF;
	Pushdown(root);
	maxx=max(maxx,Ask_max(l,r,L(root)));
	maxx=max(maxx,Ask_max(l,r,R(root)));
	Pushup(root);
	return maxx;
}
int Ask_min(int l,int r,int root)
{
	if(tree[root].l>r||tree[root].r<l)return INF;
	if(l<=tree[root].l&&tree[root].r<=r)return tree[root].minn;
	int minn=INF;
	Pushdown(root);
	minn=min(minn,Ask_min(l,r,L(root)));
	minn=min(minn,Ask_min(l,r,R(root)));
	Pushup(root);
	return minn;
}
void Doneg(int x,int y)
{
	int t1=top[x],t2=top[y];
	while(t1!=t2)
	{
		if(dep[t1]<dep[t2]){swap(x,y);swap(t1,t2);}
		Neg(tid[t1],tid[x],1);
		x=fa[t1];t1=top[x];
	}
	if(x==y)return;
	if(dep[x]>dep[y])swap(x,y);
	Neg(tid[son[x]],tid[y],1);
}
int Getsum(int x,int y)
{
	int t1=top[x],t2=top[y],ans=0;
	while(t1!=t2)
	{
		if(dep[t1]<dep[t2]){swap(x,y);swap(t1,t2);}
		ans+=Ask_sum(tid[t1],tid[x],1);
		x=fa[t1];t1=top[x];
	}
	if(x==y)return ans;
	if(dep[x]>dep[y])swap(x,y);
	ans+=Ask_sum(tid[son[x]],tid[y],1);
	return ans;
}
int Getmax(int x,int y)
{
	int t1=top[x],t2=top[y],maxx=-INF;
	while(t1!=t2)
	{
		if(dep[t1]<dep[t2]){swap(x,y);swap(t1,t2);}
		maxx=max(maxx,Ask_max(tid[t1],tid[x],1));
		x=fa[t1];t1=top[x];
	}
	if(x==y)return maxx;
	if(dep[x]>dep[y])swap(x,y);
	maxx=max(maxx,Ask_max(tid[son[x]],tid[y],1));
	return maxx;
}
int Getmin(int x,int y)
{
	int t1=top[x],t2=top[y],minn=INF;
	while(t1!=t2)
	{
		if(dep[t1]<dep[t2]){swap(x,y);swap(t1,t2);}
		minn=min(minn,Ask_min(tid[t1],tid[x],1));
		x=fa[t1];t1=top[x];
	}
	if(x==y)return minn;
	if(dep[x]>dep[y])swap(x,y);
	minn=min(minn,Ask_min(tid[son[x]],tid[y],1));
	return minn;
}
void init()
{
	scanf("%d",&n);
	for(int i=1;i<n;i++)
	{
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		x++;y++;
		T[i].x=x;T[i].y=y;T[i].val=z;
		Addarc(T[i].x,T[i].y,T[i].val);
		Addarc(T[i].y,T[i].x,T[i].val);
	}
	DFS1(1,0,0);
	DFS2(1,1);
	Build(1,n,1);
	for(int i=1;i<n;i++)
	{
		if(dep[T[i].x]>dep[T[i].y])swap(T[i].x,T[i].y);
		Change(tid[T[i].y],T[i].val,1);
	}
}
void work()
{
	int m,a,b;
	string s;
	scanf("%d",&m);
	for(int i=1;i<=m;i++)
	{
		cin>>s;
		scanf("%d%d",&a,&b);
		if(s=="C")Change(tid[T[a].y],b,1);
		if(s=="N")Doneg(++a,++b);
                if(s=="SUM")printf("%d\n",Getsum(++a,++b));
     	        if(s=="MAX")printf("%d\n",Getmax(++a,++b));
		if(s=="MIN")printf("%d\n",Getmin(++a,++b));
	}
}
int main()
{
	init();
	work();
return 0;
}
/*
3
0 1 1
1 2 2
8
SUM 0 2
MAX 0 2
N 0 1
SUM 0 2
MIN 0 2
C 1 3
SUM 0 2
MAX 0 2
*/
稍微注意下這個細節,因爲是邊權所以記錄在兒子頭上,所以最後一步修改是xxx(tid[son[x]],tid[y],1)

還有就是題目是邊權的情況下不要忘了x==y返回這個條件。

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章