日……又被傻B错坑了一整天……
原题:
SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一条从A至B的水管的路径,接着通过信息化的控制中心通知路径上的水管进入准备送水状态,等到路径上每一条水管都准备好了,供水公司就可以开始送水了。嘟嘟一次只能处理一项送水任务,等到当前的送水任务完成了,才能处理下一项。
在处理每项送水任务之前,路径上的水管都要进行一系列的准备操作,如清洗、消毒等等。嘟嘟在控制中心一声令下,这些水管的准备操作同时开始,但由于各条管道的长度、内径不同,进行准备操作需要的时间可能不同。供水公司总是希望嘟嘟能找到这样一条送水路径,路径上的所有管道全都准备就绪所需要的时间尽量短。嘟嘟希望你能帮助他完成这样的一个选择路径的系统,以满足供水公司的要求。另外,由于MY市的水管年代久远,一些水管会不时出现故障导致不能使用,你的程序必须考虑到这一点。
不妨将MY市的水管网络看作一幅简单无向图(即没有自环或重边):水管是图中的边,水管的连接处为图中的结点。
N ≤ 100000
M ≤ 1000000
Q ≤ 100000
M ≤ 1000000
Q ≤ 100000
删边咋都不好维护,离线边加边
求的是路径上最大值最小,这回不二分了,在最小生成树上亮点路径中最大边就是
然后开始kruskal,加边lct就行了。。。
题好写,但是被两个傻B错坑了整整一天
1.
之前的3/4天我一直在debug,直到我看到了一篇题解:
然后我tm才发现删过边后图可能不连通。。。。。。
2.
不解释了,心累
代码:
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
int rd(){int z=; char ch=getchar();
while(ch<''||ch>'') ch=getchar();
while(ch>=''&&ch<=''){z=(z<<)+(z<<)+ch-''; ch=getchar();}
return z;
}
struct nds{int x,y,z,id;}e[],q[];
int n,m,o;
int fth[],chd[][],mxv[],v[],mxid[],rvs[];
int stck[],tp=;
bool flg[]; int tl[],tr[];
int ans[],ast=;
inline bool isrt(int x){ return (chd[fth[x]][]!=x)&(chd[fth[x]][]!=x);}
inline void pshu(int x){
if(!x) return ;
mxv[x]=v[x],mxid[x]=x;
if(mxv[chd[x][]]>mxv[x]) mxv[x]=mxv[chd[x][]],mxid[x]=mxid[chd[x][]];
if(mxv[chd[x][]]>mxv[x]) mxv[x]=mxv[chd[x][]],mxid[x]=mxid[chd[x][]];
}
inline void pshd(int x){
if(!rvs[x]) return ;
rvs[chd[x][]]^=,rvs[chd[x][]]^=,rvs[x]=;
swap(chd[x][],chd[x][]);
}
void rtt(int x){
int y=fth[x],z=fth[fth[x]],l,r;
r=(chd[y][]==x); l=r^;
if(!isrt(y)) chd[z][chd[z][]==y]=x;
fth[x]=z,fth[y]=x,fth[chd[x][r]]=y;
chd[y][l]=chd[x][r],chd[x][r]=y;
pshu(y),pshu(x);
}
void sply(int x){
stck[tp=]=x;
for(int i=x;!isrt(i);i=fth[i]) stck[++tp]=fth[i];
while(tp) pshd(stck[tp--]);
while(!isrt(x)){
if(!isrt(fth[x])) rtt((chd[fth[x]][]==x)^(chd[fth[fth[x]]][]==fth[x])?x:fth[x]);
rtt(x);
}
}
inline void accs(int x){ for(int i=;x;sply(x),chd[x][]=i,pshu(x),x=fth[i=x]);}
inline void qdrt(int x){ accs(x),sply(x),rvs[x]^=;}
inline void lk(int x,int y){ qdrt(x),fth[x]=y,sply(x);}
inline void ct(int x,int y){ qdrt(x),accs(y),sply(y),fth[x]=chd[y][]=,pshu(y);}
inline int gtrt(int x){ while(fth[x]) x=fth[x]; return x;}
inline int qry(int x,int y){ qdrt(x),accs(y),sply(y);
return mxid[y]-n;}
bool cmp1(nds x,nds y){ return x.x==y.x ? x.y<y.y : x.x<y.x;}
bool cmp2(nds x,nds y){ return x.z<y.z;}
bool cmp3(nds x,nds y){ return x.id<y.id;}
int bnrsch(int x,int y){
int l=tl[x],r=tr[x],md;
while(l+<r) md=(l+r)>>,(e[md].y<=y?l:r)=md;
//while(l+1<r) md=(l+r)>>1,(e[md].y>=y?l:r)=md;
return e[l].y==y?l:r;
}
void krsk(){
int bwl=;
sort(e+,e+m+,cmp2);
for(int i=;i<=m;++i)if((!flg[e[i].id])&(gtrt(e[i].x)!=gtrt(e[i].y))){
lk(e[i].x,e[i].id+n),lk(e[i].y,e[i].id+n);
if(++bwl==n-) break;
}
}
void ist(int x){
if(gtrt(e[x].x)!=gtrt(e[x].y)){
lk(e[x].x,x+n),lk(e[x].y,x+n);
return ;
}
int tmp=qry(e[x].x,e[x].y);
if(e[tmp].z<e[x].z) return ;
ct(e[tmp].x,tmp+n),ct(e[tmp].y,tmp+n);
mxv[tmp+n]=v[tmp+n]=;
lk(e[x].x,x+n),lk(e[x].y,x+n);
}
int main(){//freopen("ddd.in","r",stdin);
//freopen("ddd.out","w",stdout);
cin>>n>>m>>o;
for(int i=;i<=m;++i){
e[i].x=rd(),e[i].y=rd(),e[i].z=rd();
e[i].id=i,mxv[i+n]=v[i+n]=e[i].z,mxid[i+n]=i+n;
}
sort(e+,e+m+,cmp1);
for(int i=;i<=m+;++i)if(e[i].x!=e[i-].x) tr[e[i-].x]=i-,tl[e[i].x]=i;
for(int i=;i<=o;++i){
q[i].z=rd(),q[i].x=rd(),q[i].y=rd();
if(q[i].z==) flg[q[i].id=e[bnrsch(q[i].x,q[i].y)].id]=;
}
krsk();
sort(e+,e+m+,cmp3);
for(int i=o;i>=;--i){
if(q[i].z==) ist(q[i].id);
else ans[++ast]=e[qry(q[i].x,q[i].y)].z;
}
while(ast) printf("%d\n",ans[ast--]);
return ;
}