最大流(Dinic临接表)

这个算法用的蛮多,不过我不太喜欢。。有点复杂

但还是要掌握的

没太看懂,当模板记下

#include <stdio.h>

#include <string.h>

#include <queue>

using namespace std;

#define MAXN 210

#define INF 0x3f3f3f3f


struct Edge

{

    int st, ed;

    int c;

    int next;

}edge[MAXN << 1];

int n, m;

int s, t;

int ans;

int e = 0;

int head[MAXN];//

int d[MAXN];


int min(int a, int b)

{

    return a < b ? a : b;

}


void init()

{

    int i;

    int a, b, c;

    s = 1;

    t = m;

    e = 0;

    ans = 0;

    memset(head, -1, sizeof(head));

    for(i = 1; i <= n; i++)

    {

        scanf("%d%d%d", &a, &b, &c);

        edge[e].st = a;

        edge[e].ed = b;

        edge[e].c = c;

        edge[e].next = head[a];

        head[a]= e++;

        edge[e].st = b;

        edge[e].ed = a;

        edge[e].next = head[b];

        head[b] = e++;

    }

}


int bfs()

{

    memset(d, -1, sizeof(d));

    queue<int> q;

    d[s] = 0;

    q.push(s);

    int i;

    int cur;

    while(!q.empty())

    {

        cur = q.front();

        q.pop();

        for(i = head[cur]; i != -1; i = edge[i].next)

        {

            if(d[edge[i].ed] == -1 && edge[i].c > 0)

            {

                d[edge[i].ed] = d[cur] + 1;

                q.push(edge[i].ed);

            }

        }

    }

    if(d[t] < 0)

        return 0;

    return 1;

}


int dinic(int x, int flow)

{

    if(x == t)

        return flow;

    int i, a;

    for(i = head[x]; i != -1; i = edge[i].next)

    {

        if(d[edge[i].ed] == d[x] + 1 && edge[i].c > 0 &&

           (a = dinic(edge[i].ed, min(flow, edge[i].c))))

        {

            edge[i].c -= a;

            edge[i^1].c += a;

            return a;

        }

    }

    return 0;

}


void solve()

{

    while(scanf("%d%d", &n, &m) != EOF)

    {

        init();

        while(bfs())

        {

            int increment;

            increment = dinic(1, INF);

            ans +=  increment;

        }

        printf("%d\n", ans);

    }

}


int main()

{

    //#ifdef LOCAL

    //  freopen("poj1273.txt", "r", stdin);

    // freopen(".txt", "w", stdout);

    //#endif

    solve();

    return 0;

}

/*

 6 5

 1 2 20

 1 3 40

 2 4 10

 2 5 10

 3 5 20

 4 5 10

 */


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值