BZOJ1070修车

1070: [SCOI2007]修车
Time Limit: 1 Sec Memory Limit: 162 MB
Submit: 3771 Solved: 1532
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)
zkw费用流。。
建图比较奇特。。于是求教了hzwer学长。。
把每个工人拆成N个点。记为A[i,j]表示第i个工人修倒数第j辆车。
每个车跟所有N*M个工人拆出的点连边。流量为1，费用为time[i,j]*k。
源和每辆车连边，N*M个点和汇连边，流量都为1，费用同为0。
为什么这么构图呢？
考虑第i个工人，他修第j辆车只对后面要修的车有影响，而前面修过的车已经对当前没有影响了。
而这个影响就是后面每个将要修理的车都多等待了time的时间。
其他边流量都为1是显然的，每辆车修一次，每个工人一个时段只能修理一辆车。
然而读入优化的时候没加负数的读入WA了无数次，至今不知为何。。
感谢DaD3zZ&hzwer。。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m,ans,sum,cnt=1,head,tail,h[500001],q[500001],dis[500001],T,a[101][10];
bool mark[500001],vis[500001];
struct node
{
    int to,next,v,c;
};
node edge[500001];

int read()
{
    int w=0,c=1;
    char ch=getchar();
    while (ch<'0' || ch>'9')
      {
        if (ch=='-')
          c=-1;
        ch=getchar();
      }
    while (ch>='0' && ch<='9')
      {
        w=w*10+ch-'0';
        ch=getchar();
      }
    return w*c;
}

void add(int u,int v,int w,int cost)
{
    cnt++;
    edge[cnt].next=h[u];
    h[u]=cnt;
    edge[cnt].to=v;
    edge[cnt].v=w;
    edge[cnt].c=cost;
}

void insert(int u,int v,int w,int cost)
{
    add(u,v,w,cost);
    add(v,u,0,-cost);
}

bool spfa()
{
    memset(vis,false,sizeof(vis));
    for (int i=0;i<=T;i++)
      dis[i]=0x7fffffff;
    head=0;
    tail=1;
    q[0]=T;
    vis[T]=true;
    dis[T]=0;
    while (head<tail)
      {
        int now=q[head];
        head++;
        vis[now]=false;
        for (int i=h[now];i;i=edge[i].next)
          if (edge[i^1].v && dis[now]-edge[i].c<dis[edge[i].to])
            {
              dis[edge[i].to]=dis[now]-edge[i].c;
              if (!vis[edge[i].to])
                {
                  vis[edge[i].to]=true;
                  q[tail++]=edge[i].to;
                }
            }
      }
    return dis[0]!=0x7fffffff;
}

int dfs(int x,int f)
{
    int w,used=0;
    mark[x]=true;
    if (x==T)
      return f;
    for (int i=h[x];i;i=edge[i].next)
      if (dis[edge[i].to]==dis[x]-edge[i].c && edge[i].v && !mark[edge[i].to])
        {
            w=f-used;
            w=dfs(edge[i].to,min(w,edge[i].v));
            ans+=w*edge[i].c;
            edge[i].v-=w;
            edge[i^1].v+=w;
            used+=w;
            if (used==f)
              return f;
        }
    return used;
}

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++)
        a[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,a[k][i]*j);
    while (spfa())
      {
        mark[T]=true;
        while (mark[T])
          {
            memset(mark,false,sizeof(mark));
            sum+=dfs(0,0x7fffffff);
          }
      }
    printf("%.2lf",(double)ans/m);
    return 0;
}