BZOJ1070 SCOI2007 修车 费用流

题意：给定N辆车和M个人，每个人修每辆车的时间不同，安排修理顺序使得修理这N辆车时间总和最小

题解：

将每个工作人员拆点，某个车b向ai连边表示工作人员a倒数第i个修理该车，显然其对答案的贡献（该边费用）就是w[b][a]*i（其后所有等待的车都要受其影响，而之前的车都不会）。

构图两排点，左边是车，右边是拆点后的工作人员，建立S和T跑最小费用流即可。

printf坏掉了死活保留不了两位小数，度娘了一种输出方法……

#include <queue>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <climits>
#include <iomanip>
#include <iostream>
#include <algorithm>
using namespace std;

const int MAXN=3000+2;
const int MAXM=80000+2;
struct HASH{
    int u;
    HASH *next;
    HASH(){}
    HASH(int _u,HASH *_next):u(_u),next(_next){}
}*table[MAXN],mem[MAXM];
struct EDGE{
    int u,v,w,c;
    EDGE(){}
    EDGE(int _u,int _v,int _c,int _w):u(_u),v(_v),c(_c),w(_w){}
}e[MAXM];
int N,M,cnt=2,ans,cur[MAXN],d[MAXN];
bool flag[MAXN];
queue<int> q;

void Insert(int u,int v,int c,int w){
    table[u]=&(mem[cnt]=HASH(cnt,table[u])),e[cnt++]=EDGE(u,v,c,w);
    table[v]=&(mem[cnt]=HASH(cnt,table[v])),e[cnt++]=EDGE(v,u,0,-w);
}

bool SPFA(int s,int t){
    for(int i=s;i<=t;i++) d[i]=INT_MAX;
    d[s]=0,flag[s]=1,q.push(s);

    int x;
    while(!q.empty()){
        x=q.front(),q.pop();
        for(HASH *p=table[x];p;p=p->next)
            if(e[p->u].c && d[e[p->u].v]>d[x]+e[p->u].w){
                d[e[p->u].v]=d[x]+e[p->u].w,cur[e[p->u].v]=p->u;
                if(!flag[e[p->u].v]) flag[e[p->u].v]=1,q.push(e[p->u].v);
            }
        flag[x]=0;
    }
    return d[t]<INT_MAX;
}

int Find(int s,int t){
    int ret=0,c=INT_MAX;
    for(int i=cur[t];i;i=cur[e[i].u]) c=min(c,e[i].c);
    for(int i=cur[t];i;i=cur[e[i].u]){
        e[i].c-=c,e[i^1].c+=c;
        ret+=e[i].w*c;
    }
    return ret;
}

int main(){
    scanf("%d %d",&M,&N);
    for(int i=1;i<=N;i++) Insert(0,i,1,0);
    for(int i=1,w;i<=N;i++)
        for(int j=1;j<=M;j++){
            scanf("%d",&w);
            for(int k=1;k<=N;k++) Insert(i,N+(k-1)*M+j,1,k*w);
        }
    for(int i=1;i<=M;i++)
        for(int j=1;j<=N;j++)
            Insert(N+(j-1)*M+i,N+N*M+1,1,0);

    while(SPFA(0,N+N*M+1)) ans+=Find(0,N+N*M+1);
    cout << setiosflags(ios::fixed);
    cout.precision(2);
    cout << (double)ans/N << endl;

    return 0;
}