poj 1459（最大流入门）

题意：
在一个网络图中有n个点,其中有np个发电站,nc个用户,m条电线;每个发电站,用户,和电线都对应有一个最大的电流;让求出该网络中最大的电流;

添加一个源点为(0)和一个汇点(n+1),源点与发电站相连,电流量为发电站的最大电流;用户与汇点相连,电流量为用户需要的最大电流。

EK算法模版：

10944604

xinghan0219

1459

Accepted

768K

625MS

G++

1290B

2012-10-22 23:17:50

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>

using namespace std;

#define inf 0x7ffffff
#define Min(x,y) (x>y?y:x)
#define MAXN 105
int n,np,nc,m;
int flow[MAXN][MAXN];
int maxflow;
int st,en;
int a[MAXN];
int cap[MAXN][MAXN];
int pre[MAXN];

int EK()
{
    queue<int>q;
    while(!q.empty()) q.pop();
    memset(flow,0,sizeof(flow));
    maxflow=0;
    while(1)
    {
        memset(a,0,sizeof(a));
        a[st]=inf;
        q.push(st);
        while(!q.empty())
        {
            int u=q.front();q.pop();
            for(int v=0;v<=n+1;v++)
                if(!a[v] && cap[u][v]>flow[u][v])
                {
                    pre[v]=u;q.push(v);
                    a[v]=Min(a[u],cap[u][v]-flow[u][v]);
                }
        }
        if(a[en]==0) break;
        for(int u=en;u!=st;u=pre[u])
        {
            flow[pre[u]][u]+=a[en];
            flow[u][pre[u]]-=a[en];
        }
        maxflow+=a[en];
    }
    return maxflow;
}

int main()
{
    while(scanf("%d%d%d%d",&n,&np,&nc,&m) != EOF)
    {
        int u,v,c;
        memset(cap,0,sizeof(cap));
        for(int i=0;i<m;i++)
        {
            scanf(" (%d,%d)%d",&u,&v,&c);
            cap[u+1][v+1]+=c;
        }
        for(int i=0;i<np;i++)
        {
            scanf(" (%d)%d",&v,&c);
            cap[0][v+1]+=c;
        }
        for(int i=0;i<nc;i++)
        {
            scanf(" (%d)%d",&u,&c);
            cap[u+1][n+1]+=c;
        }
        st=0,en=n+1;
        printf("%d\n",EK());
    }
    return 0;
}

dinic算法模版：

10945264

xinghan0219

1459

Accepted

720K

422MS

G++

1669B

2012-10-23 10:11:55

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>

using namespace std;

#define MAXN 105
#define inf 0x7ffffff
#define Min(x,y) (x>y?y:x)
int cap[MAXN][MAXN];//记录图的邻接矩阵
int dis[MAXN];//记录图中各点的层次
int n,m;//点,边
int st,en;//源点，汇点

bool bfs()
{
    queue<int>q;
    while(!q.empty()) q.pop();
    memset(dis,-1,sizeof(dis));
    dis[st]=0;//如果处理不慎就很容易死循环
    q.push(st);
    while(!q.empty())
    {
        int u=q.front();q.pop();
        for(int v=0;v<=n+1;v++)//注意此处是所有顶点(0,n+1是超级源点,汇点)
            if(dis[v]==-1 && cap[u][v])//此处应是cap[v][i]!=0,因为dfs是可能会走错，这样可以对其进行修正。
            {
                dis[v]=dis[u]+1;
                q.push(v);
            }
    }
    return dis[en]!=-1;//汇点的层次不等于-1说明找到
}

int dfs(int u,int cur_flow)
{
    int dt=cur_flow;
    if(u==en) return cur_flow;
    for(int v=0;v<=n+1;v++)//注意此处是所有顶点
    {
        if(cap[u][v]>0 && dis[u]+1==dis[v])
        {
            int flow=dfs(v,Min(dt,cap[u][v]));
            cap[u][v]-=flow;
            cap[v][u]+=flow;
            dt-=flow;
        }
    }
    return cur_flow-dt;
}
int dinic()
{
    int cur_flow,ans=0;
    while(bfs())//一次bfs可以找到几条增广路
    {
        while(cur_flow=dfs(st,inf))
            ans+=cur_flow;
    }
    return ans;
}
 
int main()
{
    int np,nc;
    while(scanf("%d%d%d%d",&n,&np,&nc,&m) != EOF)
    {
        int u,v,c;
        memset(cap,0,sizeof(cap));
        for(int i=0;i<m;i++)
        {
            scanf(" (%d,%d)%d",&u,&v,&c);
            cap[u+1][v+1]+=c;
        }
        for(int i=0;i<np;i++)
        { 
            scanf(" (%d)%d",&v,&c);
            cap[0][v+1]+=c;
        }
        for(int i=0;i<nc;i++)
        {
            scanf(" (%d)%d",&u,&c);
            cap[u+1][n+1]+=c;
        }
        st=0,en=n+1;
        printf("%d\n",dinic());
     }
     return 0;
 }