KM poj 2195

题意：给出一个地图，地图上有人和房子，问如何分配哪个人去哪个房子，走的路最短？

这道题是个完备匹配的情况下，问怎么才能走的路最少，可以用KM来做。

只不过KM算法是用来求解最大最优值，所以我们得改一下数据，将每个人去房子的路程都改为负数。

最后再得出 -KM（）即刻

 1 #include<cstdio>
 2 #include<algorithm>
 3 #include<string.h>
 4 #include<math.h>
 5 #include<iostream>
 6 using namespace std;
 7 const int maxn=1e2+10;
 8 const int inf=0x3f3f3f3f;
 9 int lx[maxn],ly[maxn];
10 int match[maxn];
11 int visx[maxn],visy[maxn];
12 struct node
13 {
14     int x,y;
15 }a[maxn],b[maxn];
16 int G[maxn][maxn];
17 int numa,numb;
18 int dfs(int k)
19 {
20     visx[k]=1;
21     for(int i=1;i<=numa;i++){
22         if(!visy[i]&&G[k][i]==lx[k]+ly[i]){
23             visy[i]=1;
24             if(!match[i]||dfs(match[i])){
25                 match[i]=k;
26                 return 1;
27             }
28         }
29     }
30     return 0;
31 }
32 int KM()
33 {
34     for(int i=1;i<=numa;i++){
35         lx[i]=-inf,ly[i]=0;
36         for(int j=1;j<=numa;j++)
37             lx[i]=max(lx[i],G[i][j]);
38     }
39     for(int k=1;k<=numa;k++){
40         while(1){
41             memset(visx,0,sizeof(visx));
42             memset(visy,0,sizeof(visy));
43             if(dfs(k)) break;
44             int mn=inf;
45 
46             for(int i=1;i<=numa;i++) if(visx[i])
47             for(int j=1;j<=numa;j++) if(!visy[j])
48             mn=min(mn,lx[i]+ly[j]-G[i][j]);
49             if(mn==inf) return -1;
50             for(int i=1;i<=numa;i++) if(visx[i]) lx[i]-=mn;
51             for(int i=1;i<=numa;i++) if(visy[i]) ly[i]+=mn;
52         }
53     }
54     int ans=0;
55     for(int i=1;i<=numa;i++)
56     if(match[i])
57     ans+=G[match[i]][i];
58     return ans;
59 }
60 void init()
61 {
62     numa=numb=0;
63     memset(match,0,sizeof(match));
64     memset(G,0,sizeof(G));
65 }
66 int main()
67 {
68     int n,m;
69     while(scanf("%d%d",&n,&m)!=EOF){
70         if(n==0&&m==0) break;
71         init();
72         char t;
73         for(int i=0;i<n;i++){
74             for(int j=0;j<m;j++){
75                 cin>>t;
76                 if(t=='H'){
77                     numa++;
78                     a[numa].x=i;
79                     a[numa].y=j;
80                 }
81                 else if(t=='m'){
82                     numb++;
83                     b[numb].x=i;
84                     b[numb].y=j;
85                 }
86             }
87         }
88         for(int i=1;i<=numa;i++){
89             int x1=a[i].x,y1=a[i].y;
90             for(int j=1;j<=numb;j++){
91                 int x2=b[j].x,y2=b[j].y;
92                 int w=fabs(x1-x2)+fabs(y1-y2);
93                 G[i][j]=-w;
94             }
95         }
96         printf("%d\n",-KM());
97     }
98     return 0;
99 }