#include<stdio.h>
#include<queue>
#include<string.h>
using namespace std;
#define inf 0x3fffffff
#define N 200
struct node {
int v,w,next;
}bian[N*N*2],fbian[N*N*2];
int head[N],yong,tt;
int deep[N];//深度保留层次图
void addedge(int u,int v,int w) {
bian[yong].v=v;
bian[yong].w=w;
bian[yong].next=head[u];
head[u]=yong++;
}
int bfs(int s) {//寻找层次图
int i,cur;
queue<int>q;
memset(deep,-1,sizeof(deep));
q.push(s);
deep[s]=1;
while(!q.empty()) {
cur=q.front();
q.pop();
for(i=head[cur];i!=-1;i=bian[i].next)
if(bian[i].w&&deep[bian[i].v]==-1) {
deep[bian[i].v]=deep[cur]+1;
q.push(bian[i].v);
}
}
if(deep[tt]!=-1)return 1;
return 0;
}
int Min(int a,int b) {
return a>b?b:a;
}
int dfs(int s,int limit) {//多次dfs
int i,cost=0,flow,v;
// printf("%d %d\n",s,head[s]);
if(s==tt) {return limit;}
for(i=head[s];i!=-1;i=bian[i].next) {
v=bian[i].v;
// printf("%d %d %d\n",s,v,bian[i].w);
if(bian[i].w&&deep[v]==deep[s]+1) {
// printf("%d\n",v);
flow=dfs(v,Min(limit-cost,bian[i].w));
if(flow>0) {
bian[i].w-=flow;
bian[i^1].w+=flow;
cost+=flow;
if(limit==cost)break;
}
else
deep[v]=-1;//如果再次遇到v就不会进行下去
}
}
return cost ;
}
int dinic_maxflow(int s) {//dinic算法模板
int sum=0;
while(bfs(1)) {
//for(i=1;i<=tt;i++)
// printf("%d\n",deep[i]);
sum+=dfs(1,inf);
}
return sum;
}
/*int visit[N];
void dfs1(int u) {
int i;
visit[u]=1;
for(i=head[u];i!=-1;i=bian[i].next) {
if(bian[i].w&&visit[bian[i].v]==0)
dfs1(bian[i].v);
}
}*/
int ds[N],dds,dt[N],ddt;
int fhead[N],fyong;
int main() {
int t,m,n,i,j,k,maxflow,ma,fma,fa,fb;
scanf("%d",&t);
while(t--) {
scanf("%d%d",&n,&m);
yong=0;
memset(head,-1,sizeof(head));
while(m--) {
scanf("%d%d%d",&i,&j,&k);
addedge(i,j,k);//加入边
addedge(j,i,0);//退边
}
tt=n;
maxflow=dinic_maxflow(1);//最大流模板
dds=0;ddt=0;
/* memset(visit,0,sizeof(visit));
dfs1(1);*/
for(i=2;i<n;i++) {//根据deep是否为-1可以得到源点集和汇点集
if(deep[i]!=-1)//源点集
ds[++dds]=i;
else
dt[++ddt]=i;//汇点集
}
// printf("%d %d\n",dds,ddt);
/* for(i=2;i<n;i++) {
if(visit[i])
ds[++dds]=i;
else
dt[++ddt]=i;
}*/
for(i=0;i<yong;i++)//保留残留网络
fbian[i]=bian[i];
for(i=1;i<=n;i++)//保留残留网络
fhead[i]=head[i];
ma=0;
fyong=yong;//保留残留网络
for(i=1;i<=dds;i++)//枚举
for(j=1;j<=ddt;j++) {
fa=ds[i];fb=dt[j];
for(k=0;k<yong;k++)//加入了新边,必须跟新为原来的求最大流一次的残留网络
bian[k]=fbian[k];
memset(head,0,sizeof(head));
for(k=1;k<=n;k++) //因为加入了新边,所以必须跟新head的值为原来的残留网络
head[k]=fhead[k];
yong=fyong;//
addedge(fa,fb,inf);
addedge(fb,fa,0);
fma=dinic_maxflow(1);
if(fma>ma)//求出最大值
ma=fma;
}
printf("%d\n",ma+maxflow);
}
return 0;
}

hdu 2435 dinic算法模板+最小割性质的更多相关文章

  1. hdu 2435dinic算法模板+最小割性质

    hdu2435最大流最小割 2014-03-22 我来说两句 来源:hdu2435最大流最小割 收藏 我要投稿 2435 There is a war 题意: 给你一个有向图,其中可以有一条边是无敌的 ...

  2. POJ 1815 - Friendship - [拆点最大流求最小点割集][暴力枚举求升序割点] - [Dinic算法模板 - 邻接矩阵型]

    妖怪题目,做到现在:2017/8/19 - 1:41…… 不过想想还是值得的,至少邻接矩阵型的Dinic算法模板get√ 题目链接:http://poj.org/problem?id=1815 Tim ...

  3. POJ 1273 Drainage Ditches(网络流dinic算法模板)

    POJ 1273给出M条边,N个点,求源点1到汇点N的最大流量. 本文主要就是附上dinic的模板,供以后参考. #include <iostream> #include <stdi ...

  4. HDU1532最大流 Edmonds-Karp,Dinic算法 模板

    Drainage Ditches Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) To ...

  5. 洛谷.4897.[模板]最小割树(Dinic)

    题目链接 最小割树模板.具体见:https://www.cnblogs.com/SovietPower/p/9734013.html. ISAP不知为啥T成0分了.. Dinic: //1566ms ...

  6. POJ 3469.Dual Core CPU 最大流dinic算法模板

    Dual Core CPU Time Limit: 15000MS   Memory Limit: 131072K Total Submissions: 24830   Accepted: 10756 ...

  7. POJ 1273 - Drainage Ditches - [最大流模板题] - [EK算法模板][Dinic算法模板 - 邻接表型]

    题目链接:http://poj.org/problem?id=1273 Time Limit: 1000MS Memory Limit: 10000K Description Every time i ...

  8. HDU 3526 Computer Assembling(最小割)

    http://acm.hdu.edu.cn/showproblem.php?pid=3526 题意:有个屌丝要配置电脑,现在有n个配件需要购买,有两家公司出售这n个配件,还有m个条件是如果配件x和配件 ...

  9. [模板]最小割树(Gomory-Hu Tree)(luogu4897)

    给定一个\(n\)个点\(m\)条边的无向连通图,多次询问两点之间的最小割 两点间的最小割是这样定义的:原图的每条边有一个割断它的代价,你需要用最小的代价使得这两个点不连通 Input 第一行两个数\ ...

随机推荐

  1. 19 个必须知道的 VS 快捷键

    本文将为大家列出在Visual Studio中常用的快捷键,正确熟练地使用快捷键,将大大提高你的编程工作效率. 项目相关的快捷键 Ctrl + Shift + B = 生成项目 Ctrl + Alt ...

  2. Codeforces Round #253 (Div. 1)-A,B

    A题: 由题意可知,最多翻10次就能够(事实上8次就够了).那么我们就用状态压缩表示状态. 对于某种状态,假设某一位为0,那么代表这一位不翻,否则代表这一位翻. 对于某一种翻的状态: 假设牌中有G3, ...

  3. Foundation框架和文件操作

    NSString --实例化方法-------------- NSString *str = [[NSString alloc] init]; NSString *str = [[[NSString ...

  4. [iOS]UITableViewController完毕收回键盘操作

    UITableViewController 本身可以实现键盘适配(cell中控件焦点会移动到键盘上方 在做键盘收回的时候思考过例如以下方案 1.tableview加入点击事件 结果:点击事件和tabl ...

  5. iOS 获取当前时间格式化字符串

    iOS 获取当前时间格式化字符串 太阳火神的漂亮人生 (http://blog.csdn.net/opengl_es) 本文遵循"署名-非商业用途-保持一致"创作公用协议 转载请保 ...

  6. jstl 标签 循环 序号

     大家好: 今天搜了一天.最终找到它了,尽管不是我想要的,可是为了我辛苦的一天.我也要记录下: <c:forEach items="${signBusList}" var ...

  7. vim设置solarized配色

    最近在学习使用VIM.我选择Solarized配色.相信很多人也都在用. 官网地址: http://ethanschoonover.com/solarized 看这配色:八卦乾坤,赏心悦目,高端大气上 ...

  8. perl的安装

    http://www.activestate.com/activeperl/downloads 安装的时候,默认把perl放置到环境变量的PATH中 之后,需要重启电脑,确保环境变量生效 执行perl ...

  9. Asura监控---AsuraMonitor,阿修罗监控开源

    阿修罗Monitor是一个功能强大.灵活的监控系统. 系统安装简单,配置简单,相比zabbix, nagios,cacti,小米监控等都使用相当简单.只需要会写脚本,语言不限就可以实现任意监控需求. ...

  10. 一个php+jquery+json+ajax实例

    json.php <!DOCTYPE html Public "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://ww ...