POJ 1273 Drainage Ditches(最大流入门)
题目链接:poj1273
题意:
现在有m个池塘(从1到m开始编号,1为源点,m为汇点),及n条有向水渠,给出这n条水渠所连接的点和所能流过的最大流量,求从源点到汇点能流过的最大流量。
分析:
裸的求1点到M点最大流的问题,直接用Dinic算法即可.
AC代码:
#include<cstdio>
#include<cstring>
#include<queue>
#define INF 1e9
using namespace std;
const int maxn=200+20;//之前这里只写10+5,一直TLE,真是悲剧
struct Edge
{
Edge(){}
Edge(int from,int to,int cap,int flow):from(from),to(to),cap(cap),flow(flow){}
int from,to,cap,flow;
//cap是容量 flow是流量
};
struct Dinic
{
int n,m,s,t; //结点数,边数(包括反向弧),源点与汇点编号
vector<Edge> edges; //边表 edges[e]和edges[e^1]互为反向弧
vector<int> G[maxn]; //邻接表,G[i][j]表示结点i的第j条边在e数组中的序号
bool vis[maxn]; //BFS使用,标记一个节点是否被遍历过
int d[maxn]; //从起点到i点的距离
int cur[maxn]; //当前弧下标
void init(int n,int s,int t)
{
this->n=n,this->s=s,this->t=t;
for(int i=1;i<=n;i++) G[i].clear();
edges.clear();
}
void AddEdge(int from,int to,int cap)
{
edges.push_back( Edge(from,to,cap,0) );
edges.push_back( Edge(to,from,0,0) );//这条是反向边
m = edges.size();//总共边的数量
G[from].push_back(m-2);//正向边在edges数组中的编号
G[to].push_back(m-1);//反向边在edges数组中的编号
}
bool BFS()
{
memset(vis,0,sizeof(vis));
queue<int> Q;//用来保存节点编号的
Q.push(s);
d[s]=0;
vis[s]=true;
while(!Q.empty())
{
int x=Q.front(); Q.pop();
for(int i=0; i<G[x].size(); i++)//x所连的边
{
Edge& e=edges[G[x][i]];//取出所连的那条边
if(!vis[e.to] && e.cap>e.flow)//如果还没有经过 且流量小于容量 还可以再通
{
vis[e.to]=true;
d[e.to] = d[x]+1;
Q.push(e.to);
}
}
}
return vis[t];//表示能否到大t点 不能的话返回false
}
int DFS(int x,int a)
{
if(x==t || a==0)return a;//到大汇点或者流量为0时返回
int flow=0,f;//flow用来记录从x到t的最小残量
for(int& i=cur[x]; i<G[x].size(); i++)
{
Edge& e=edges[G[x][i]];
if(d[x]+1==d[e.to] && (f=DFS( e.to,min(a,e.cap-e.flow) ) )>0 )
//判断是否是可行边 并对其进行深搜 是否有结果
{
e.flow +=f;
edges[G[x][i]^1].flow -=f;//与其对应的反向边流量减少 以便于后悔恢复
flow += f;
a -= f;
if(a==0) break;
}
}
return flow;
}
int Maxflow()
{
int flow=0;
while(BFS())
{
memset(cur,0,sizeof(cur));
flow += DFS(s,INF);
}
return flow;
}
}DC;
int main()
{
int n,m;
while(~scanf("%d%d",&m,&n))
{
DC.init(n,1,n);
while(m--)
{
int u,v,w;
scanf("%d%d%d",&u,&v,&w);
DC.AddEdge(u,v,w);
}
printf("%d\n",DC.Maxflow());
}
return 0;
}