【网络流#1】hdu 3549

因为坑了无数次队友 要开始学习网络流了，先从基础的开始，嗯~

这道题是最大流的模板题，用来测试模板好啦~

Edmonds_Karp模板 with 前向星

时间复杂度o(V*E^2)

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<queue>
#define eps 0.000001
#define MAXN 20
#define MAXM 2005
#define INF (1<<30)
using namespace std;
int i,j,k,n,m,x,y,T,head[MAXN],num,w,pre[MAXN],lin[MAXN];
struct edgenode
{
    int to,next,w;
} map[MAXM];

void add_edge(int id,int x,int y,int w)
{
    map[id].to=y;
    map[id].w=w;
    map[id].next=head[x];
    head[x]=id;
}

bool bfs(int s,int t)
{
    int u,v;
    queue <int> q;
    memset(pre,-1,sizeof(pre));
    pre[s]=s;
    q.push(s);
    while(!q.empty())
    {
        u=q.front();
        q.pop();
        for(int i=head[u];i!=-1;i=map[i].next)
        {
            v=map[i].to;
            if(pre[v]==-1&&map[i].w>0)
            {
                pre[v]=u;
                lin[v]=i;
                if (v==t) return true;
                q.push(v);
            }
        }
    }
    return false;
}

int Edmonds_Karp(int s,int t)
{
    int flow=0,d,i;
    while(bfs(s,t))
    {
        d=INF;
        for(i=t;i!=s;i=pre[i])
            d=min(d,map[lin[i]].w);
        for(i=t;i!=s;i=pre[i])
        {
            map[lin[i]].w-=d;
            map[lin[i]^1].w+=d;
        }
        flow+=d;
    }
    return flow;
}

void init()
{
    memset(head,-1,sizeof(head));
    scanf("%d%d",&n,&m);
    num=0;
    for (i=0;i<m;i++)
    {
        scanf("%d%d%d",&x,&y,&w);
        add_edge(i*2,x,y,w);
        add_edge(i*2+1,y,x,0);
    }
}

int main()
{
    scanf("%d",&T);
    for (int cas=1;cas<=T;cas++)
    {
        init();
        printf("Case %d: %d
",cas,Edmonds_Karp(1,n));
    }
}

 最大流的优化还有dinic和isap

dinic传送：dinic