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

时间:2021-06-30 04:23:52


   给定一个有向网络,边权为拆掉边的代价,现在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;
}