[BZOJ 1070][SCOI2007]修车：费用流

点击这里查看原题

费用流经典好题。将每个维修人员拆成n个点，第i个维修人员的第j个点与第k辆车连边表示k是i修的倒数第j辆车，因此费用为v[k][i]*j。
将源点与所有车连容量为1费用为0的边，车与维修人员拆成的所有点连容量为1费用为v[k][i]*j的边，维修人员拆成的所有点与汇点连容量为1费用为0的边
（这道题我分别WA、RE了一次，WA是没注意到读入的顺序是m,n而不是n,m；RE是因为对边的数量的估算过小

/*
User:Small
Language:C++
Problem No.:1070
*/
#include<bits/stdc++.h>
#define ll long long
#define inf 999999999
using namespace std;
const int M=805;
int n,m,v[65][15],fir[M],tot=1,ans,dis[M];
bool vis[M];
struct edge{
    int v,w,c,nex;
}e[200005];
void add(int u,int v,int w,int c){
    e[++tot]=(edge){v,w,c,fir[u]};
    fir[u]=tot;
    e[++tot]=(edge){u,0,-c,fir[v]};
    fir[v]=tot;
}
bool spfa(int s,int t){
    memset(vis,0,sizeof(vis));
    memset(dis,127,sizeof(dis));
    vis[t]=1;
    dis[t]=0;
    queue<int> q;
    q.push(t);
    while(!q.empty()){
        int u=q.front();
        q.pop();
        vis[u]=0;
        for(int i=fir[u];i;i=e[i].nex){
            int v=e[i].v;
            if(dis[u]+e[i^1].c<dis[v]&&e[i^1].w){
                dis[v]=dis[u]+e[i^1].c;
                if(!vis[v]){
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
    return dis[s]<inf;
}
int dfs(int u,int f,int t){
    vis[u]=1;
    if(u==t) return f;
    int sum=0;
    for(int i=fir[u];i;i=e[i].nex){
        int v=e[i].v;
        if(dis[u]-e[i].c==dis[v]&&e[i].w&&!vis[v]){
            int d=dfs(v,min(f,e[i].w),t);
            ans+=d*e[i].c;
            e[i].w-=d;
            e[i^1].w+=d;
            f-=d;
            sum+=d;
            if(!f) return sum;
        }
    }
    return sum;
}
void zkw(int s,int t){
    while(spfa(s,t)){
        vis[t]=1;
        while(vis[t]){
            memset(vis,0,sizeof(vis));
            dfs(s,inf,t);
        }
    }
}
int main(){
    freopen("data.in","r",stdin);//
    scanf("%d%d",&m,&n);
    for(int i=1;i<=n;i++){
        add(1,i+1,1,0);
        for(int j=1;j<=m;j++){
            scanf("%d",&v[i][j]);
            for(int k=1;k<=n;k++)
                add(i+1,(j-1)*n+k+n+1,1,v[i][j]*k);
        }
    }
    for(int i=1;i<=m;i++){
        for(int j=1;j<=n;j++){
            add((i-1)*n+j+n+1,m*n+n+2,1,0);
        }
    }
    zkw(1,m*n+n+2);
    printf("%.2lf",(double)ans/n);
    return 0;
}