bzoj 1070: [SCOI2007]修车（费用流）

1070: [SCOI2007]修车

Time Limit: 1 Sec   Memory Limit: 128 MB
Submit: 5818   Solved: 2467
[ Submit][ Status][ Discuss]

Description

　　同一时刻有N位车主带着他们的爱车来到了汽车维修中心。维修中心共有M位技术人员，不同的技术人员对不同的车进行维修所用的时间是不同的。现在需要安排这M位技术人员所维修的车及顺序，使得顾客平均等待的时间最小。 说明：顾客的等待时间是指从他把车送至维修中心到维修完毕所用的时间。

Input

　　第一行有两个m,n，表示技术人员数与顾客数。 接下来n行，每行m个整数。第i+1行第j个数表示第j位技术人员维修第i辆车需要用的时间T。

Output

　　最小平均等待时间，答案精确到小数点后2位。

Sample Input

2 2

3 2

1 4

Sample Output

1.50

将每个工人都拆成n个点，源点为0，汇点为n+n*m+1

源点与所有的车(n个)连接一条费用为0，流量为1的边

所有的工人(n*m个)与汇点连接一条费用为0，流量为1的边

所有的车与对应的所有工人连接一条费用为对应修车时间，流量为1的边（这一行连接是不对的）

但是题目中工人修完一辆车，才能修下一辆车，这就是要将每个工人拆成n个点的原因

应该所有的车与对应的所有工人拆的n个点分别连接费用为对应修车时间*k的边

(1<=k<=n表示工人倒数第k个修这个车，因为这后面修的所有车都会慢一个单位的修车时间)

样例建图如下：

之后最小费用最大流就好了

#include<stdio.h>
#include<string.h>
#include<queue>
#include<algorithm>
using namespace std;
typedef struct
{
	int cost;
	int flow;
}Res;
Res road[707][707];
int t, a[66][12], link[707], vis[707], best[707];
int SPFA()
{
	int now, i;
	memset(best, 2, sizeof(best));
	memset(link, -1, sizeof(link));
	memset(vis, 0, sizeof(vis));
	queue<int> q;
	best[0] = 0, q.push(0);
	while(q.empty()==0)
	{
		now = q.front();
		q.pop();
		vis[now] = 0;
		for(i=1;i<=t;i++)
		{
			if(road[now][i].flow>0 && best[now]+road[now][i].cost<best[i])
			{
				best[i] = best[now]+road[now][i].cost;
				link[i] = now;
				if(vis[i]==0)
				{
					q.push(i);
					vis[i] = 1;
				}
			}
		}
	}
	if(best[t]<=1000000)
		return 1;
	return 0;
}
int main(void)
{
	int m, n, i, j, k, ans, x, pre;
	while(scanf("%d%d", &m, &n)!=EOF)
	{
		memset(road, 0, sizeof(road));
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=m;j++)
				scanf("%d", &a[i][j]);
		}
		t = n+n*m+1;
		for(i=1;i<=n;i++)
			road[0][i].flow = 1;
		for(i=n+1;i<=n+n*m;i++)
			road[i][t].flow = 1;
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=m;j++)
			{
				for(k=1;k<=n;k++)
				{
					road[i][j*n+k].flow = 1;
					road[i][j*n+k].cost = k*a[i][j];
					road[j*n+k][i].cost = -k*a[i][j];
				}
			}
		}
		ans = 0;
		while(SPFA())
		{
			x = t;
			while(x!=0)
			{
				pre = link[x];
				road[pre][x].flow -= 1;
				road[x][pre].flow += 1;
				ans += road[pre][x].cost;
				x = pre;
			}
		}
		printf("%.2f\n", ans*1.0/n);
	}
	return 0;
}