题解 [SCOI2007]修车

题面

解析

这题要拆点。。

首先，证明一个式子：

设修理员M修了N辆车，

且修每辆车的时间为W1，W2....W_N。

那么，这个修理员一共花的时间就为：W1*N+W2*(N-1)+...+W_N*1。

因此，若i号修理员修第j辆车的时间为c[i][j]，

将c[i][j]拆成1...n个点，

将其中第k个点与j号顾客（j号车）相连，

费用为c[i][j]*k，

最后将修理员与源点相连，顾客与汇点相连（反过来也可以），

流量为1，费用为0。

跑费用流就行了。

上AC代码：

//#include<bits/stdc++.h>
#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;

inline int read(){
    int sum=,f=;char ch=getchar();
    while(ch>'' || ch<''){if(ch=='-')f=-;ch=getchar();}
    while(ch>='' && ch<=''){sum=sum*+ch-'';ch=getchar();}
    return f*sum;
}

const int INF=0x3f3f3f3f;
struct node{
    int to,next,v,w;
}e[];
struct hh{
    int fa,edge;
}pre[];
int n,m,s,t;
int c[][]/*第i辆车给第j个人修的时间*/;
int head[],cnt=;
int d[],inq[],mi[];

void add(int x,int y,int v,int w){
    //printf("%d -> %d : %d\n",x,y,w);
     e[++cnt].to=head[x];
    e[cnt].next=y;
    e[cnt].v=v;e[cnt].w=w;
    head[x]=cnt;
    e[++cnt].to=head[y];
    e[cnt].next=x;
    e[cnt].v=;e[cnt].w=-w;
    head[y]=cnt;
}

bool spfa(){
    memset(d,0x3f,sizeof(d));
    memset(inq,,sizeof(inq));
    memset(mi,0x3f,sizeof(mi));
    queue <int> que;
    que.push(s);
    d[s]=;
    while(!que.empty()){
        int x=que.front();
        que.pop();
        inq[x]=;
        for(int i=head[x];i;i=e[i].to){
            int k=e[i].next;
            if(!e[i].v||d[k]<=d[x]+e[i].w) continue;
            d[k]=d[x]+e[i].w;
            pre[k].fa=x;pre[k].edge=i;
            mi[k]=min(mi[x],e[i].v);
            if(!inq[k]) que.push(k);
            inq[k]=;
        }
    }
    return d[t]!=INF;
}

void EK(){
    int ans=;
    while(spfa()){
        for(int i=t;i!=s;i=pre[i].fa){
            e[pre[i].edge].v-=mi[t];
            e[pre[i].edge^].v+=mi[t];
        }
        ans+=d[t]*mi[t];
    }
    printf("%.2lf\n",(double)((double)ans/(double)n));
}

int main(){
//    freopen("fixed.in","r",stdin);
//    freopen("fixed.out","w",stdout);
    m=read();n=read();
    s=+(m+)*n;t=+(m+)*n;
    for(int i=;i<=n;i++){
        for(int j=;j<=m;j++){
            c[i][j]=read();
        }
    }
    for(int i=;i<=m;i++){
        for(int j=;j<=n;j++){
            add(s,i+m*j,,);
            for(int k=;k<=n;k++){
                add(i+m*j,k+m*n+n,,c[k][i]*j);
            }
        }
    }
    for(int i=;i<=n;i++){
        add(i+m*n+n,t,,);
    }
    EK();
    return ;
}