hdu 2255 奔小康赚大钱二分权值最大匹配

最新推荐文章于 2024-10-18 11:16:53 发布

永远热血沸腾

最新推荐文章于 2024-10-18 11:16:53 发布

阅读量272

点赞数

CC 4.0 BY-SA版权

分类专栏：图论 hdu

本文链接：https://blog.youkuaiyun.com/xinag578/article/details/49981713

hdu 同时被 2 个专栏收录

273 篇文章

订阅专栏

图论

146 篇文章

订阅专栏

本文介绍了一种使用最小费用最大流算法解决实际问题的方法，并通过实例展示了算法的应用过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<stdio.h>
#include<string.h>
#include<queue>
#include<vector>
using namespace std;
const int N=1024;
const int inf=100000000;

int n,m,w[N][N],aug,mat[N],lx[N],ly[N];
int visx[N],visy[N],slack[N];

int dfs(int u)
{
    visx[u]=1;
    for(int v=1;v<=m;v++)
    {
        int t=lx[u]+ly[v]-w[u][v];
        if(t) slack[v]=min(slack[v],t);
        else
        {
            if(!visy[v])
            {
                visy[v]=1;
                if(mat[v]==-1||dfs(mat[v]))
                {
                    mat[v]=u;
                    return 1;
                }
            }
        }
    }
    return 0;
}

int km()
{
    int i,j,ret=0;
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++)
        lx[i]=max(lx[i],w[i][j]);
    for(int x=1;x<=n;x++)
    {
        for(i=0;i<=m;i++) slack[i]=inf;
        while(1)
        {
            memset(visx,0,sizeof(visx));
            memset(visy,0,sizeof(visy));
            if(dfs(x)) break;
            aug=inf;
            for(i=1;i<=m;i++)
                if(!visy[i]) aug=min(aug,slack[i]);
            for(i=1;i<=n;i++)
                if(visx[i]) lx[i]-=aug;
            for(i=1;i<=m;i++)
                if(visy[i]) ly[i]+=aug;
                else slack[i]-=aug;
        }
    }
    int cnt=0;
    for(i=1;i<=m;i++)
    {
        if(mat[i]==-1||w[mat[i]][i]==-inf) continue;
        cnt++;
        ret+=w[mat[i]][i];
    }
    if(cnt!=n) return -1;
    return ret;
}

void init()
{
    memset(mat,-1,sizeof(mat));
    memset(lx,-0x7f,sizeof(lx));
    memset(ly,0,sizeof(ly));
    memset(w,-0x7f,sizeof(w));
}

int main()
{
    int i,j;
    while(~scanf("%d",&n))
    {
        init();
        m=n;
        for(i=1;i<=n;i++)
        {
            for(j=1;j<=n;j++)
            {
                scanf("%d",&w[i][j]);
            }
        }
        printf("%d\n",km());
    }
}