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;
}