[bzoj1070][SCOI2007]修车 费用流

1070: [SCOI2007]修车

Time Limit: 1 Sec  Memory Limit:
128 MB

[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

HINT

数据范围: (2<=M<=9,1<=N<=60), (1<=T<=1000)

Source

太晚了，不想写题解了

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#define INF 0x7fffffff
using namespace std;
const int N = 605;
int n,m,cnt=1,dis
,q
,last
,S,T,ans,t[65][10];
bool mark
;
struct Edge{
int to,next,v,c;
}e[100005];
void insert( int u, int v, int w, int c ){
e[++cnt].to = v; e[cnt].next = last[u]; last[u] = cnt; e[cnt].v = w; e[cnt].c = c;
e[++cnt].to = u; e[cnt].next = last[v]; last[v] = cnt; e[cnt].v = 0; e[cnt].c = -c;
}
bool spfa(){
memset(mark,0,sizeof(mark));
for( int i = S; i <= T; i++ ) dis[i] = INF;
int head = 0, tail = 1;
dis[T] = 0; mark[T] = 1; q[0] = T;
while( head != tail ){
int now = q[head++]; if( head == T ) head = 0;
for( int i = last[now]; i; i = e[i].next )
if( e[i^1].v && dis[e[i].to] > dis[now] - e[i].c ){
dis[e[i].to] = dis[now] - e[i].c;
if( !mark[e[i].to] ){
mark[e[i].to] = 1; q[tail++] = e[i].to; if( tail == T ) tail = 0;
}
}
mark[now] = 0;
}
return dis[S] != INF;
}
int dfs( int x, int f ){
mark[x] = 1;
int w, used = 0;
if( x == T ) return f;
for( int i = last[x]; i; i = e[i].next )
if( !mark[e[i].to] && e[i].v && dis[e[i].to] == dis[x] - e[i].c ){
w = dfs( e[i].to, min(f-used,e[i].v) );
ans += w*e[i].c; e[i].v -= w; e[i^1].v += w; used += w;
if( used == f ) return f;
}
return used;
}
void zkw(){
while( spfa() ){
mark[T] = 1;
while( mark[T] ){
memset(mark,0,sizeof(mark));
dfs(S,INF);
}
}
}
int main(){
scanf("%d%d", &n, &m);
S = 0; T = n*m+m+1;
for( int i = 1; i <= m; i++ )
for( int j = 1; j <= n; j++ )
scanf("%d", &t[i][j]);
for( int i = 1; i <= n*m; i++ ) insert( S, i, 1, 0 );
for( int i = n*m+1; i <= n*m+m; i++ ) insert( i, T, 1, 0 );
for( int i = 1; i <= n; i++ )
for( int j = 1; j <= m; j++ )
for( int k = 1; k <= m; k++ )
insert( (i-1)*m+j, n*m+k, 1, t[k][i]*j );
zkw();
printf("%.2lf", (double)ans/m);
return 0;
}