poj2112--Optimal Milking（二分多重匹配）

题目来源：http://poj.org/problem?id=2112

题意

k个挤奶机和c头牛分布在牧场上，各自的距离以（k+c）*（k+c）的二维矩阵呈现。每个机器招待牛儿上限是m头。。。为了使走的最远的牛行走的路径最短。。。。写程序进行分配。。输出最大最小值。。。

思路

依旧是二分多重匹配。。。
需要floyed思想去初始化机器与牛的距离。
然后就是模板。，。。。
floyed我写错了。。，。。改了好几次。。。。都没发现。。。TΛT。。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=300;
const int INF=0x3f3f3f;
int mp[maxn][maxn],group[maxn][maxn],dist[maxn][maxn];
bool vis[maxn];
int k,c,m,top_limit;

void init()
{
    scanf("%d%d%d",&k,&c,&m);
    int n=k+c;
    for(int i=1; i<=n; i++)
    {
        for(int j=1; j<=n; j++)
        {
            scanf("%d",&dist[i][j]);
            if(!dist[i][j])
                dist[i][j]=INF;
        }
    }
    for(int t=1; t<=n; t++)
    {
        for(int i=1; i<=n; i++)
        {
            for(int j=1; j<=n; j++)
            {
                if(dist[i][j]>dist[i][t]+dist[t][j])
                    dist[i][j]=dist[i][t]+dist[t][j];
            }
        }
    }
}

bool dfs(int i)
{
    for(int j=1; j<=k; j++)
    {
        if(mp[i][j]&&!vis[j])
        {
            vis[j]=1;
            if(group[j][0]<m)
            {
                group[j][++group[j][0]]=i;
                return 1;
            }
            for(int p=1; p<=group[j][0]; p++)
            {
                if(dfs(group[j][p]))
                {
                    group[j][p]=i;
                    return 1;
                }
            }
        }
    }
    return 0;
}

bool judge()
{
    memset(group,0,sizeof(group));
    for(int i=1; i<=c; i++)
    {
        memset(vis,0,sizeof(vis));
        if(!dfs(i))
            return 0;
    }
    return 1;
}
void solve()
{
    int left=1,right=INF,res;
    while(left<=right)
    {
        top_limit=(left+right)>>1;
        memset(mp,0,sizeof(mp));
        for(int i=k+1; i<=k+c; i++)
            for(int j=1; j<=k; j++)
                if(dist[i][j]<=top_limit)
                    mp[i-k][j]=1;
        if(judge())
        {
            res=top_limit;
            right=top_limit-1;
        }
        else
        {
            left=top_limit+1;
        }
    }
    printf("%d\n",res);
}

int main()
{
    init();
    solve();
    return 0;
}