1070: [SCOI2007]修车

题目链接

题目大意：m个工人，n个顾客，并给出每个工人洗每辆车的时间。n个顾客同时到达，问平均最小等待时间（每个人从到达到洗完）是多少？

题解：看上去比较像贪心……我按照二分图最大权匹配建了图，然后愉快地过了样例，然后就WA了……
仔细分析了一下，发现这样搞没有累加上仍在等待的顾客的时间。如两个技术人员，三个顾客，那么必定有一个顾客会多等一会儿，这一段时间没有算进来

把每个工人分别拆成N个点，记为A[i,j]表示第i个工人修倒数第j辆车。由于工人修车只对未修理的顾客影响，可以处理每个点对全局的影响。连(i,A[j,k],1,k*T[i][j])，其中i为每辆车。费用则表示会让k个人每人等待T[i][j]的时间

然后按照刚才的二分图建法把连一下源和汇就好了。

我的收获：拆点劲啊

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
const int M=5050;
#define INF 0x3f3f3f3f
int n,m,t,st,ed,ans,tot,nod;
bool vis[M];
int head[M],d[M],pre[M],T[65][65],cnt[65][65];

queue<int> q;

struct edge{int fro,to,nex,c,val;}e[M*20];

void add(int i,int j,int ca,int co){e[t].fro=i,e[t].to=j,e[t].nex=head[i],e[t].c=ca,e[t].val=co;head[i]=t++;}
void insert(int i,int j,int w,int z){add(i,j,w,z),add(j,i,0,-z);}

bool spfa()
{
    for(int i=0;i<=ed;i++) vis[i]=0,d[i]=INF;
    d[st]=0;q.push(st);
    while(!q.empty())
    {
        int u=q.front();q.pop();vis[u]=false;
        for(int i=head[u];i!=-1;i=e[i].nex){
            int v=e[i].to;
            if(e[i].c&&d[v]>d[u]+e[i].val){
                d[v]=d[u]+e[i].val;
                pre[v]=i;
                if(!vis[v]){
                    vis[v]=true;
                    q.push(v);
                }
            }
        }
    }
    return d[ed]!=INF;
}
void flow()
{
    int mx=INF;
    for(int u=ed;u!=st;u=e[pre[u]].fro)
        mx=min(mx,e[pre[u]].c);
    for(int u=ed;u!=st;u=e[pre[u]].fro){
        e[pre[u]].c-=mx;e[pre[u]^1].c+=mx;
        ans+=mx*e[pre[u]].val;//费用 
    }
}

void work()
{
    while(spfa()) flow();
    printf("%.2lf",double(ans)/n);
}

void build()
{
    for(int i=1;i<=n;i++)
        insert(st,++nod,1,0);
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n;j++)
            cnt[i][j]=++nod,insert(nod,ed,1,0);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            for(int k=1;k<=n;k++)
                insert(i,cnt[j][k],1,k*T[i][j]);
}

void init()
{
    int x,y,z,w;t=0;
    memset(head,-1,sizeof(head));
    cin>>m>>n;st=0;ed=n*m+n+1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            scanf("%d",&T[i][j]);
    build();
}

int main()
{
    init();
    work();
    return 0;
}