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)

题解

神奇的费用流，将m个工人复制成n份，第i个工人的第j个分身表示这辆车是i工人修的倒数第j辆车（网上的题解说的到底是这个工人修的倒数第j辆车还是所有人修的倒数第j辆车，让人非常懵逼），从这个点向汇点连一条容量为1的边，因为工人修的倒数第j辆车只能有一辆（非常得显然），费用为j乘以这个工人修理这辆车所需的时间，因为除了这辆车的车主以外，还有j-1个人也在等这辆车修完（这辆车是这个工人修的倒数第j辆车嘛）。从每辆车向这n*m个分身都连一条边，然后跑费用流就可以了。
然而这为什么是对的呢？（我为了搞明白至少看了20个题解）
按照平时的思路，所有的从车流出的流量都会最终流向修这辆车最快的那个人，而这肯定是错误的。但是在这道题的做法中不会出现这种现象，因为这个工人修理的倒数第一辆车只能有一辆，其他的车所花费的时间都已经被翻了k倍，而费用流的性质保证取出最优解。

ps：注意点的个数和连边！点的个数不是（n+1）*m，是（m+1）*n！

CODE：

#include<queue>
#include<cstdio>
#include<cstring>
using namespace std;
const int INF=1e9;
struct edge
{
    int next,to,remain,cost;
}a[200000];
int head[100000];
int dis[100000];
int u[100000];
int f[100000];
bool b[100000];
int n,m,x,S,T,num=1,ans;
queue<int>q;
inline int min(int a,int b){return a<b?a:b;}
inline void add(int x,int y,int z,int w)
{
    a[++num].next=head[x],a[num].to=y,a[num].remain=z,a[num].cost=w,head[x]=num;
    a[++num].next=head[y],a[num].to=x,a[num].cost=-w,head[y]=num;
}
inline bool SPFA()
{
    memset(dis,0x3f,sizeof(dis));
    dis[S]=0;q.push(S);b[S]=1;u[S]=INF;
    while(!q.empty())
    {
        int tmp=q.front();q.pop();
        b[tmp]=0;
        for(int i=head[tmp];i;i=a[i].next)
          if(a[i].remain&&dis[a[i].to]>dis[tmp]+a[i].cost)
          {
            dis[a[i].to]=dis[tmp]+a[i].cost;
            f[a[i].to]=i;
            u[a[i].to]=min(u[tmp],a[i].remain);
            if(!b[a[i].to]) b[a[i].to]=1,q.push(a[i].to);
          }
    }
    if(dis[T]>INF) return 0;
    ans+=dis[T]*u[T];
    int p=T;
    while(p!=S)
    {
        a[f[p]].remain-=u[T];
        a[f[p]^1].remain+=u[T];
        p=a[f[p]^1].to;
    }
    return 1;
}
inline void mcmf(){while(SPFA());}
int main()
{
    scanf("%d%d",&m,&n);
    S=(m+1)*n+1,T=S+1;
    for(int i=1;i<=n;i++)
      add(S,i,1,0);
    for(int i=1;i<=n;i++)
      for(int j=1;j<=m;j++)
      {
        scanf("%d",&x);
        for(int k=1;k<=n;k++)
          add(i,n*j+k,1,k*x);
      }
    for(int i=n+1;i<=(m+1)*n;i++)
      add(i,T,1,0);
    mcmf();
    printf("%.2lf",ans*1.0/n);
    return 0;
}