本题n个工作分给n个人,求最优,最差方案,所以是KM算法。当然也可以用费用流建模。
建模方法如下:
建立源点ST和汇点ED,由ST向n个工作连费用为0,流量为1的边,由每个人向ED连费用为0,流量为1的边,由n个工作向n个人连流量为无穷大,费用为wij的边,最后跑一遍最小费用最大流和最大费用最大流即可。
代码如下:
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <vector>
#include <queue>
#define M 1001011
#define INF 0x7fffffff/3
#define ST 2001
#define ED 2002
using namespace std;
int n,m;
int num=0;
int ans=0;
int p[M];
int fmap[201][201];
int cost[M],flag[M],pre[M];
int q[4*M];
struct node
{
int a;
int b;
int w;
int f;
int next;
}str1[M];
void init()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
int x;
scanf("%d",&x);
fmap[i][j]=x;
}
}
}
void add(int x,int y,int z,int fl)
{
str1[num].a=x;
str1[num].b=y;
str1[num].w=z;
str1[num].f=fl;
str1[num].next=p[x];
p[x]=num++;
swap(x,y);
str1[num].a=x;
str1[num].b=y;
str1[num].w=-z;
str1[num].f=0;
str1[num].next=p[x];
p[x]=num++;
}
bool spfa(int src)
{
memset(flag,0,sizeof(flag));
for(int i=0;i<=ED;i++)
cost[i]=INF;
cost[src]=0;
pre[src]=-1;
flag[src]=1;
int l=0,r=0;
q[r++]=src;
while(l<r)
{
int k=q[l++];
flag[k]=0;
for(int i=p[k];i!=-1;i=str1[i].next)
{
int v=str1[i].b;
if(cost[v]>cost[k]+str1[i].w&&str1[i].f>0)
{
cost[v]=cost[k]+str1[i].w;
pre[v]=i;
if(!flag[v])
{
flag[v]=1;
q[r++]=v;
}
}
}
}
if(cost[ED]==INF)
return false;
return true;
}
void start()
{
int minx=INF;
for(int i=pre[ED];i!=-1;i=pre[str1[i^1].b])
minx=min(minx,str1[i].f);
for(int i=pre[ED];i!=-1;i=pre[str1[i^1].b])
{
str1[i].f-=minx;
str1[i^1].f+=minx;
}
ans+=cost[ED]*minx;
}
void spfaflow()
{
ans=0;
while(spfa(ST))
start();
}
void start1()
{
num=0;
memset(p,-1,sizeof(p));
for(int i=1;i<=n;i++)
add(ST,i,0,1);
for(int i=1;i<=n;i++)
add(i+n,ED,0,1);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
add(i,j+n,fmap[i][j],INF);
}
}
spfaflow();
printf("%d\n",ans);
}
void start2()
{
num=0;
memset(p,-1,sizeof(p));
for(int i=1;i<=n;i++)
add(ST,i,0,1);
for(int i=1;i<=n;i++)
add(i+n,ED,0,1);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
add(i,j+n,-fmap[i][j],INF);
}
}
spfaflow();
printf("%d\n",-ans);
}
int main()
{
freopen("T18.in","r",stdin);
freopen("T18.out","w",stdout);
init();
start1();
start2();
return 0;
}