【BZOJ1969】[Ahoi2005]LANE 航線規劃 離線+樹鏈剖分+線段樹
阿新 • • 發佈:2017-07-16
個數 wap 鎖定 樹邊 mes zoj swap 相同 swa
第一行有兩個整數N,M。表示有N個星球(1< N < 30000),初始時已經有M條航線(1 < M < 100000)。隨後有M行,每行有兩個不相同的整數A、B表示在星球A與B之間存在一條航線。接下來每行有三個整數C、A、B。C為1表示詢問當前星球A和星球B之間有多少條關鍵航線;C為0表示在星球A和星球B之間的航線被破壞,當後面再遇到C為1的情況時,表示詢問航線被破壞後,關鍵路徑的情況,且航線破壞後不可恢復; C為-1表示輸入文件結束,這時該行沒有A,B的值。被破壞的航線數目與詢問的次數總和不超過40000。
5 5
1 2
1 3
3 4
4 5
4 2
1 1 5
0 4 2
1 5 1
-1
3
【BZOJ1969】[Ahoi2005]LANE 航線規劃
Description
對Samuel星球的探險已經取得了非常巨大的成就,於是科學家們將目光投向了Samuel星球所在的星系——一個巨大的由千百萬星球構成的Samuel星系。 星際空間站的Samuel II巨型計算機經過長期探測,已經鎖定了Samuel星系中許多星球的空間坐標,並對這些星球從1開始編號1、2、3……。 一些先遣飛船已經出發,在星球之間開辟探險航線。 探險航線是雙向的,例如從1號星球到3號星球開辟探險航線,那麽從3號星球到1號星球也可以使用這條航線。 例如下圖所示: 在5個星球之間,有5條探險航線。 A、B兩星球之間,如果某條航線不存在,就無法從A星球抵達B星球,我們則稱這條航線為關鍵航線。 顯然上圖中,1號與5號星球之間的關鍵航線有1條:即為4-5航線。 然而,在宇宙中一些未知的磁暴和行星的沖撞,使得已有的某些航線被破壞,隨著越來越多的航線被破壞,探險飛船又不能及時回復這些航線,可見兩個星球之間的關鍵航線會越來越多。 假設在上圖中,航線4-2(從4號星球到2號星球)被破壞。此時,1號與5號星球之間的關鍵航線就有3條:1-3,3-4,4-5。 小聯的任務是,不斷關註航線被破壞的情況,並隨時給出兩個星球之間的關鍵航線數目。現在請你幫助完成。Input
Output
對每個C為1的詢問,輸出一行一個整數表示關鍵航線數目。 註意:我們保證無論航線如何被破壞,任意時刻任意兩個星球都能夠相互到達。在整個數據中,任意兩個星球之間最多只可能存在一條直接的航線。Sample Input
1 2
1 3
3 4
4 5
4 2
1 1 5
0 4 2
1 5 1
-1
Sample Output
13
題解:既然能離線搞肯定要離線啦~由於保證了最後時刻圖是聯通的,所以我們先求出一棵生成樹,然後反著往樹上加邊。每加一條邊,就相當於將這條邊覆蓋的所有樹邊都打上標記,意味著它們不再是關鍵路徑。查詢時,看一下路徑上有多少條未標記的邊就行了。
#include <cstdio> #include <iostream> #include <cstring> #include <map> #include <utility> #define MP(A,B) make_pair(A,B) #define lson x<<1 #define rson x<<1|1 using namespace std; const int maxm=200010; const int maxn=30010; typedef pair<int,int> pii; int n,m,q,cnt; int to[maxm],next[maxm],head[maxn],val[maxm],dep[maxn],fa[maxn],top[maxn],son[maxn],siz[maxn],p[maxn]; int s[maxn<<2],tag[maxn<<2],qa[maxm],qb[maxm],op[maxm],ans[maxm],pa[maxm],pb[maxm]; map<pii,int> mp; int rd() { int ret=0,f=1; char gc=getchar(); while(gc<‘0‘||gc>‘9‘) {if(gc==‘-‘)f=-f; gc=getchar();} while(gc>=‘0‘&&gc<=‘9‘) ret=ret*10+gc-‘0‘,gc=getchar(); return ret*f; } void build(int l,int r,int x) { if(l==r) { s[x]=(l!=1); return ; } int mid=l+r>>1; build(l,mid,lson),build(mid+1,r,rson); s[x]=s[lson]+s[rson]; } void pushdown(int x) { if(tag[x]) tag[lson]=tag[rson]=1,s[lson]=s[rson]=0,tag[x]=0; } void updata(int l,int r,int x,int a,int b) { if(a<=l&&r<=b) { s[x]=0,tag[x]=1; return ; } pushdown(x); int mid=l+r>>1; if(a<=mid) updata(l,mid,lson,a,b); if(b>mid) updata(mid+1,r,rson,a,b); s[x]=s[lson]+s[rson]; } int query(int l,int r,int x,int a,int b) { if(a<=l&&r<=b) return s[x]; pushdown(x); int mid=l+r>>1; if(b<=mid) return query(l,mid,lson,a,b); if(a>mid) return query(mid+1,r,rson,a,b); return query(l,mid,lson,a,b)+query(mid+1,r,rson,a,b); } void add(int a,int b) { mp[MP(max(a,b),min(a,b))]=cnt,to[cnt]=b,next[cnt]=head[a],head[a]=cnt++; to[cnt]=a,next[cnt]=head[b],head[b]=cnt++; } void dfs1(int x) { siz[x]=1; for(int i=head[x];i!=-1;i=next[i]) { if(!dep[to[i]]&&!val[i]) { dep[to[i]]=dep[x]+1,fa[to[i]]=x,dfs1(to[i]),siz[x]+=siz[to[i]]; if(siz[to[i]]>siz[son[x]]) son[x]=to[i]; } else { if(to[i]!=fa[x]&&!val[i]&&!val[i^1]) pa[++pa[0]]=x,pb[pa[0]]=to[i]; val[i]=1; } } } void dfs2(int x,int tp) { top[x]=tp,p[x]=++p[0]; if(son[x]) dfs2(son[x],tp); for(int i=head[x];i!=-1;i=next[i]) if(!val[i]&&to[i]!=son[x]) dfs2(to[i],to[i]); } void modify(int x,int y) { while(top[x]!=top[y]) { if(dep[top[x]]<dep[top[y]]) swap(x,y); updata(1,n,1,p[top[x]],p[x]),x=fa[top[x]]; } if(x==y) return ; if(dep[x]>dep[y]) swap(x,y); updata(1,n,1,p[x]+1,p[y]); } int ask(int x,int y) { int ret=0; while(top[x]!=top[y]) { if(dep[top[x]]<dep[top[y]]) swap(x,y); ret+=query(1,n,1,p[top[x]],p[x]),x=fa[top[x]]; } if(x==y) return ret; if(dep[x]>dep[y]) swap(x,y); return ret+query(1,n,1,p[x]+1,p[y]); } int main() { n=rd(),m=rd(); int i,a,b; memset(head,-1,sizeof(head)); for(i=1;i<=m;i++) a=rd(),b=rd(),add(a,b); while(1) { op[++q]=rd(); if(op[q]==-1) break; qa[q]=rd(),qb[q]=rd(); if(!op[q]) { a=mp[MP(max(qa[q],qb[q]),min(qa[q],qb[q]))]; val[a]=val[a^1]=1; } } dep[1]=1,dfs1(1),dfs2(1,1),build(1,n,1); for(i=1;i<=pa[0];i++) modify(pa[i],pb[i]); for(i=q-1;i>=1;i--) { if(!op[i]) modify(qa[i],qb[i]); else ans[i]=ask(qa[i],qb[i]); } for(i=1;i<q;i++) if(op[i]) printf("%d\n",ans[i]); return 0; } //5 5 1 2 1 3 3 4 4 5 4 2 1 1 5 0 4 2 1 5 1 -1
【BZOJ1969】[Ahoi2005]LANE 航線規劃 離線+樹鏈剖分+線段樹