hdu 2426 Interesting Housing Problem (KM算法)

Interesting Housing Problem

Time Limit: 10000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2388    Accepted Submission(s): 879

Problem Description

For any school, it is hard to find a feasible accommodation plan with every student assigned to a suitable apartment while keeping everyone happy, let alone an optimal one. Recently the president of University ABC, Peterson, is facing a similar problem. While Peterson does not like the idea of delegating the task directly to the class advisors as so many other schools are doing, he still wants to design a creative plan such that no student is assigned to a room he/she dislikes, and the overall quality of the plan should be maximized. Nevertheless, Peterson does not know how this task could be accomplished, so he asks you to solve this so-called "interesting" problem for him.
Suppose that there are N students and M rooms. Each student is asked to rate some rooms (not necessarily all M rooms) by stating how he/she likes the room. The rating can be represented as an integer, positive value meaning that the student consider the room to be of good quality, zero indicating neutral, or negative implying that the student does not like living in the room. Note that you can never assign a student to a room which he/she has not rated, as the absence of rating indicates that the student cannot live in the room for other reasons.
With limited information available, you've decided to simply find an assignment such that every student is assigned to a room he/she has rated, no two students are assigned to the same room, and the sum of rating is maximized while satisfying Peterson's requirement. The question is … what exactly is the answer?

 

Input

There are multiple test cases in the input file. Each test case begins with three integers, N, M, and E (1 <= N <= 500, 0 <= M <= 500, 0 <= E <= min(N * M, 50000)), followed by E lines, each line containing three numbers, S _i, R _i, V _i, (0 <= S _i < N, 0 <= R _i < M, |V _i| <= 10000), describing the rating V _i given by student S _ifor room R _i. It is guaranteed that each student will rate each room at most once.
Each case is followed by one blank line. Input ends with End-of-File.

 

Output

For each test case, please output one integer, the requested value, on a single line, or -1 if no solution could be found. Use the format as indicated in the sample output.

 

Sample Input

  
  

   
   3 5 5
0 1 5
0 2 7
1 1 6
1 2 3
2 4 5

1 1 1
0 0 0

1 1 0
  
  

 

Sample Output

  
  

   
   Case 1: 18
Case 2: 0
Case 3: -1
  
  

 

Source

2008 Asia Hangzhou Regional Contest Online

m!=n,把边权为负数的不读入就AC了

#include"stdio.h"
#include"string.h"
#define N 505
#define mmax(a,b) ((a)>(b)?(a):(b))
#define mmin(a,b) ((a)<(b)?(a):(b))
const int inf=(int)1e8;
int g[N][N],slack[N],n,m;
int link[N],lx[N],ly[N];
int visx[N],visy[N];
int find(int k)
{
    int i;
    visx[k]=1;
    for(i=0;i<m;i++)
    {
        if(visy[i])
            continue;
        int d=lx[k]+ly[i]-g[k][i];
        if(d==0)
        {
            visy[i]=1;
            if(link[i]==-1||find(link[i]))
            {
                link[i]=k;
                return 1;
            }
        }
        else
            slack[i]=mmin(slack[i],d);
    }
    return 0;
}
int KM()
{
    int i,j;
    memset(ly,0,sizeof(ly));
    memset(link,-1,sizeof(link));
    for(i=0;i<n;i++)
    {
        lx[i]=-inf;
        for(j=0;j<m;j++)
        {
            lx[i]=mmax(lx[i],g[i][j]);
        }
    }
    for(i=0;i<n;i++)
    {
        for(j=0;j<m;j++)
        {
            slack[j]=inf;
        }
        while(1)
        {
            memset(visx,0,sizeof(visx));
            memset(visy,0,sizeof(visy));
            if(find(i))
                break;
            else
            {
                int d=inf;
                for(j=0;j<m;j++)
                {
                    if(!visy[j])
                        d=mmin(slack[j],d);
                }
                for(j=0;j<n;j++)
                {
                    if(visx[j])
                        lx[j]-=d;
                }
                for(j=0;j<m;j++)
                {
                     if(visy[j])
                        ly[j]+=d;
                }
            }
        }
    }
    int ans=0,cnt=0;
    for(i=0;i<m;i++)
    {
        if(link[i]==-1)
            continue;
		if(g[link[i]][i]<0)
			break;
		cnt++;
        ans+=g[link[i]][i];
    }
    if(cnt<n)
        return -1;
    return ans;
}
int main()
{
    int i,j,u,v,w,e,cnt=1;
    while(scanf("%d%d%d",&n,&m,&e)!=-1)
    {
        int max=mmax(n,m);
        for(i=0;i<max;i++)
        {
            for(j=0;j<max;j++)
            {
                g[i][j]=-inf;
            }
        }
        while(e--)
        {
            scanf("%d%d%d",&u,&v,&w);
            if(w>=0)
				g[u][v]=w;
        }
		printf("Case %d: %d\n",cnt++,KM());
    }
    return 0;
}