【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)

题解
对于第i个人，拆成n个点，第k个点表示修倒数第k辆车，每个拆点向车连流量为1，费用为k*time[i][j]

代码

#include<bits/stdc++.h>
#define inf 0x7fffffff
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m,tot=1,ans,t[70][10],T;
int dis[1005],from[1005],head[1005];
bool inq[1005];
struct edge{int from,to,next,c,w;}e[100001];
inline void ins(int u,int v,int w,int c)
{
    tot++;
    e[tot].to=v;e[tot].from=u;
    e[tot].w=w;e[tot].c=c;
    e[tot].next=head[u];head[u]=tot;
}
inline void insert(int u,int v,int w,int c){ins(u,v,w,c);ins(v,u,0,-c);}
bool spfa()
{
    for (int i=1;i<=T;i++) dis[i]=inf;queue<int>q;
    dis[0]=0;q.push(0);inq[0]=1;
    while (!q.empty())
    {
        int now=q.front();q.pop();
        for (int i=head[now];i;i=e[i].next)
            if (dis[e[i].to]>dis[now]+e[i].c&&e[i].w)
            {
                dis[e[i].to]=dis[now]+e[i].c;from[e[i].to]=i;
                if(!inq[e[i].to]){inq[e[i].to]=1;q.push(e[i].to);} 
            }
        inq[now]=0;
    }
    if (dis[T]==inf) return 0;return 1;
}
void mcf()
{
    int x=inf,i=from[T];
    while (i){x=min(x,e[i].w);i=from[e[i].from];}
    i=from[T];
    while (i){ans+=x*e[i].c;e[i].w-=x;e[i^1].w+=x;i=from[e[i].from];}
}
int main()
{
    n=read();m=read();T=n*m+m+1;
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n;j++)
            t[i][j]=read();
    for(int i=1;i<=n*m;i++)
        insert(0,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);
    while (spfa()) mcf();
    printf("%.2lf",(double)ans/m);
    return 0;
}