poj1087解题报告(最大流 EK算法)

题目大意:题意:有n个不同的插座,有m台不同的机器需要m种插头,有k组转换:插头A能由插头B转换而来。问这些机器最少有几台不能插上插座。

解题思路:使用网络流,以设备的插头类型和插座的插头类型组成一个图,建立超级源点,超级汇点,源点到设备的流量为1,插座到汇点的流量为1,设备和流量之间流量为无限大(可以任取),最后用设备数减去最大流量。(利用FLOYD算法进行闭包传递,根据转化器的关系)

#include<iostream>
#include<map>
#include<cstdio>
#include<string>
#include<vector>
#include<cstring>
#include<queue>
using namespace std;
const int maxn=400+10;
vector<string> names;
int ID(const string &s)
{
        for(int i=0;i<names.size();i++)
                if(s==names[i])  return i;
        names.push_back(s);
        return names.size() - 1;
}


const int INF = 1000000000;

struct Edge {
  int from, to, cap, flow;
  Edge(int u, int v, int c, int f):from(u),to(v),cap(c),flow(f) {}
};

struct EdmondsKarp {
  int n, m;
  vector<Edge> edges;    // 边数的两倍
  vector<int> G[maxn];   // 邻接表,G[i][j]表示结点i的第j条边在e数组中的序号
  int a[maxn];           // 当起点到i的可改进量
  int p[maxn];           // 最短路树上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));
    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;
  }
};
EdmondsKarp g;
int n, m, k;       // 插座个数,设备个数,转换器个数
int d[maxn][maxn]; // d[i][j]=1表示插头类型i可以转化为插头类型j
int target[maxn];  // 各个插座的类型
int device[maxn];  // 各个设备的类型
int main()
{
        while(scanf("%d",&n)!=EOF)
        {
                names.clear();
                string s1, s2;
                for(int i=0;i<n;i++)
                {
                        cin>>s1;
                        target[i]=ID(s1);
                }
                scanf("%d",&m);
                for(int i=0;i<m;i++)
                {
                        cin>>s1>>s2;
                        device[i]=ID(s2);
                }
                memset(d, 0, sizeof(d));
                 scanf("%d",&k);
                for(int i=0;i<k;i++)
                {
                        cin>>s1>>s2;
                        d[ID(s1)][ID(s2)]=1;
                }
                 int V = names.size();
                 for(int k=0;k<V;k++)
                        for(int i=0;i<V;i++)
                               for(int j=0;j<V;j++)
                     d[i][j] |= d[i][k] && d[k][j];
                 g.init(V+2);
                     //构造网络
                 for(int i = 0; i < m; i++)
      g.AddEdge(V, device[i], 1); // 源点->设备
    for(int i = 0; i < n; i++)
      g.AddEdge(target[i], V+1, 1); // 插座->汇点
      for(int i = 0; i < m; i++)
      for(int j = 0; j < n; j++)
        if(d[device[i]][target[j]]) g.AddEdge(device[i], target[j], INF); // 设备->插座
    int r = g.Maxflow(V, V+1);
    cout << m-r << "\n";

        }
        return 0;
}


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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值