奔小康赚大钱
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 8258 Accepted Submission(s): 3665
Problem Description
传说在遥远的地方有一个非常富裕的村落,有一天,村长决定进行制度改革:重新分配房子。
这可是一件大事,关系到人民的住房问题啊。村里共有n间房间,刚好有n家老百姓,考虑到每家都要有房住(如果有老百姓没房子住的话,容易引起不安定因素),每家必须分配到一间房子且只能得到一间房子。
另一方面,村长和另外的村领导希望得到最大的效益,这样村里的机构才会有钱.由于老百姓都比较富裕,他们都能对每一间房子在他们的经济范围内出一定的价格,比如有3间房子,一家老百姓可以对第一间出10万,对第2间出2万,对第3间出20万.(当然是在他们的经济范围内).现在这个问题就是村领导怎样分配房子才能使收入最大.(村民即使有钱购买一间房子但不一定能买到,要看村领导分配的).
这可是一件大事,关系到人民的住房问题啊。村里共有n间房间,刚好有n家老百姓,考虑到每家都要有房住(如果有老百姓没房子住的话,容易引起不安定因素),每家必须分配到一间房子且只能得到一间房子。
另一方面,村长和另外的村领导希望得到最大的效益,这样村里的机构才会有钱.由于老百姓都比较富裕,他们都能对每一间房子在他们的经济范围内出一定的价格,比如有3间房子,一家老百姓可以对第一间出10万,对第2间出2万,对第3间出20万.(当然是在他们的经济范围内).现在这个问题就是村领导怎样分配房子才能使收入最大.(村民即使有钱购买一间房子但不一定能买到,要看村领导分配的).
Input
输入数据包含多组测试用例,每组数据的第一行输入n,表示房子的数量(也是老百姓家的数量),接下来有n行,每行n个数表示第i个村名对第j间房出的价格(n<=300)。
Output
请对每组数据输出最大的收入值,每组的输出占一行。
Sample Input
2 100 10 15 23
Sample Output
123
Source
题目大意:
n个村名选n个房子,每个村名都有自己喜欢的房子并且有个好感值,问怎么安排才能使好感值和最大
题目思路:
带权二分图,解决该类问题的算法主要是KM算法,也可以用费用流但时间复杂度更高
关于KM算法我推荐两个写的好的博客:
AC代码:
#include<cstring>
#include<cstdio>
#define min(x,y) (x<y?x:y)
#define max(x,y) (x>y?x:y)
#define abs(x) ((x)<0?-(x):(x))
const int maxn = 1e2+20;
const int inf = 0x3f3f3f3f;
class KM{
public:
int n,m,numh,numm;
int ph[maxn][2],pm[maxn][2];
int mp[maxn][maxn],slack[maxn],link[maxn],ex_l[maxn],ex_r[maxn];
bool vis_l[maxn],vis_r[maxn];
bool dfs(int u){
vis_l[u]=true;
for(int i=1;i<=n;i++){
if(!vis_r[i]){
if(ex_l[u]+ex_r[i]==mp[u][i]){
vis_r[i]=true;
if(link[i]==-1||dfs(link[i])){
link[i]=u;
return true;
}
}else {
slack[i]=min(slack[i],ex_l[u]+ex_r[i]-mp[u][i]);
}
}
}
return false;
}
int sove(){
for(int i=1;i<=n;i++){
memset(slack,0x3f,sizeof(slack));
while(1){
memset(vis_l,false,sizeof(vis_l));
memset(vis_r,false,sizeof(vis_r));
if(dfs(i))break;
int d = inf;
for(int j=1;j<=n;j++){
if(!vis_r[j])d=min(d,slack[j]);
}
for(int j=1;j<=n;j++){
if(vis_l[j])ex_l[j]-=d;
if(vis_r[j])ex_r[j]+=d;
else slack[j]-=d;
}
}
}
int res = 0;
for(int i=1;i<=n;i++){
if(link[i]!=-1){
res+=mp[link[i]][i];
}
}
return res;
}
void star(){
while(~scanf("%d%d",&n,&m),n+m){
init();
char str[maxn];
for(int i=1;i<=n;i++){
scanf("%s",str+1);
for(int j=1;j<=m;j++){
if(str[j]=='H')ph[++numh][0]=i,ph[numh][1]=j;
if(str[j]=='m')pm[++numm][0]=i,pm[numm][1]=j;
}
}
n=numh;
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
int x = abs(pm[i][0]-ph[j][0])+abs(pm[i][1]-ph[j][1]);
mp[i][j]=-x;
ex_l[i]=max(ex_l[i],mp[i][j]);
}
}
printf("%d\n",-sove());
}
}
void init(){
memset(link,-1,sizeof(link));
memset(ex_r,0,sizeof(ex_r));
memset(ex_l,-0x3f,sizeof(ex_l));
numh=numm=0;
}
}km;
int main()
{
km.star();
return 0;
}