bzoj1070 [SCOI2007]修车

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

HINT

数据范围: ((2le Mle 9,1le Nle 60), (1le Tle 1000))

Solution

啊♂费用流

把 (n) 和 (m) 调换一下，符合习惯。

• 把每个工人拆成 (n) 个点。记为 (a[i,j]) 表示第 (i) 个工人修倒数第 (j) 辆车。

• 每个车跟所有 (n imes m) 个工人拆出的点连边。流量为 (1) ，费用为 (time[i,j] imes k) 。

• 源和每辆车连边， (n imes m) 个点和汇连边，流量都为 (1) ，费用为 (0) 。

#include<bits/stdc++.h>
using namespace std;

#define rep(i, a, b) for (int i = a; i <= b; i++)
#define N 1010
#define INF 0x7fffffff

inline int read() {
	int x = 0, flag = 1; char ch = getchar(); while (!isdigit(ch)) { if (!(ch ^ '-')) flag = -1; ch = getchar(); }
	while (isdigit(ch)) x = (x << 1) + (x << 3) + ch - '0', ch = getchar(); return x * flag;
}

int n, m, t[61][10];
struct edge { int u, v, next, w, c; }eg[100001]; int tot = 1;
int head[N], from[N];
int ans = 0;

queue<int> q; bool inq[N];
int dis[N];
bool spfa() {
	rep(i, 0, 1001) dis[i] = INF; dis[0] = 0; q.push(0), inq[0] = 1;
	while(!q.empty()) {
		int u = q.front(), v; q.pop(); //printf("%d
", dis[u]);
		for(int i = head[u]; i; i = eg[i].next) if(eg[i].w && dis[(v = eg[i].v)] > dis[u] + eg[i].c) {
			dis[v] = dis[u] + eg[i].c, from[v] = i; if(!inq[v]) q.push(v), inq[v] = 1;
		}
		inq[u] = 0;
	}
	return dis[1001] != INF;
}

inline void update() {
	int tmp = INF;
	for(int i = from[1001]; i; i = from[eg[i].u]) tmp = min(tmp, eg[i].w);
    for(int i = from[1001]; i; i = from[eg[i].u]) eg[i].w -= tmp, eg[i ^ 1].w += tmp, ans += eg[i].c * tmp;
}

inline void ins(int u, int v, int w, int c) { eg[++tot] = edge { u, v, head[u], w, c }; head[u] = tot; }
inline void add(int u, int v, int w, int c) { ins(u, v, w, c); ins(v, u, 0, -c); }
int main() {
//	freopen("c.in", "r", stdin), freopen("c.out", "w", stdout);
	n = read(), m = read(); rep(i, 1, m) rep(j, 1, n) t[i][j] = read();
	rep(i, 1, n * m) add(0, i, 1, 0); rep(i, n * m + 1, n * m + m) add(i, 1001, 1, 0);
	rep(i, 1, n) rep(j, 1, m) rep(k, 1, m) add(m * (i - 1) + j, k + n * m, 1, t[k][i] * j);
//	rep(i, 1, tot) printf("%d %d %d %d %d
", eg[i].u, eg[i].v, eg[i].next, eg[i].c, eg[i].w);
	while(spfa()) update(); printf("%.2lf", 1.0 * ans / m); return 0;
}