bzoj1070: [SCOI2007]修车

网络流——最小费用流。

好久没写了板子都快忘了。将m个工人拆成n个点，点（i，j）表示让i工人在倒数第j辆车时去修连向这个点的车，那么显而易见代价为j*ti（后面的人要等）。

#include<bits/stdc++.h>
using namespace std;
#define N 70
#define M 15
#define INF 1e9
int n,m,cnt=2,ti[M][N],head[5005],p[5005],a[5005],ans,d[5005],S,T;
bool vis[5005];
queue<int>q;
struct edges{
    int fr,to,cap,flow,cost,next;
}e[100005];
void inser(int u,int v,int c,int co){
    e[cnt]=(edges){u,v,c,0,co,head[u]};head[u]=cnt++;
    e[cnt]=(edges){v,u,0,0,-co,head[v]};head[v]=cnt++;
}
bool spfa(){
    memset(d,0x3f,sizeof(d));
    q.push(S); d[S]=0; a[S]=INF;
    while(!q.empty()){
        int x=q.front(); q.pop(); vis[x]=0;
        for(int i=head[x];i;i=e[i].next){
            if(d[e[i].to]>d[x]+e[i].cost && e[i].cap>e[i].flow){
                d[e[i].to]=d[x]+e[i].cost; p[e[i].to]=i; a[e[i].to]=min(a[x],e[i].cap-e[i].flow);
                if(!vis[e[i].to]) vis[e[i].to]=1,q.push(e[i].to);
            }
        }
    }
    return d[T]<INF;
}
void minc(){
    ans+=a[T]*d[T];
    int u=T;
    while(u!=S){
        e[p[u]].flow+=a[T];
        e[p[u]^1].flow-=a[T];
        u=e[p[u]].fr;
    }
}
int main(){
    scanf("%d%d",&m,&n);
    S=0,T=n*m*2;
    for(int i=1;i<=n;i++)
       for(int j=1;j<=m;j++)
          scanf("%d",&ti[j][i]);
    for(int i=1;i<=m;i++)
       for(int j=1;j<=n;j++)
          inser((i-1)*n+j,T,1,0);
    for(int i=1;i<=n;i++) inser(S,n*m+i,1,0);
    for(int i=1;i<=n;i++)
       for(int j=1;j<=m;j++)
          for(int k=1;k<=n;k++)
          inser(n*m+i,(j-1)*n+k,1,k*ti[j][i]);
    while(spfa()) minc();
    printf("%.2lf
",(double)ans/n);
    return 0;
}