网络流

真的是浅浅的了解了一下网络流。
其中增广路径和残余网络,需要理解清楚,这对于算法的理解是有帮助的。
可以看看别人上传的资料,资料(我也想传,不知道是不是侵权)。
也推荐一个图解比较详细的博客,图解
接下去就看模板题和Dinic模板吧。


poj3281

/*
 *Memory 0.8
 *Time 79
 *本题居然和网络流有关,我的天啊,原谅我的见识短浅
 *全程套用《挑战程序设计竞赛第二版》模板
 *此题用到了Dinic算法和拆点技巧
 *边的容量都是 1 ,如果牛能按照题意选择正确的食物和饮料
 *那么从S到T的流量增加 1
*/
#include <iostream>
#include <vector>
#include <queue>
#include <algorithm>
using namespace std;
const int MAXN = 408;
const int MAXF = 102;
const int MAXD = 102;
const int INF = 1<<30;
struct Edge{
    int to, cap, rev;
};//下一节点,容量,反向边
vector <Edge> G[MAXN];
int level[MAXN], iter[MAXN];
bool likeF[MAXN][MAXF], likeD[MAXN][MAXD];
void addEdge(int u, int v, int cap) {
    G[u].push_back((Edge) {v, cap, G[v].size()});
    G[v].push_back((Edge) {u, 0, G[u].size() - 1});
}
void bfs(int s, int n) {
    for(int i = 0; i <= n; ++i) level[i] = -1;
    queue<int> Q;
    level[s] = 0;
    Q.push(s);
    while(!Q.empty()) {
        int v = Q.front(); Q.pop();
        for(int i = 0; i < G[v].size(); ++i) {
            Edge &e = G[v][i];
            if(e.cap > 0 && level[e.to] < 0) {
                level[e.to] = level[v] + 1;
                Q.push(e.to);
            }
        }
    }
}
int dfs(int v, int t, int f) {
    if(v == t) return f;
    for(int &i = iter[v]; i < G[v].size(); ++i) {
        Edge &e = G[v][i];
        if(e.cap > 0 && level[v] < level[e.to]) {
            int d = dfs(e.to, t, min(f, e.cap));
            if(d > 0) {
                e.cap -= d;
                G[e.to][e.rev].cap += d;
                return d;
            }
        }
    }
    return 0;
}
int maxFlow(int s, int t, int n) {
    int flow = 0;
    while(true) {
        bfs(s, n);
        if(level[t] < 0) return flow;
        for(int i = 0; i <= n; ++i) iter[i] = 0;
        int f;
        while((f = dfs(s, t, INF)) > 0) flow += f;
    }
}
void solve(int n, int f, int d) {
    //0 - n-1:表示食物一侧的牛
    //n - 2n-1:表示饮料一侧的牛
    //2n - 2n+f-1:表示食物
    //2n+f - 2n+f+d-1:表示饮料
    int s = 2 * n + f + d, t = s + 1;
    for(int i = 0; i < f; ++i) addEdge(s, 2 * n + i, 1);
    for(int i = 0; i < d; ++i) addEdge(2 * n + f + i, t, 1);
    for(int i = 0; i < n; ++i) {
        addEdge(i, n + i, 1);
        for(int j = 0; j < f; ++j)
            if(likeF[i][j]) addEdge(2 * n + j, i, 1);
        for(int j = 0; j < d; ++j)
            if(likeD[i][j]) addEdge(n + i, 2 * n + f + j, 1);
    }
    cout<<maxFlow(s, t, t)<<endl;
}
int main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    int N, F, D;
    int f, d, fi, di;
    while(cin>>N>>F>>D) {
        for(int i = 0; i <= N; ++i)
            for(int j = 0; j <= N; ++j)
                likeD[i][j] = likeF[i][j] = false;
        for(int i = 0; i <= N; ++i) G[i].clear();

        for(int i = 0; i < N; ++i) {
            cin>>f>>d;
            for(int j = 0; j < f; ++j) {
                cin>>fi;
                likeF[i][fi - 1] = true;
            }
            for(int j = 0; j < d; ++j) {
                cin>>di;
                likeD[i][di - 1] = true;
            }
        }
        solve(N, F, D);
    }
    return 0;
}

Dinic算法模板:

struct Edge{
    int to, cap, rev;
};
vector <Edge> G[MAXN];
int level[MAXN], iter[MAXN];

void addEdge(int u, int v, int cap) {
    G[u].push_back((Edge) {v, cap, G[v].size()});
    G[v].push_back((Edge) {u, 0, G[u].size() - 1});
}
void bfs(int s, int n) {
    for(int i = 0; i <= n; ++i) level[i] = -1;
    queue<int> Q;
    level[s] = 0;
    Q.push(s);
    while(!Q.empty()) {
        int v = Q.front(); Q.pop();
        for(int i = 0; i < (int)G[v].size(); ++i) {
            Edge &e = G[v][i];
            if(e.cap > 0 && level[e.to] < 0) {
                level[e.to] = level[v] + 1;
                Q.push(e.to);
            }
        }
    }
}
int dfs(int v, int t, int f) {
    if(v == t) return f;
    for(int &i = iter[v]; i < G[v].size(); ++i) {
        Edge &e = G[v][i];
        if(e.cap > 0 && level[v] < level[e.to]) {
            int d = dfs(e.to, t, min(f, e.cap));
            if(d > 0) {
                e.cap -= d;
                G[e.to][e.rev].cap += d;
                return d;
            }
        }
    }
    return 0;
}
int maxFlow(int s, int t, int n) {
    int flow = 0;
    while(true) {
        bfs(s, n);
        if(level[t] < 0) return flow;
        for(int i = 0; i <= n; ++i) iter[i] = 0;
        int f;
        while((f = dfs(s, t, INF)) > 0) flow += f;
    }
}

如果以后有时间在怼它(~.~)

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值