Jamie's Contact Groups

最新推荐文章于 2020-04-21 10:02:42 发布

LeeBooL

最新推荐文章于 2020-04-21 10:02:42 发布

阅读量821

点赞数

分类专栏：图论网络流

本文链接：https://blog.csdn.net/u013570474/article/details/44685601

版权

图论同时被 2 个专栏收录

47 篇文章 0 订阅

订阅专栏

网络流

9 篇文章 0 订阅

订阅专栏

二分 + 多值匹配

思路：一对多的二分图的多重匹配。二分图的多重匹配算法的实现类似于匈牙利算法，对于集合C中的元素xi，找到一个与其相连的元素yi后，检查匈牙利算法的两个条件是否成立，若yi未被匹配，则将
xi，yi匹配。否则，如果与yi匹配的元素已经达到上限，那么在所有与yi匹配的元素中选择一个元素，检查是否能找到一条增广路径，如果能，则让出位置，让xi与yi匹配。
二分求出limit，知道找到可以构成多重匹配的最小限制limit，在main函数中二分搜索。

#include <iostream>
#include <string.h>
#include <stdlib.h>
#include <stdio.h>
#include <math.h>

using namespace std;
#define MAXN 1111
int cx[MAXN][MAXN];
bool gra[MAXN][MAXN], vis[MAXN];
int n, m, limit;

bool dfs( int u )
{
    for( int v = 0; v < m; v++)
    {
        if(gra[u][v] && !vis[v])
        {
            vis[v] = 1;
            if(cx[v][0] < limit)
            {
                cx[v][++cx[v][0]] = u;
                return true;
            }
            for( int i = 1; i <= cx[v][0]; i++)
              if(dfs(cx[v][i]))
              {
                  cx[v][i] = u;
                  return true;
              }

        }
    }
    return false;
}
bool check( )
{
    int res = 0;
    for( int i = 0; i < m; i++)
         cx[i][0] = 0;
    for( int i = 0; i < n; i++)
    {
        memset(vis, 0, sizeof(vis));
        if(!dfs(i)) return false;
    }
    return true;
}

void solve( )
{
    int l, r;
    l = 0, r = n;
    int ans = 0;
    while(l < r)
    {
        limit = ( l + r ) >> 1;
        if(check())
            r = limit;
        else
            l = limit + 1;
    }
    printf("%d\n",r);
}

int main()
{
    char str[2100];
    while(scanf("%d %d",&n, &m) != EOF &&( n + m))
    {
        memset(gra, 0, sizeof(gra));
        getchar();
        for( int i = 0; i < n; i++)
        {
            gets(str);
            int len  = strlen(str);
            for( int j = 0; j < len; j++)
            {
                if(str[j] >= '0' && str[j] <= '9')
                {
                    int v = 0;
                    while(str[j] >= '0' && str[j] <= '9')
                        {
                            v = v*10 + (str[j] - '0');
                            j++;
                        }
                     gra[i][v] = 1;
                }
            }
        }
        solve();
    }
    return 0;
}