ccz181078

  博客园 :: 首页 :: 博问 :: 闪存 :: 新随笔 :: 联系 :: :: 管理 ::

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 *++ptr
const 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;
}

 

posted on 2017-03-23 20:51  nul  阅读(266)  评论(0编辑  收藏  举报