您的位置：首页 > 其它

poj 1459 最大流 Dinic模板题

2016-07-09 18:15 387 查看

题意：

给几个发电站，给几个消耗站，再给几个转发点。

发电站只发电，消耗站只消耗电，转发点只是转发电，再给各个传送线的传电能力。

问你消耗站能获得的最多电是多少。

思路：这题建图是比较显然的，超源点和发电战之间连边，容量是发电量，消耗站和超汇点之间连边，容量是消耗量，然后输电线本来就是边，这样就建好图了，剩下的就是套最大流的模板了。

此题n<=100，m<=n^2，边比较多，所以采用DInic算法，时间复杂度O(mn^2)。如果用EK算法，并且输入用cin的话可能会超时，这题用scanf会节省很多时间。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
using namespace std;
const int INF=0x3f3f3f3f;
const int N=100+9;
struct Edge
{
int from,to,cap,flow;
};
vector<Edge>edges;
vector<int>G
;
int d
; //从起点到i点的距离
int cur
; //当前弧下标
bool vis
;
int s,t;
bool bfs()
{
memset(vis,0,sizeof(vis));
queue<int>q;
q.push(s);
d[s]=0;
vis[s]=1;
while(!q.empty()){
int x=q.front();q.pop();
for(int i=0;i<G[x].size();i++){
Edge& e=edges[G[x][i]];
if(!vis[e.to]&&e.cap>e.flow){
vis[e.to]=1;
d[e.to]=d[x]+1;
q.push(e.to);
}
}
}
return vis[t];
}
int dfs(int x,int a)
{
if(x==t||a==0)return a;
int flow=0,f;
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;
}
void addedge(int from,int to,int cap)
{
edges.push_back((Edge){from,to,cap,0});
edges.push_back((Edge){to,from,0,0});
int m=edges.size();
G[from].push_back(m-2);
G[to].push_back(m-1);
}

int main()
{
int n,np,nc,m;
while (~scanf("%d%d%d%d",&n,&np,&nc,&m)){
s=0,t=n+1;
int a, b, c;
for (int i = 0; i <= n+1; i++)
G[i].clear();
edges.clear();
for (int i = 0; i < m; i++) {
scanf(" (%d,%d)%d",&a,&b,&c);
//cin >> t >> a >> t >> b >> t >> c;
addedge(a + 1, b + 1, c);
}
for (int i = 0; i < np; i++) {
//cin >> t >> a >> t >> b;
scanf(" (%d)%d",&a,&b);
addedge(s, a + 1, b);
}
for (int i = 0; i < nc; i++) {
//cin >> t >> a >> t >> b;
scanf(" (%d)%d",&a,&b);
addedge(a + 1, t, b);
}
printf("%d\n", Maxflow());
}
return 0;
}

内容来自用户分享和网络整理，不保证内容的准确性，如有侵权内容，可联系管理员处理

标签：

相关文章推荐

新的分享

章节导航