【BZOJ】【2594】【WC2006】水管局长

LCT

  动态维护MST嘛……但是有删边= =好像没法搞的样子

  离线记录所有修改&询问,倒序处理,就可以变删边为加边了~

  论如何用LCT维护最小生成树:先搞出一棵最小生成树,然后每次加边(u,v)时,在LCT上询问u->v这条链上权值最大的边,如果这条边权值比新加的边权值要小,则忽略这条新加的边,否则断掉这条权值最大的边(cut),加入这条新边(link)。

  然后……我不会捉……又去orz了一下Hzwer,学习了一下怎么维护边(拆边为点)……

  1 /**************************************************************
  2     Problem: 2594
  3     User: Tunix
  4     Language: C++
  5     Result: Accepted
  6     Time:21616 ms
  7     Memory:108244 kb
  8 ****************************************************************/
  9
 10 //BZOJ 2594
 11 #include<cstdio>
 12 #include<cstdlib>
 13 #include<cstring>
 14 #include<iostream>
 15 #include<algorithm>
 16 #define rep(i,n) for(int i=0;i<n;++i)
 17 #define F(i,j,n) for(int i=j;i<=n;++i)
 18 #define D(i,j,n) for(int i=j;i>=n;--i)
 19 using namespace std;
 20
 21 int getint(){
 22     int v=0,sign=1; char ch=getchar();
 23     while(ch<‘0‘||ch>‘9‘) {if (ch==‘-‘) sign=-1; ch=getchar();}
 24     while(ch>=‘0‘&&ch<=‘9‘) {v=v*10+ch-‘0‘; ch=getchar();}
 25     return v*=sign;
 26 }
 27 /*******************tamplate********************/
 28 const int N=1500005;
 29 struct LCT{
 30     int c[N][2],fa[N],v[N],size[N],mx[N];
 31     bool rev[N];
 32     int st[N],top;
 33     #define L c[x][0]
 34     #define R c[x][1]
 35     void Push_up(int x){
 36         mx[x]=x;
 37         if (v[mx[L]]>v[mx[x]]) mx[x]=mx[L];
 38         if (v[mx[R]]>v[mx[x]]) mx[x]=mx[R];
 39     }
 40     void Push_down(int x){
 41         if (rev[x]) rev[x]=0,rev[L]^=1,rev[R]^=1,swap(L,R);
 42     }
 43     bool not_root(int x){
 44         return c[fa[x]][0]==x || c[fa[x]][1]==x;
 45     }
 46     void rotate(int x){
 47         int y=fa[x],z=fa[y],l=c[y][1]==x,r=l^1;
 48         if (not_root(y)) c[z][c[z][1]==y]=x;
 49         fa[x]=z; fa[y]=x; fa[c[x][r]]=y;
 50         c[y][l]=c[x][r]; c[x][r]=y;
 51         Push_up(y);
 52     }
 53     void preview(int x){
 54         top=0; st[++top]=x;
 55         for(;not_root(x);x=fa[x])
 56             st[++top]=fa[x];
 57         D(i,top,1) Push_down(st[i]);
 58     }
 59     void splay(int x){
 60         int y=0;
 61         for(preview(x);not_root(x);rotate(x))
 62             not_root(y=fa[x]) ? rotate( (c[y][1]==x^c[fa[y]][1]==y ? x : y)),1 : 1;
 63         Push_up(x);
 64     }
 65     void access(int x,int las=0){
 66         for(;x;splay(x),c[x][1]=las,las=x,x=fa[x]);
 67     }
 68     void makeroot(int x){
 69         access(x),splay(x),rev[x]^=1;
 70     }
 71     void link(int x,int y){
 72         makeroot(x),fa[x]=y;
 73     }
 74     void cut(int x,int y){
 75         makeroot(x),access(y),splay(y);
 76         if (c[y][0]==x) c[y][0]=fa[x]=0;
 77     }
 78     int query(int x,int y){
 79         makeroot(x),access(y),splay(y);
 80         return mx[y];
 81     }
 82 }t;
 83 /*********************LCT***********************/
 84 int n,m,Q;
 85 struct edge{
 86     int x,y,w,id;
 87     bool d;
 88 }e[N];
 89 struct que{
 90     int x,y,k,id,ans;
 91 }q[N];
 92 bool cmp(edge a,edge b){ return a.w<b.w; }
 93 bool cmp2(edge a,edge b){ return a.x<b.x || (a.x==b.x && a.y<b.y);}
 94 bool cmp3(edge a,edge b){ return a.id<b.id; }
 95 int find(int x,int y){
 96     int l=1,r=m,mid;
 97     while(l<=r){
 98         mid=l+r>>1;
 99         if (e[mid].x<x || (e[mid].x==x && e[mid].y<y))l=mid+1;
100         else if(e[mid].x==x && e[mid].y==y) return mid;
101         else r=mid-1;
102     }
103 }
104 /********************edge&ques******************/
105 int fa[N];
106 int find(int x){ return fa[x]==x ? x : fa[x]=find(fa[x]); }
107 /********************并查集*********************/
108 int main(){
109 #ifndef ONLINE_JUDGE
110     freopen("2594.in","r",stdin);
111     freopen("2594.out","w",stdout);
112 #endif
113     n=getint(); m=getint(); Q=getint();
114     F(i,1,n) fa[i]=i;
115     F(i,1,m){
116         e[i].x=getint(); e[i].y=getint(); e[i].w=getint();
117         if (e[i].x>e[i].y) swap(e[i].x,e[i].y);
118     }
119     sort(e+1,e+m+1,cmp);//边权序
120     F(i,1,m){
121         e[i].id=i;
122         t.v[n+i]=e[i].w;
123         t.mx[n+i]=n+i;
124     }
125     sort(e+1,e+m+1,cmp2);//字典序
126     F(i,1,Q){
127         q[i].k=getint(); q[i].x=getint(); q[i].y=getint();
128         if (q[i].k==2){
129             if (q[i].x>q[i].y) swap(q[i].x,q[i].y);
130             int t=find(q[i].x,q[i].y);
131             e[t].d=1; q[i].id=e[t].id;
132             //找到每次删除的边在边权序中的位置
133         }
134     }
135     sort(e+1,e+m+1,cmp3);//边权序
136     int tot=0;
137     F(i,1,m)
138         if (!e[i].d){
139             int f1=find(e[i].x),f2=find(e[i].y);
140             if (f1!=f2){
141                 fa[f1]=f2;
142                 t.link(e[i].x,i+n); t.link(e[i].y,i+n);
143                 tot++;
144                 if (tot==n-1) break;
145             }
146         }
147     D(i,Q,1){
148         if (q[i].k==1)
149             q[i].ans=t.v[t.query(q[i].x,q[i].y)];
150         else{
151             int k=q[i].id;
152             int tmp=t.query(q[i].x,q[i].y);
153             if (e[k].w<t.v[tmp]){
154                 t.cut(e[tmp-n].x,tmp); t.cut(e[tmp-n].y,tmp);
155                 t.link(q[i].x,k+n); t.link(q[i].y,k+n);
156             }
157         }
158     }
159     F(i,1,Q) if(q[i].k==1)
160         printf("%d\n",q[i].ans);
161     return 0;
162 }

时间: 2024-10-06 05:37:07

【BZOJ】【2594】【WC2006】水管局长的相关文章

bzoj 2594: [Wc2006]水管局长数据加强版 动态树

2594: [Wc2006]水管局长数据加强版 Time Limit: 25 Sec  Memory Limit: 128 MBSubmit: 934  Solved: 291[Submit][Status] Description SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一条从A至B的水管的路径,接着通过信息化的控制中心通知路径上的水管进入准备送水状态,等到路径

BZOJ 2594: [Wc2006]水管局长数据加强版(kruskal + LCT)

Description SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一条从A至B的水管的路径,接着通过信息化的控制中心通知路径上的水管进入准备送水状态,等到路径上每一条水管都准备好了,供水公司就可以开始送水了.嘟嘟一次只能处理一项送水任务,等到当前的送水任务完成了,才能处理下一项. 在处理每项送水任务之前,路径上的水管都要进行一系列的准备操作,如清洗.消毒等等.嘟嘟在

BZOJ 2594 Wc2006 水管局长数据加强版 Link-Cut-Tree

题目大意:给定一个无向图,多次删除某条边,多次查询两点之间路径上边权最大值的最小值 Link-Cut-Tree维护动态最小生成树 首先倒着做 将所有被删除的边标记(找边我用的排序+二分) 将没标记的边跑一遍Kruskal 求出最小生成树 然后每次加边和查询正常维护即可 LInk-Cut-Tree一气呵成写完,Kruskal尼玛写挂了-- 居然忘记把并查集连边 这我也是醉了 顺便吐槽一下题干上给的读入优化真尼玛弱--自己随便写一个都可以优化到RANK前十-- #include<cstdio> #

BZOJ 2594 [Wc2006]水管局长数据加强版 LCT

题意:链接 方法: LCT 解析: 搞了一个上午加1个小时的题,TM最后大错误居然是排序元素太多排不回原来的样子! 我要重新学排序! 这题是用LCT维护动态最小生成树,但是最小生成树上删边应该是做不到的,所以我们可以离线操作,之后先把所有该删的边删了然后倒着搞所有询问,这样删边就变成了加边,之后询问就是x到y路径上的最大边权. 图是动态的,所以想到LCT,但是LCT不能搞最大边权怎么办! 把每个边看做一个点. 假设这是第i个边,那么把他看做第i+n个点. 显然点权就是边权,然后将这个边连接的两个

【BZOJ 2594】 [Wc2006]水管局长数据加强版

2594: [Wc2006]水管局长数据加强版 Time Limit: 25 Sec  Memory Limit: 128 MB Submit: 1138  Solved: 364 [Submit][Status] Description SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一条从A至B的水管的路径,接着通过信息化的控制中心通知路径上的水管进入准备送水状态,等

[bzoj2594][Wc2006]水管局长数据加强版

论蒟蒻的自我修养T_T.. 和noi2014魔法森林基本一样...然而数据范围大得sxbk...100w你告诉我(n+m)log(n+m)可过?[掀桌] 蒟蒻又蠢了..复杂度应该是O((n+q)log(n+m))吧.. 一开始数组开太小re了两发(要开到maxn+maxm),然后又开太大mle一发,然后无限tle...把记录类型全改成数组还是tle.... 最后把非lct部分改得和黄学长全部一样终于20+s卡过去了......... 然后发现自己原来是有个地方写萎了..一开始把没被删的边做kru

BZOJ_2594_[Wc2006]水管局长数据加强版_LCT

Description SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一条从A至B的水管的路径,接着通过信息化的控制中心通知路径上的水管进入准备送水状态,等到路径上每一条水管都准备好了,供水公司就可以开始送水了.嘟嘟一次只能处理一项送水任务,等到当前的送水任务完成了,才能处理下一项. 在处理每项送水任务之前,路径上的水管都要进行一系列的准备操作,如清洗.消毒等等.嘟嘟在

P4172 [WC2006]水管局长

[LCT]P4172水管局长 \(Solution\) 如果没有删除,那么就是维护一个最小生成树,然后倍增求两点之间的最大边权(货车运输). 因为有删边操作,想到LCT,但这是删除,最大值不满足减法,所以不好搞. 但注意到只有删除没有添加,所以我们可以倒过来处理,一条一条边link维护最小生成树以及两点间最大值,这样最大值就变成可加的了. 倒着处理,如果要加一条边u->v,那么先查询u->v的最大边权(假如这条边为x->y),如果要加的边权比这还大,那么就忽略(贪心),否则就删掉x-&g

[BZOJ2594] [Wc2006]水管局长数据加强版(LCT + kruskal + 离线)

传送门 WC这个题真是丧心病狂啊,就是想学习一下怎么处理边权,给我来了这么一个破题! ORZ hzwer 临摹黄学长代码233 但还是复杂的一匹 理一下思路吧 题目大意:给定一个无向图,多次删除图中的某一条边,求两点间路径最大值的最小值 求两点间的路径最大值的最小值的话,可以求最小生成树,那么这个值就是最小生成树上两点间路径上的最大值 但是题目要求是删除边,LCT维护最小生成树不支持删边操作,那么就离线处理,倒着加边,用LCT维护. 就是这个离线处理是最恶心的. 来说说如何处理边权,把边也抽象成

【bzoj2594】[Wc2006]水管局长数据加强版

真是神题 当时调了几天没调出来 后来没管了 当时把fread去掉就TLE,加上就RE 一直在底下跟网上的程序拍,尝试各种优化常数都没用 拍出几组不一样的,发现我是对的,醉了,网上那个是怎么过的 记一下这蛋疼的代码 1 #include<iostream> 2 #include<cstring> 3 #include<cstdlib> 4 #include<algorithm> 5 #include<cstdio> 6 #include<ma