bzoj1070: [SCOI2007]修车

传送门
网络流。
把每个工人拆成N个点。记为A[i,j]表示第i个工人修倒数第j辆车。
每个车跟所有N*M个工人拆出的点连边。流量为1，费用为time[i,j]*k。
源和每辆车连边，N*M个点和汇连边，流量都为1，费用同为0。
考虑第i个工人，他修第j辆车只对后面要修的车有影响，而前面修过的车已经对当前没有影响了。
而这个影响就是后面每个将要修理的车都多等待了time的时间。
其他边流量都为1是显然的，每辆车修一次，每个工人一个时段只能修理一辆车。
所以构图是正确的。
只要跑一遍网络流就行了。

#include<cstdio>
#include<cmath>
#include<ctime>
#include<cstring>
#include<iostream>
#include<algorithm>
#define T 1001
#define mo 300005
using namespace std;
struct orz{
    int to,u,v,ne,fr;
}b[100005];
int d[10005],q[10005],fr[10005],c[10005],v[10005];
int n,m,tot=1,ans,a[61][10];
void ins(int u,int v,int w,int x){
    tot++;
    b[tot].fr=u;
    b[tot].to=v;
    b[tot].u=w;
    b[tot].v=x;
    b[tot].ne=c[u];
    c[u]=tot;
}
void inss(int u,int v,int w,int c){
    ins(u,v,w,c); ins(v,u,0,-c);
}
int bfs(){
    int h=0,t=1,x;
    memset(d,100,sizeof(d));
    memset(v,0,sizeof(v));
    memset(fr,0,sizeof(fr));
    v[0]=1;
    d[0]=q[1]=0;
    while (h!=t){
        h=h%mo+1;
        x=q[h];
        v[x]=0;
        for (int i=c[x];i;i=b[i].ne)
            if (b[i].u&&d[b[i].to]>d[x]+b[i].v){
                d[b[i].to]=d[x]+b[i].v;
                fr[b[i].to]=i;
                if (!v[b[i].to]){
                    t=t%mo+1;
                    q[t]=b[i].to;
                    v[b[i].to]=1;
                } 
            }
    }
    if (d[T]>1000000000) return 0;
    return 1;
}
void del(){
    int x=1000000000;
    for (int i=fr[T];i;i=fr[b[i].fr]) x=min(x,b[i].u);
    for (int i=fr[T];i;i=fr[b[i].fr]){
        b[i].u-=x;
        b[i^1].u+=x;
        ans+=b[i].v*x;
    }
} 
int main(){
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++)
        for (int j=1;j<=n;j++) scanf("%d",&a[i][j]);
     for (int i=1;i<=n*m;i++) inss(0,i,1,0);
     for (int i=n*m+1;i<=n*m+m;i++) inss(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++)
                inss((i-1)*m+j,n*m+k,1,a[k][i]*j);
    while (bfs()) del();
    printf("%.2lf",(double)ans/m);
    return 0;
}