[BZOJ]2594 水管局长数据加强版(Wc2006)
失踪人口回归。
LCT一直是小C的弱项,特别是这种维护链的信息的,写挂了就会调代码调到心态爆炸。
不过还好这一次的模板练习没有出现太多的意外。
Description
SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一条从A至B的水管的路径,接着通过信息化的控制中心通知路径上的水管进入准备送水状态,等到路径上每一条水管都准备好了,供水公司就可以开始送水了。嘟嘟一次只能处理一项送水任务,等到当前的送水任务完成了,才能处理下一项。
在处理每项送水任务之前,路径上的水管都要进行一系列的准备操作,如清洗、消毒等等。嘟嘟在控制中心一声令下,这些水管的准备操作同时开始,但由于各条管道的长度、内径不同,进行准备操作需要的时间可能不同。供水公司总是希望嘟嘟能找到这样一条送水路径,路径上的所有管道全都准备就绪所需要的时间尽量短。嘟嘟希望你能帮助他完成这样的一个选择路径的系统,以满足供水公司的要求。另外,由于MY市的水管年代久远,一些水管会不时出现故障导致不能使用,你的程序必须考虑到这一点。
不妨将MY市的水管网络看作一幅简单无向图(即没有自环或重边):水管是图中的边,水管的连接处为图中的结点。
Input
输入文件第一行为3个整数:N, M, Q分别表示管道连接处(结点)的数目、目前水管(无向边)的数目,以及你的程序需要处理的任务数目(包括寻找一条满足要求的路径和接受某条水管坏掉的事实)。
以下M行,每行3个整数x, y和t,描述一条对应的水管。x和y表示水管两端结点的编号,t表示准备送水所需要的时间。我们不妨为结点从1至N编号,这样所有的x和y都在范围[1, N]内。
以下Q行,每行描述一项任务。其中第一个整数为k:若k=1则后跟两个整数A和B,表示你需要为供水公司寻找一条满足要求的从A到B的水管路径;若k=2,则后跟两个整数x和y,表示直接连接x和y的水管宣布报废(保证合法,即在此之前直接连接x和y尚未报废的水管一定存在)。
Output
按顺序对应输入文件中每一项k=1的任务,你需要输出一个数字和一个回车/换行符。该数字表示:你寻找到的水管路径中所有管道全都完成准备工作所需要的时间(当然要求最短)。
Sample Input
4 4 3
1 2 2
2 3 3
3 4 2
1 4 2
1 1 4
2 1 4
1 1 4
Sample Output
2
3
HINT
N ≤ 100000
M ≤ 1000000
Q ≤ 100000
任何时候我们考虑的水管网络都是连通的,即从任一结点A必有至少一条水管路径通往任一结点B。
Solution
很容易看出题目要求我们维护无向图的最小生成树,查询两点树上路径最大权值,并支持删边操作。
删边不好处理,而且题目没有要求我们在线,我们很容易想到倒过来做,把删边改为加边。
于是用hash判断删去了哪些边(注意hash可能会撞),对删去所有边之后的图做最小生成树。
这样就是LCT维护最小生成树的裸题了。
对于每条新加入的边,查找两端点之间路径最大权值,如果小于那个权值,就删去那条边,并加入这条边。
#include <cstdio> #include <cstring> #include <algorithm> #define l(a) son[a][0] #define r(a) son[a][1] #define nrt(a) (l(fa[a])==a || r(fa[a])==a) #define MM 1000005 #define MN 100005 #define mod 10000007 using namespace std; struct meg1{int x,y,z;}a[MM]; struct meg2{int g,x,y,z,pos;}b[MN]; struct edge{int nex,to,wt,pos;}e[MN<<1]; struct dmp{dmp *nex; int x,y,val;}*mp[mod]; int son[MN][2],fa[MN],tfa[MN],tson[MN],tg[MN],mx[MN]; int hr[MN],hua[MN],pin; int n,m,p; inline int read() { int n=0,f=1; char c=getchar(); while (c<'0' || c>'9') {if(c=='-')f=-1; c=getchar();} while (c>='0' && c<='9') {n=n*10+c-'0'; c=getchar();} return n*f; } void revs(int x) {tg[x]^=1; swap(l(x),r(x)); swap(tfa[x],tson[x]);} void down(int x) {if (tg[x]) {revs(l(x)); revs(r(x)); tg[x]=0;}} int L(int x) {down(x); return l(x)?L(l(x)):x;} void aldown(int x) {if (nrt(x)) aldown(fa[x]); down(x);} void upd(int &x,int y) {if (a[y].z>a[x].z) x=y;} void update(int x) {mx[x]=0; upd(mx[x],tfa[x]); upd(mx[x],tson[x]); upd(mx[x],mx[l(x)]); upd(mx[x],mx[r(x)]);} void rotate(int x) { register int y,z,l,r; y=fa[x]; z=fa[y]; l=r(y)==x; r=l^1; if (nrt(y)) son[z][r(z)==y]=x; fa[y]=x; fa[x]=z; fa[son[x][r]]=y; son[y][l]=son[x][r]; son[x][r]=y; update(y); } void splay(int x) { aldown(x); register int y,z; for (;nrt(x);rotate(x)) {y=fa[x]; if (nrt(y)) z=fa[y],rotate(r(z)==y^r(y)==x?x:y);} update(x); } void access(int x) { for (int i=0;x;x=fa[i=x]) { splay(x); tson[x]=tfa[L(i)]; r(x)=i; update(x); } } int getmx(int x,int y) { access(x); splay(x); revs(x); access(y); splay(y); return mx[y]; } void link(int x,int y,int z) { access(x); splay(x); access(y); splay(y); revs(y); fa[y]=x; r(x)=y; tfa[y]=tson[x]=z; update(y); update(x); } void cut(int x,int y) { access(x); splay(x); revs(x); access(y); splay(y); l(y)=fa[l(y)]=0; tson[x]=tfa[y]=0; update(x); update(y); } inline void ins(int x,int y,int z,int w) { e[++pin]=(edge){hr[x],y,z,w}; hr[x]=pin; e[++pin]=(edge){hr[y],x,z,w}; hr[y]=pin; } void dfs(int x,int fat) { fa[x]=fat; for (register int i=hr[x];i;i=e[i].nex) if (e[i].to!=fat) mx[e[i].to]=tfa[e[i].to]=e[i].pos,dfs(e[i].to,x); } inline bool gather(int x,int y) { if (x==y) return false; hua[x]=y; return true; } int getf(int x) {return hua[x]?hua[x]=getf(hua[x]):x;} bool cmp(const meg1& a,const meg1& b) {return a.z<b.z;} int seamp(int h,int x,int y) { for (dmp* q=mp[h];q;q=q->nex) if (q->x==x && q->y==y) return q->val; return 0; } void insmp(int h,int x,int y,int z) { dmp* q=new(dmp); q->nex=mp[h]; q->x=x; q->y=y; q->val=z; mp[h]=q; } int geth(int x,int y) {return (1LL*x*300007+y)%mod;} int main() { register int i,lt; n=read(); m=read(); p=read(); for (i=1;i<=m;++i) { a[i].x=read(); a[i].y=read(); a[i].z=read(); if (a[i].x>a[i].y) swap(a[i].x,a[i].y); } for (i=1;i<=p;++i) { b[i].g=read(); b[i].x=read(); b[i].y=read(); if (b[i].x>b[i].y) swap(b[i].x,b[i].y); if (b[i].g==2) insmp(geth(b[i].x,b[i].y),b[i].x,b[i].y,i); } sort(a+1,a+m+1,cmp); for (i=1;i<=m;++i) { if (lt=seamp(geth(a[i].x,a[i].y),a[i].x,a[i].y)) {b[lt].pos=i; continue;} if (gather(getf(a[i].x),getf(a[i].y))) ins(a[i].x,a[i].y,a[i].z,i); } dfs(1,0); for (i=p;i;--i) { if (b[i].g==1) b[i].z=a[getmx(b[i].x,b[i].y)].z; else if (b[i].g==2) {if (a[b[i].pos].z<a[lt=getmx(b[i].x,b[i].y)].z) cut(a[lt].x,a[lt].y),link(b[i].x,b[i].y,b[i].pos);} } for (i=1;i<=p;++i) if (b[i].g==1) printf("%d\n",b[i].z); }
Last Word
由于link和cut部分可能有点难以理解,LCT一直是令小C头疼的算法之一,但是小C还是成功把模板默出来了。
然后迷之RE,并且发现树上部分节点的关系违背了伦理道德。调完hash之后,结果发现是在找最左边的节点的节点时忘记down操作了。
果然LCT在维护信息时,考虑哪里需要update,哪里需要down,顺序是什么是最烦的。
不得不说,D-map真的好用。