srupのメモ帳

競プロで解いた問題や勉強したことを記録していくメモ帳

aoj GRL_7 A Matching - Bipartite Matching

問題

問題概要

省略

解法

2部マッチングの確認問題

ミス

蟻本のやつを写経した。

コード

#include <iostream>
#include <cstdio>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
#define rep(i,n) for(int i=0;i<(n);i++)
const int INF = 1e9;
const int MAX_V = 300;

int X, Y;
struct Flow{
    struct edge{
        int to, cap, rev;
    };
    vector<edge> G[MAX_V];//隣接リスト
    bool used[MAX_V];

    void add_edge(int from, int to, int cap){
        G[from].push_back((edge){to, cap, (int)G[to].size()});//from -> to
        G[to].push_back((edge){from, 0, (int)G[from].size() - 1});//to -> from
    }

    //増加パスを探す
    int dfs(int v, int t, int f){
        if(v == t) return f;
        used[v] = true;
        for (int i = 0; i < G[v].size(); ++i){
            edge &e = G[v][i];
            if(!used[e.to] && e.cap > 0){
                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;
    }
    //sからtへの最大流
    int max_flow(int s, int t){
        int flow = 0;
        while(1){
            memset(used, 0, sizeof(used));
            int f = dfs(s, t, INF);
            if(f == 0) return flow;
            flow += f;
        }
    }
};

int main(void){
    int e; cin >> X >> Y >> e;
    int s = X + Y;
    int t = X + Y + 1;
    Flow mf;
    //s -> x
    rep(i, X){
        mf.add_edge(s, i, 1);
    }
    // y -> t
    rep(i, Y){
        mf.add_edge(X + i, t, 1);
    }
    // x -> y
    rep(i, e){
        int tx, ty; cin >> tx >> ty;
        mf.add_edge(tx, X + ty, 1);
    }

    printf("%d\n", mf.max_flow(s, t));
}