【bzoj 1070】修车（费用流）

传送门biu~
把每个工人拆成$n$个点，源点连每个$n\times m$的点，流量为$1$，费用为$0$。$n\times m$的点连每辆车，流量为$1$，费用为$time_{i,j}\times k$。每辆车连汇点，流量为$1$，费用为$0$。因为每个工人修车所带来的影响是它和它后面的车的等待时间都增加了$time_{i,j}$。

#include<bits/stdc++.h>
using namespace std;
const int INF=1e9;
int n,m,S,T,cost,t[65][15];
int head[1005],nex[100005],to[100005],cap[100005],val[100005],tp=1;
int fir[1005],dis[1005];
bool b[1005];
inline void add(int x,int y,int c,int v){
    nex[++tp]=head[x];
    head[x]=tp;
    to[tp]=y;
    cap[tp]=c;
    val[tp]=v;
}
inline void Insert(int x,int y,int c,int v){add(x,y,c,v);add(y,x,0,-v);}
inline bool spfa(){
    queue<int>q;
    for(int i=S;i<=T;++i)       dis[i]=INF,b[i]=false;
    dis[S]=0;q.push(S);
    while(!q.empty()){
        int x=q.front();q.pop();b[x]=false;
        for(int i=head[x];i;i=nex[i]){
            if(cap[i] && dis[x]+val[i]<dis[to[i]]){
                dis[to[i]]=dis[x]+val[i];
                if(!b[to[i]]){
                    q.push(to[i]);
                    b[to[i]]=true;
                }
            }
        }
    }
    return dis[T]^INF;
}
int dfs(int x,int now){
    if(x==T || !now){
        cost+=now*dis[T];
        return now;
    }   
    int c=0;
    b[x]=true;
    for(int &i=fir[x];i;i=nex[i]){
        if(!b[to[i]] && cap[i] && dis[to[i]]==dis[x]+val[i]){
            int f=dfs(to[i],min(now,cap[i]));
            now-=f;
            cap[i]-=f;
            cap[i^1]+=f;
            c+=f;
            if(!now)    break;
        }
    }
    return c;
}
inline int Dinic(){
    int c=0;
    while(spfa()){
        for(int i=S;i<=T;++i)       fir[i]=head[i];
        c+=dfs(S,INF);
    }
    return c;
}
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;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);
    Dinic();
    printf("%.2lf",(double)cost/m);
    return 0;           
}