hdu 2435(最小割 --dinic模板-- 最大流)

   给定一个有向网络，边权为拆掉边的代价，现在1要到n去，n试图阻止1到达，它至少花多大代价。    条件是1可以在任意两点（不含1和n）加入一条边（此边不可被拆除），求n要花费的最小代价最大值。    如果没条件，那么就是求最小割。其实最小割就是求最大流。    可以枚举符合条件的所有的边，但是消耗太大了~这是题目要卡的地方。    所以用求最小割的时候，把源集记录下来，剩下的就是汇集了，这也源集是DINIC算法中BFS中找不到汇点时的queue中的点。然后枚举源集到汇集所有的边就可以了，，算出附加代价。   这题可以作为DINIC的模板题。#include<iostream>#include<cstdio>#include<cstring>#include<set>#include<cstdio>#include<string>using namespace std;const int maxn=105;//const int maxm=300006;const int inf=1<<30;struct node{    int v,next;    int val;}s[maxn*maxn*2];int level[maxn];//顶点的层次int p[maxn];int que[maxn*10];//BFS中用于遍历的顶点，DFS求增广中记录边int out[10*maxn];//DFS用于几乎定点的分支int ind;int cop_poit[maxn];node cop_mp[maxn*maxn*2];int topset;void init(){    ind=0;    memset(p,-1,sizeof(p));}inline void insert(int x,int y,int z){    s[ind].v=y;    s[ind].val=z;    s[ind].next=p[x];    p[x]=ind++;    s[ind].v=x;    s[ind].val=0;    s[ind].next=p[y];    p[y]=ind++;}int max_flow(int n,int source,int sink){    int ret=0;    int h=0,r=0;    while(1)//DFS    {        int i;        for(i=0;i<=n;++i)        level[i]=0;        h=0,r=0;        level[source]=1;        que[0]=source;        while(h<=r)//BFS        {            int  t=que[h++];            for(i=p[t];i!=-1;i=s[i].next)            {                if(s[i].val&&level[s[i].v]==0)                {                    level[s[i].v]=level[t]+1;                    que[++r]=s[i].v;                }            }        }        topset=r;//记录原点的集合个数        if(level[sink]==0)break;//找不到汇点        for(i=0;i<=n;++i)        out[i]=p[i];        int  q=-1;        while(1)        {            if(q<0)            {                int  cur=out[source];                for(;cur!=-1;cur=s[cur].next)                {                    if(s[cur].val&&out[s[cur].v]!=-1&&level[s[cur].v]==2)                    {                        break;                    }                }                if(cur>=0)                {                    que[++q]=cur;                    out[source]=s[cur].next;                }                else                {                    break;                }            }            int  u=s[que[q]].v;            if(u==sink)//一条增广路            {                int  dd=inf;                int  index=-1;                for(i=0;i<=q;i++)                {                    if(dd>s[que[i]].val)                    {                        dd=s[que[i]].val;                        index=i;                    }                }                ret+=dd;                //cout<<ret<<endl;                for(i=0;i<=q;i++)                {                    s[que[i]].val-=dd;                    s[que[i]^1].val+=dd;                }                for(i=0;i<=q;i++)                {                    if(s[que[i]].val==0)                    {                        q=index-1;                        break;                    }                }            }            else            {                long cur=out[u];                for(;cur!=-1;cur=s[cur].next)                {                    if(s[cur].val&&out[s[cur].v]!=-1&&level[u]+1==level[s[cur].v])                    {                        break;                    }                }                if(cur!=-1)                {                    que[++q]=cur;                    out[u]=s[cur].next;                }                else                {                    out[u]=-1;                    q--;                }            }        }    }    return ret;}int m,n;int main(){    int cas;    cin>>cas;    set<int>ss;    int huijie[maxn];    while(cas--)    {        init();        scanf("%d%d",&n,&m);        for(int i=0;i<m;i++)        {            long from,to,cost;            scanf("%d %d %d",&from,&to,&cost);            insert(from,to,cost);        }        long Start,End;        Start=1;End=n;        int tmp=max_flow(n,Start,End);//计算一次最大流（最小割）        for(int i=1;i<=n;i++)//复制网络        cop_poit[i]=p[i];        for(int i=0;i<ind;i++)        cop_mp[i]=s[i];        int ans=0;        int r=0;        for(int j=2;j<n;j++)//求出汇集---因为bfs找不到一条增广路时，qu中就是原集，剩下来就是汇集，也可以是leve数组便利不到的顶点       {            if(level[j]==0)            huijie[r++]=j;        }        for(int i=0;i<=topset;i++)//枚举两顶点        {            if(que[i]!=1)            {                for(int j=0;j<r;j++)                {                    for(int k=1;k<=n;k++)                    p[k]=cop_poit[k];                    for(int k=0;k<ind;k++)                    s[k]=cop_mp[k];                    insert(que[i],huijie[j],inf);//加边                    int g=max_flow(n,1,n);//格外的代价                    if(ans<g)ans=g;                    //删边                    p[que[i]]=s[ind-2].next;                    p[huijie[j]]=s[ind-1].next;                    ind-=2;                }            }        }        printf("%d\n",ans+tmp);    }    return 0;}