Dinic邻接矩阵版本
2015-08-07 11:26
127 查看
#include<iostream>
#include<cstdio>
#include<string.h>
#include<string>
#include<stack>
#include<set>
#include<algorithm>
#include<cmath>
#include<vector>
#include<map>
#include<sstream>
#include<queue>
#define ll __int64
#define lll unsigned long long
#define MAX 1000009
#define MAXN 2009
#define eps 1e-8
#define INF 0x7fffffff
#define mod 1000000007
#define clr(a) memset(a,0,sizeof(a))
#define clr1(a) memset(a,-1,sizeof(a))
#define lson l , m , rt << 1
#define rson m + 1 , r , rt << 1 | 1
using namespace std;
inline ll Max(ll a,ll b)
{
return a>b?a:b;
}
inline ll Min(ll a,ll b)
{
return a<b?a:b;
}
//M为边数 N为点数 点标从1-n
const int N = 1000;
const int M = 100000;
int cap
;
int flow
;
int dis
,cur
;//距离起点的距离 cur[i]表示i点正在考虑的边 优化不再考虑已经用过的点 初始化为head
bool vis
;
bool BFS(int Start,int End)
{
memset(vis,0,sizeof(vis));
memset(dis,-1,sizeof(dis));
queue<int>Q;
while(!Q.empty())Q.pop();
Q.push(Start);
dis[Start]=0;
vis[Start]=1;
while(!Q.empty())
{
int u = Q.front();
Q.pop();
for(int v = 1; v<=N; v++)
{
if(!vis[v] && cap[u][v]>flow[u][v])
{
vis[v]=1;
dis[v]=dis[u]+1;
if(v==End)return true;
Q.push(v);
}
}
}
return false;
}
int DFS(int x, int a,int End) //流入x 的流量是a
{
if(x==End || a==0)return a;
int Flow = 0, f;
for(int v = 1; v<=N; v++)
{
if(dis[x]+1 == dis[v] && (f = DFS(v , Min(a,cap[x][v]-flow[x][v]), End))>0 )
{
flow[x][v] += f;
flow[v][x] -= f;//反向边要减掉
Flow += f;
a -= f;
if(a==0)break;
}
}
return Flow;
}
int Dinic(int Start,int End)
{
int Flow=0;
while(BFS(Start,End))
{
//cout<<Start<<" "<<End<<endl;
Flow += DFS(Start, INF , End);
}
return Flow;
}
int main()
{
#ifdef ONLINE_JUDGE
#else
freopen("an.txt","r", stdin);
#endif
int m,n;
int u,v,w;
while(~scanf("%d%d",&n,&m))
{
clr(cap);
clr(flow);
for(int i = 0;i<n;i++)
{
scanf("%d%d%d",&u,&v,&w);
cap[u][v]+=w;
}
cout<<Dinic(1,m)<<endl;
}
return 0;
}
#include<cstdio>
#include<string.h>
#include<string>
#include<stack>
#include<set>
#include<algorithm>
#include<cmath>
#include<vector>
#include<map>
#include<sstream>
#include<queue>
#define ll __int64
#define lll unsigned long long
#define MAX 1000009
#define MAXN 2009
#define eps 1e-8
#define INF 0x7fffffff
#define mod 1000000007
#define clr(a) memset(a,0,sizeof(a))
#define clr1(a) memset(a,-1,sizeof(a))
#define lson l , m , rt << 1
#define rson m + 1 , r , rt << 1 | 1
using namespace std;
inline ll Max(ll a,ll b)
{
return a>b?a:b;
}
inline ll Min(ll a,ll b)
{
return a<b?a:b;
}
//M为边数 N为点数 点标从1-n
const int N = 1000;
const int M = 100000;
int cap
;
int flow
;
int dis
,cur
;//距离起点的距离 cur[i]表示i点正在考虑的边 优化不再考虑已经用过的点 初始化为head
bool vis
;
bool BFS(int Start,int End)
{
memset(vis,0,sizeof(vis));
memset(dis,-1,sizeof(dis));
queue<int>Q;
while(!Q.empty())Q.pop();
Q.push(Start);
dis[Start]=0;
vis[Start]=1;
while(!Q.empty())
{
int u = Q.front();
Q.pop();
for(int v = 1; v<=N; v++)
{
if(!vis[v] && cap[u][v]>flow[u][v])
{
vis[v]=1;
dis[v]=dis[u]+1;
if(v==End)return true;
Q.push(v);
}
}
}
return false;
}
int DFS(int x, int a,int End) //流入x 的流量是a
{
if(x==End || a==0)return a;
int Flow = 0, f;
for(int v = 1; v<=N; v++)
{
if(dis[x]+1 == dis[v] && (f = DFS(v , Min(a,cap[x][v]-flow[x][v]), End))>0 )
{
flow[x][v] += f;
flow[v][x] -= f;//反向边要减掉
Flow += f;
a -= f;
if(a==0)break;
}
}
return Flow;
}
int Dinic(int Start,int End)
{
int Flow=0;
while(BFS(Start,End))
{
//cout<<Start<<" "<<End<<endl;
Flow += DFS(Start, INF , End);
}
return Flow;
}
int main()
{
#ifdef ONLINE_JUDGE
#else
freopen("an.txt","r", stdin);
#endif
int m,n;
int u,v,w;
while(~scanf("%d%d",&n,&m))
{
clr(cap);
clr(flow);
for(int i = 0;i<n;i++)
{
scanf("%d%d%d",&u,&v,&w);
cap[u][v]+=w;
}
cout<<Dinic(1,m)<<endl;
}
return 0;
}
相关文章推荐
- eclipse中常用快捷键
- oracle转Mysql中,varchar2(10)和number应该转换为什么类型?
- git笔记 常规使用
- Sort函数的相关知识
- 【暑假】[实用数据结构]UVAlive 3135 Argus
- 双缓冲解决VC++绘图时屏幕闪烁
- iOS 使用FMDB进行数据库操作
- Ubuntu 14.04 /etc/default/grub文件
- javax.el.PropertyNotFoundException: 'newsLine' is an unknown bean property...
- Tomcat 自动部署
- 41. Element getElementsByTagName() 方法
- 求Sn = a+aa+aaa+·······+a(n个a)a之值
- javax.el.PropertyNotFoundException: 'newsLine' is an unknown bean property...
- iOS 第六课 自己常用方法集合写法
- 隐藏文字
- Eclipse背景和匹配出现单词的一些设置
- 生活感受
- php生成随机密码的几种方法
- centOS 7.1.1503 搭建Apache+Mariadb+PHP环境
- 【dfs/bfs+set+快速幂】swjtuOJ 2094