POJ-1087(最大流+EK算法)

A Plug for UNIX

POJ-1087

• 这一题也是最大流的应用，但是题目有点绕。
• 题目的意思就是有插座和插头，只不过这里使用设备取代了插头。只有插座和设备进行匹配了。
• 题目要注意的是那个适配器的问题，适配器给出的是 A,B 的模式，但是这里并不是指A可以取代B而是指B可以取代A，所以建图的时候要注意这个细节。
• 还有一个问题就是题目的规模是未知的，所以maxn尽量设置的较大，这样才保险。

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<string>
#include<vector>
#include<queue>
#include<map>
using namespace std;
const int maxn=1003;
const int INF=0X3F3F3F3F;
int n,m;
int s,t;
int m1;
int k;
map<string,int> ma;//设备
struct Edge {
    int from, to, cap, flow;
    Edge(int u, int v, int c, int f) : from(u), to(v), cap(c), flow(f) {}
};
vector<Edge> edges;//这里的边是实际边数的两倍，包括反向边
vector<int> G[maxn];//邻接表，G[i][j]表示结点i的第j条边在edges数组中的序号
int a[maxn];//a[i]表示起点到i的可改进量
int p[maxn];//edges数组中的编号，最短路图上p的入弧编号
void init(int n) {
    for (int i = 0; i <= n; i++) 
        G[i].clear();
    edges.clear();
}
void AddEdge(int from, int to, int cap) {
    edges.push_back(Edge(from, to, cap, 0));
    edges.push_back(Edge(to, from, 0, 0));//反向弧，容量为0
    m = edges.size();
    G[from].push_back(m - 2);
    G[to].push_back(m - 1);
}
int Maxflow(int s, int t) {
    int flow = 0;
    for (;;) {
        memset(a, 0, sizeof(a));
        queue<int> Q;
        Q.push(s);
        a[s] = INF;
        while (!Q.empty()) {
            int x = Q.front();
            Q.pop();
            for (int i = 0; i < G[x].size(); i++) {
                Edge& e = edges[G[x][i]];
                if (!a[e.to] && e.cap > e.flow) {
                    p[e.to] = G[x][i];
                    a[e.to] = min(a[x], e.cap - e.flow);
                    Q.push(e.to);
                }
            }
            if (a[t]) break;
        }
        if (!a[t]) break;
        for (int u = t; u != s; u = edges[p[u]].from) {
            edges[p[u]].flow += a[t];
            edges[p[u] ^ 1].flow -= a[t];
        }
        flow += a[t];
    }
    return flow;
}
int main(){
    cin>>n;
    s=0;
    t=maxn-1;
    int numn=0;
    for(int i=0;i<n;i++){
        string a;
        cin>>a;
        ma[a]=++numn;
        AddEdge(s,numn,1);
    }
    cin>>m1;
    for(int i=0;i<m1;i++){
        string a,b;
        cin>>a>>b;
        ma[a]=++numn;
        if(!ma.count(b))
            ma[b]=++numn;
        AddEdge(ma[b],ma[a],1);
        AddEdge(ma[a],t,1);
    }
    cin>>k;
    for(int i=0;i<k;i++){
        string a,b;
        cin>>a>>b;
        if(!ma.count(a))
            ma[a]=++numn;
        if(!ma.count(b))
            ma[b]=++numn;
        AddEdge(ma[b],ma[a],INF);
    }
    int ans=Maxflow(s,t);
    cout<<m1-ans<<endl;
    //system("pause");
    return 0;
}