【BZOJ-2325】道馆之战 树链剖分 + 线段树
2325: [ZJOI2011]道館之戰(zhàn)
Time Limit:?40 Sec??Memory Limit:?256 MBSubmit:?1153??Solved:?421
[Submit][Status][Discuss]
Description
口袋妖怪(又名神奇寶貝或?qū)櫸镄【`)紅/藍(lán)/綠寶石中的水系道館需要經(jīng)過(guò)三個(gè)冰地才能到達(dá)館主的面前,冰地中的每一個(gè)冰塊都只能經(jīng)過(guò)一次。當(dāng)一個(gè)冰地上的所有冰塊都被經(jīng)過(guò)之后,到下一個(gè)冰地的樓梯才會(huì)被打開(kāi)。三個(gè)冰地分別如下: 當(dāng)走出第三個(gè)冰地之后,就可以與館主進(jìn)行道館戰(zhàn)了。館主發(fā)現(xiàn)這個(gè)難度太小,導(dǎo)致經(jīng)常有挑戰(zhàn)者能通過(guò),為了加大難度,將道館分成了n個(gè)房間,每個(gè)房間中是兩個(gè)冰塊或障礙,表示一列冰地。任意兩個(gè)房間之間均有且僅有一條路徑相連,即這n個(gè)房間構(gòu)成一個(gè)樹(shù)狀結(jié)構(gòu)。每個(gè)房間分成了A和B兩個(gè)區(qū)域,每一區(qū)域都是一個(gè)薄冰塊或者障礙物。每次只能移動(dòng)到相鄰房間的同一類區(qū)域(即若你現(xiàn)在在這個(gè)房間的A區(qū)域,那么你只能移動(dòng)到相鄰房間的A區(qū)域)或這個(gè)房間的另一區(qū)域。現(xiàn)在挑戰(zhàn)者從房間u出發(fā),館主在房間v,那么挑戰(zhàn)者只能朝接近館主所在房間的方向過(guò)去。一開(kāi)始挑戰(zhàn)者可以在房間u的任意一個(gè)冰塊區(qū)域內(nèi)。如果挑戰(zhàn)者踩過(guò)的冰塊數(shù)達(dá)到了最大值(即沒(méi)有一種方案踩過(guò)的冰塊數(shù)更多了),那么當(dāng)挑戰(zhàn)者走到最后一個(gè)冰塊上時(shí),他會(huì)被瞬間傳送到館主面前與館主進(jìn)行道館戰(zhàn)。自從館主修改規(guī)則后已經(jīng)經(jīng)過(guò)了m天,每天要么是有一個(gè)挑戰(zhàn)者來(lái)進(jìn)行挑戰(zhàn),要么就是館主將某個(gè)房間進(jìn)行了修改。對(duì)于每個(gè)來(lái)的挑戰(zhàn)者,你需要計(jì)算出他若要和館主進(jìn)行戰(zhàn)斗需要經(jīng)過(guò)的冰塊數(shù)。Input
第一行包含兩個(gè)正整數(shù)n和m。第2行到第n行,每行包含兩個(gè)正整數(shù)x和y,表示一條連接房間x和房間y的邊。房間編號(hào)為1…n。接下來(lái)n行,每行包含兩個(gè)字符。第n + k行表示房間k的兩個(gè)區(qū)域,第一個(gè)字符為A區(qū)域,第二個(gè)字符為B區(qū)域。其中“.”(ASCII碼為46)表示是薄冰塊,“#”(ASCII碼為35)表示是障礙物。最后的m行,每行一個(gè)操作: l C u s:將房間u里的兩個(gè)區(qū)域修改為s。 l Q u v:詢問(wèn)挑戰(zhàn)者在房間u,館主在房間v時(shí),挑戰(zhàn)者能與館主進(jìn)行挑戰(zhàn)需要踩的冰塊數(shù)。如果房間u的兩個(gè)區(qū)域都是障礙物,那么輸出0。 N≤ 30 000 M ≤ 80 000Output
包含若干行,每行一個(gè)整數(shù)。即對(duì)于輸入中的每個(gè)詢問(wèn),依次輸出一個(gè)答案。
Sample Input
5 31 2
2 3
2 4
1 5
.#
..
#.
.#
..
Q 5 3
C 1 ##
Q 4 5
Sample Output
63
HINT
Source
Day2
Solution
樹(shù)鏈剖分+線段樹(shù)維護(hù)(類似)連通性。 ?和 堵塞的交通 維護(hù)方法類似。但是并沒(méi)有做過(guò)。
題目大意:每個(gè)節(jié)點(diǎn)分為ab連個(gè)塊,這些節(jié)點(diǎn)之間聯(lián)通呈樹(shù)形,每個(gè)節(jié)點(diǎn)的ab可能為障礙或者是空地,從一個(gè)節(jié)點(diǎn)上的a可以到這個(gè)節(jié)點(diǎn)的b,或者和這個(gè)節(jié)點(diǎn)聯(lián)通的所有a(前提是路徑上無(wú)障礙)。詢問(wèn)節(jié)點(diǎn)x到節(jié)點(diǎn)y的方向上最多走多少步,支持單點(diǎn)修改.
這樣如果是鏈上的情況,可以等價(jià)為一個(gè)2*N的矩陣,上面一行表示所有的a,下面一行表示所有的b,等效到樹(shù)上,就是多個(gè)2*1個(gè)矩陣連成的樹(shù)。
維護(hù)方法就是,對(duì)于一個(gè)區(qū)間,維護(hù)8個(gè)量分別表示:
dis[0][0]表示從區(qū)間左端點(diǎn)上到右端點(diǎn)上,dis[0][1]表示從區(qū)間左端點(diǎn)下到右端點(diǎn)上,dis[1][0]表示從區(qū)間左端點(diǎn)上到右端點(diǎn)下,dis[1][1]表示從區(qū)間左端點(diǎn)下到右端點(diǎn)下
Dis[0][0]表示從區(qū)間左端點(diǎn)上最多移動(dòng)數(shù),Dis[0][1]表示從區(qū)間左端點(diǎn)下最多移動(dòng)數(shù),Dis[1][0]表示從區(qū)間右端點(diǎn)上最多移動(dòng)數(shù),Dis[1][1]表示從區(qū)間右端點(diǎn)下最多移動(dòng)數(shù)
合并的時(shí)候就是討論一下,畫(huà)畫(huà)圖就可以得到。
有一個(gè)問(wèn)題,在詢問(wèn)路徑x-->y時(shí),有一邊到LCA的路徑是從下到上的,所以在和另一邊路徑合并答案時(shí)要先反向,然后就可以了。
寫(xiě)的時(shí)候要注意細(xì)節(jié),這個(gè)在Codeing的時(shí)候由于思路清晰,寫(xiě)的還是十分美觀的,就是斷斷續(xù)續(xù)自己 想+寫(xiě)+調(diào) 搞了近兩個(gè)小時(shí)...(話說(shuō)題意真是不清晰)
Code
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<cmath> using namespace std; inline int read() {int x=0; char ch=getchar(); while (ch<'0' || ch>'9') ch=getchar(); while (ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar(); return x; } #define MAXN 30010 #define INF 100000000 int N,M,id[MAXN]; char mp[MAXN][3]; namespace SegmentTree {struct SgtNode{int l,r,dis[2][2],Dis[2][2]; SgtNode() { l=0,r=0; memset(dis,0,sizeof(dis)); memset(Dis,0,sizeof(Dis)); }inline bool have() {return !dis[0][0] && !dis[0][1] && !dis[1][0] && !dis[1][1] && !Dis[0][0] && !Dis[0][1] && !Dis[1][0] && !Dis[1][1];}inline void rever() { swap(dis[0][1],dis[1][0]),swap(Dis[0][0],Dis[1][0]),swap(Dis[0][1],Dis[1][1]); } inline void modify(int a,int b){dis[0][0]=a? 1:-INF; Dis[0][0]=Dis[1][0]=a; dis[1][1]=b? 1:-INF; Dis[0][1]=Dis[1][1]=b;if (a && b) Dis[0][0]=Dis[0][1]=Dis[1][0]=Dis[1][1]=dis[0][1]=dis[1][0]=2;else dis[0][1]=dis[1][0]=-INF;}//單點(diǎn)修改時(shí)的賦值 }tree[MAXN<<2];// dis[0][0]表示從區(qū)間左端點(diǎn)上到右端點(diǎn)上,dis[0][1]表示從區(qū)間左端點(diǎn)下到右端點(diǎn)上,dis[1][0]表示從區(qū)間左端點(diǎn)上到右端點(diǎn)下,dis[1][1]表示從區(qū)間左端點(diǎn)下到右端點(diǎn)下// Dis[0][0]表示從區(qū)間左端點(diǎn)上最多移動(dòng)數(shù),Dis[0][1]表示從區(qū)間左端點(diǎn)下最多移動(dòng)數(shù),Dis[1][0]表示從區(qū)間右端點(diǎn)上最多移動(dòng)數(shù),Dis[1][1]表示從區(qū)間右端點(diǎn)下最多移動(dòng)數(shù) #define ls now<<1#define rs now<<1|1inline void update(SgtNode &rt,SgtNode lson,SgtNode rson){SgtNode tmp=rt;if (lson.have()) {rt=rson; rt.l=tmp.l; rt.r=tmp.r; return;} if (rson.have()) {rt=lson; rt.l=tmp.l; rt.r=tmp.r; return;}rt.dis[0][0]=max( max( lson.dis[0][0]+rson.dis[0][0] , lson.dis[0][1]+rson.dis[1][0] ) , -INF);rt.dis[0][1]=max( max( lson.dis[0][0]+rson.dis[0][1] , lson.dis[0][1]+rson.dis[1][1] ) , -INF);rt.dis[1][0]=max( max( lson.dis[1][0]+rson.dis[0][0] , lson.dis[1][1]+rson.dis[1][0] ) , -INF);rt.dis[1][1]=max( max( lson.dis[1][1]+rson.dis[1][1] , lson.dis[1][0]+rson.dis[0][1] ) , -INF);rt.Dis[0][0]=max( lson.Dis[0][0] , max( lson.dis[0][0]+rson.Dis[0][0] , lson.dis[0][1]+rson.Dis[0][1] ) );rt.Dis[0][1]=max( lson.Dis[0][1] , max( lson.dis[1][0]+rson.Dis[0][0] , lson.dis[1][1]+rson.Dis[0][1] ) );rt.Dis[1][0]=max( rson.Dis[1][0] , max( rson.dis[0][0]+lson.Dis[1][0] , rson.dis[1][0]+lson.Dis[1][1] ) );rt.Dis[1][1]=max( rson.Dis[1][1] , max( rson.dis[0][1]+lson.Dis[1][0] , rson.dis[1][1]+lson.Dis[1][1] ) );}inline void Update(int now) {update(tree[now],tree[ls],tree[rs]);}inline void Modify(int now,int pos,char MP[]){int l=tree[now].l,r=tree[now].r;if (l==r) {tree[now].modify(MP[1]=='.',MP[2]=='.'); return;}int mid=(l+r)>>1;if (pos<=mid) Modify(ls,pos,MP); else Modify(rs,pos,MP);Update(now);}inline SgtNode Query(int now,int L,int R){int l=tree[now].l,r=tree[now].r;if (l==L && R==r) return tree[now];int mid=(l+r)>>1; SgtNode ret;if (R<=mid) return Query(ls,L,R);else if (L>mid) return Query(rs,L,R);else return update(ret,Query(ls,L,mid),Query(rs,mid+1,R)),ret;}inline void Build(int now,int l,int r){tree[now].l=l; tree[now].r=r;if (l==r) return;int mid=(l+r)>>1;Build(ls,l,mid); Build(rs,mid+1,r);} } using namespace SegmentTree; namespace Divide {struct EdgeNode{int next,to;}edge[MAXN<<1];int head[MAXN],cnt=1;inline void AddEdge(int u,int v) {cnt++; edge[cnt].to=v; edge[cnt].next=head[u]; head[u]=cnt;}inline void InsertEdge(int u,int v) {AddEdge(u,v); AddEdge(v,u);}int deep[MAXN],fa[MAXN],size[MAXN],son[MAXN],pl[MAXN],dfn,top[MAXN];inline void DFS_1(int now,int last){size[now]=1;for (int i=head[now]; i; i=edge[i].next)if (edge[i].to!=last){fa[edge[i].to]=now;deep[edge[i].to]=deep[now]+1;DFS_1(edge[i].to,now);size[now]+=size[edge[i].to];if (size[son[now]]<size[edge[i].to]) son[now]=edge[i].to;}}inline void DFS_2(int now,int chain){pl[now]=++dfn; top[now]=chain;if (son[now]) DFS_2(son[now],chain);for (int i=head[now]; i; i=edge[i].next)if (edge[i].to!=fa[now] && edge[i].to!=son[now])DFS_2(edge[i].to,edge[i].to);}inline void Modify(int pos,char MP[]) {SegmentTree::Modify(1,pl[pos],MP);}inline void GetAns(int x,int y){SgtNode la,ra,ans; ra.l=ra.r=1; la.l=la.r=1;while (top[x]!=top[y])if (deep[top[x]]<deep[top[y]]) update(ra,Query(1,pl[top[y]],pl[y]),ra),y=fa[top[y]];else update(la,Query(1,pl[top[x]],pl[x]),la),x=fa[top[x]];if (deep[x]<deep[y]) update(ra,Query(1,pl[x],pl[y]),ra);elseupdate(la,Query(1,pl[y],pl[x]),la);la.rever(); update(ans,la,ra);printf("%d\n",max(ans.Dis[0][0],ans.Dis[0][1]));}inline void BuildTree() {for (int i=1; i<=N; i++) Modify(i,mp[i]);} } using namespace Divide; int main() { // freopen("fight.in","r",stdin); freopen("fight.out","w",stdout);N=read(),M=read();for (int x,y,i=1; i<=N-1; i++) x=read(),y=read(),InsertEdge(x,y);Divide::DFS_1(1,0); Divide::DFS_2(1,1);for (int i=1; i<=N; i++) scanf("%s",mp[i]+1);SegmentTree::Build(1,1,N); Divide::BuildTree();while (M--){char opt[2],MP[3]; scanf("%s",opt); int L,R,p;switch (opt[0]){case 'Q' : L=read(),R=read(); Divide::GetAns(L,R); break;case 'C' : p=read(); scanf("%s",MP+1); Divide::Modify(p,MP); break;}}return 0; }轉(zhuǎn)載于:https://www.cnblogs.com/DaD3zZ-Beyonder/p/5978653.html
總結(jié)
以上是生活随笔為你收集整理的【BZOJ-2325】道馆之战 树链剖分 + 线段树的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 2012年6月的标致4008,电瓶什么型
- 下一篇: Codeforces Round #37