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

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

正解：网络流最小费用最大流

解题报告：建图有点神奇，对于每一个维修人员建ｎ个点，表示这是他修的倒数第几辆，那么当前的车只会影响后面的车的等待时间，假如现在是倒数第ｊ辆车，那么花费就是ｊ＊ｆ【ｉ】【ｋ】，ｆ【ｉ】【ｋ】表示第ｉ个维修工人修序号为ｋ的车的花费。然后ｎ辆车连上ｎ＊ｍ个点，容量为１，车再向源点点连容量为１的边，ｎ＊ｍ个点向汇点连容量为１的边，跑ｓｐｆａ。

#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#define MIN(a,b) a<b?a:b
#define RG register
const int N = 5000;
const int M = 1000000;
const int inf = 2147483641;

using namespace std;

int gi(){
    char ch=getchar();int x=0;
    while(ch<'0' || ch>'9') ch=getchar();
    while(ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar();
    return x;
}

int n,m,q,ans,cnt=1,S,T;
int fa[N],f[N][N],nxt[M],c[N][N],head[N],s[M],dis[M],vis[M];

struct date{
    int from,to,val,s;
}nn[M];

void link(int l,int r,int val,int s){
    nn[++cnt]=(date){l,r,val,s},nxt[cnt]=head[l],head[l]=cnt;
    nn[++cnt]=(date){r,l,0,-s},nxt[cnt]=head[r],head[r]=cnt;
    return;
}

int spfa(){
    int l=0,r=1;
    for (RG int i=S; i<=T; ++i) dis[i]=inf;
    s[1]=0,dis[0]=0,vis[0]=1;
    while(l<r){
        ++l;
        for (RG int i=head[s[l]]; i; i=nxt[i])
            if (nn[i].val && dis[nn[i].to]>dis[s[l]]+nn[i].s){
                fa[nn[i].to]=i;
                dis[nn[i].to]=dis[s[l]]+nn[i].s;
                if (!vis[nn[i].to]){
                    vis[nn[i].to]=1;s[++r]=nn[i].to;
                }
            }
        vis[s[l]]=0;
    }
    if (dis[T]==inf) return 0;
    int tmp=inf;
    for (RG int i=fa[T]; i; i=fa[nn[i].from]) tmp=MIN(tmp,nn[i].val);
    for (RG int i=fa[T]; i; i=fa[nn[i].from]) nn[i].val-=tmp,nn[i^1].val+=tmp;
    ans+=tmp*dis[T];
    return 1;
}

int main(){
    m=gi(),n=gi();
    for (RG int i=1; i<=n; ++i) link(S,i,1,0);
    for (RG int i=1; i<=n; ++i)
        for (RG int j=1; j<=m; ++j)
            f[j][i]=gi();
    ans=n,q=n;
    for (RG int i=1; i<=m; ++i)
        for (RG int j=1; j<=q; ++j)
            c[i][j]=++ans;
    T=ans+1;
    ans=0;
    for (RG int i=1; i<=n; ++i)
        for (RG int j=1; j<=m; ++j)
            for (RG int k=1; k<=q; ++k)
                link(i,c[j][k],1,f[j][i]*k);
    for (RG int j=1; j<=m; ++j)
        for (RG int k=1; k<=q; ++k)
            link(c[j][k],T,1,0);
    ans=0;
    while(spfa());
    printf("%.2lf",(ans*1.0)/n);
    return 0;
}