首页 > 代码库 > bzoj4331: JSOI2012 越狱老虎桥

bzoj4331: JSOI2012 越狱老虎桥

Description

这里,是美丽的南京;这里,是秀美的进香河;这里是安逸的老虎桥。 
如果说进香河的美,美在其秀美的风光,倒不如说是美在了那惬意的南京古典小巷式生活。如果说进香河的迷人,在其淳朴的民风,倒不如说是那被历史掩埋了的秘密吸引着人们好奇的心。 
也许很多人都还记得,老虎桥监狱,北洋时期江南最大的监狱,在近一个世纪中,面对满清、北洋、民国、新中国几朝兴衰,名称屡次更替,沧桑尽显其中。 
现在的人们,恐怕很难相信,到底有多少惊心动魄的事情曾经就在这里上演。 
那是1948年的冬天,南京地下组织的一支小分队决定偷袭老虎桥监狱,救出被困的数百名人员。那时的老虎桥监狱,被N层电网围了起来,由内而外,依次编号为1,2,…,N。第1层电网接有高压电。有M条高压线,连接了所有电网,其中第i条高压线连接了第ai和bi层电网,如果要破坏第i条高压线,需要至少动用Ti位特工。面对这么多层电网,偷袭小分队犯愁了。至少需要破坏一层电网,否则是无法偷袭成功的。 
然而,狡猾的间谍却知道了这件事情,为了破坏偷袭计划,敌人秘密地又增加了一条高压线,不让偷袭小分队的成员发现。 
为了能够偷袭成功,不论新增的这一条秘密高压线是连接哪两层电网的,小分队都必须要破坏且仅破坏一条高压线,使得至少有一层电网不通电。注意,对于新增的高压线,我们并不知道需要多少位特工才能成功破坏。现在的问题是,偷袭小分队至少需要多少名特工呢? 
决战就在今夜! 

Input

第一行有2个整数,N和M,分别表示电网层数和高压线个数。 
之后M行,每行3个整数,分别是ai, bi和Ti。 

Output

输出只有一行,包含一个整数,表示最少需要动用的特工人数。 
如果计划必然失败,则输出 -1。 
边双联通分量内的边删了不影响联通情况,所以先缩点,留下割边,问题转化为在树上选一条路径,使未被覆盖的边的最小值最大,可以直接树形dp
#include<cstdio>#define G *++ptrconst int N=500007,inf=0x3f3f3f3f;char buf[N*60],*ptr=buf-1;int _(){    int x=0,c=G;    while(c<48)c=G;    while(c>47)x=x*10+c-48,c=G;    return x;}bool ei[N*4];int n,m,es[N*6],enx[N*6],ev[N*6],e0[N],e1[N],ep=2,id[N],idp;void ae(int*e,int a,int b,int c){    es[ep]=b;enx[ep]=e[a];ev[ep]=c;e[a]=ep++;    es[ep]=a;enx[ep]=e[b];ev[ep]=c;e[b]=ep++;}int dfn[N],low[N],tk=0;void mins(int&a,int b){if(a>b)a=b;}void maxs(int&a,int b){if(a<b)a=b;}int min(int a,int b){return a<b?a:b;}int max(int a,int b){return a>b?a:b;}void tj(int w){    dfn[w]=low[w]=++tk;    for(int i=e0[w];i;i=enx[i]){        int u=es[i];        if(!u)continue;        if(!dfn[u]){            es[i^1]=0;            tj(u);            es[i^1]=w;            mins(low[w],low[u]);            if(low[u]>dfn[w])ei[i>>1]=1;        }else mins(low[w],dfn[u]);    }}void f1(int w){    id[w]=idp;    dfn[w]=0;    for(int i=e0[w];i;i=enx[i]){        int u=es[i];        if(!ei[i>>1]&&dfn[u])f1(u);    }}int v0[N],v02[N],v1[N],v2[N],vu[N],v,ans=0;void f2(int w,int pa){    v0[w]=v02[w]=v1[w]=v2[w]=vu[w]=inf;    for(int i=e1[w];i;i=enx[i]){        int u=es[i];        if(u==pa)continue;        f2(u,w);        v=min(v0[u],ev[i]);        v2[w]=max(min(v2[w],v),min(v1[w],v1[u]));        v1[w]=max(min(v1[w],v),min(v0[w],v1[u]));        if(v<=v0[w])v02[w]=v0[w],v0[w]=v;        else mins(v02[w],v);    }}void f3(int w,int pa){    v=min(vu[w],v2[w]);    maxs(ans,v);    for(int i=e1[w];i;i=enx[i]){        int u=es[i];        if(u==pa)continue;        vu[u]=min(min(ev[i],vu[w]),min(v0[u],ev[i])==v0[w]?v02[w]:v0[w]);        f3(u,w);    }}int main(){    fread(buf,1,sizeof(buf),stdin)[buf]=0;    n=_();m=_();    for(int i=0,a,b,c;i<m;++i){        a=_();b=_();c=_();        if(a==b)continue;        ae(e0,a,b,c);    }    tj(1);for(int i=1;i<=n;++i)if(dfn[i]){        ++idp;        f1(i);    }    for(int i=2;i<ep;i+=2)if(ei[i>>1])ae(e1,id[es[i]],id[es[i^1]],ev[i]);    f2(1,0);    f3(1,0);    if(ans==inf)ans=-1;    printf("%d",ans);    return 0;}

 

bzoj4331: JSOI2012 越狱老虎桥