网络流24题-1.飞行员配对方案问题

简单建图,建超源超汇
最后打印图
链接

#include <bits/stdc++.h>

#define INF 0x3f3f3f3f
using namespace std;
const int maxn = 550;
const int maxm = maxn * 50;
int n, m;
int to[maxm], nxt[maxm], w[maxm];
int fst[maxn], cnt, cur[maxn];
int s, t, d[maxn];

void init() {
    memset(fst, -1, sizeof(fst));
    cnt = 0;
}

void addedge(int u, int v, int val) {
    to[cnt] = v;
    w[cnt] = val;
    nxt[cnt] = fst[u];
    fst[u] = cnt++;

    to[cnt] = u;
    w[cnt] = 0;
    nxt[cnt] = fst[v];
    fst[v] = cnt++;
}

int bfs() {
    queue<int> Q;
    memset(d, 0, sizeof(d));
    d[s] = 1;
    Q.push(s);
    while (!Q.empty()) {
        int u = Q.front();
        Q.pop();
        for (int i = fst[u]; ~i; i = nxt[i]) {//因为初始是-1
            int v = to[i];
            if (w[i] > 0 && d[v] == 0) {
                d[v] = d[u] + 1;
                if (v == t) return 1;
                Q.push(v);
            }

        }
    }
    return d[t] != 0;
}

int dfs(int u, int flow) {
    if (u == t)
        return flow;
    int ans = 0, x = 0;
    for (int i = cur[u]; ~i; i = nxt[i]) {
        int v = to[i];
        if (w[i]>0 && d[v] == d[u]+1){
            x = dfs(v,min(flow-ans,w[i]));
            w[i]-=x;
            w[i^1]+=x;
            if (w[i]) cur[u] = i;
            ans+=x;
            if (ans == flow) return flow;
        }
    }
    if (ans == 0)d[u] =0;
    return ans;
}

int dinic(int n){
    int ans = 0;

    while (bfs()){
        for (int i=0;i<n;i++)
            cur[i]  = fst[i];
        ans+=dfs(s,INF);
    }
    return ans;
}
int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);

    cin >> m >> n;
    init();
    s = 0;t= n+1;
    for(int i =1;i<=m;i++){
        addedge(s,i,1);
    }
    for (int i=m+1;i<=n;i++)
        addedge(i,t,1);
    while (1){
        int x,y;
        cin >> x >> y;
        if (x == -1 && y == -1) break;
        addedge(x,y,1);
    }
    int ans = dinic(t);
    cout << ans << endl;
    if (ans == 0){
        cout << "No Solution!\n";
        return 0;
    }
    for (int i=m+1;i<=n;i++)
        for (int j = fst[i];~j;j = nxt[j]){
            if (to[j]!=t && w[j] == 1)  cout << to[j] <<' '<< i<<endl;
        }
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值