BZOJ1070：[SCOI2007]修车——题解

http://www.lydsy.com/JudgeOnline/problem.php?id=1070

https://www.luogu.org/problemnew/show/P2053#sub

同一时刻有N位车主带着他们的爱车来到了汽车维修中心。维修中心共有M位技术人员，不同的技术人员对不同的车进行维修所用的时间是不同的。现在需要安排这M位技术人员所维修的车及顺序，使得顾客平均等待的时间最小。

说明：顾客的等待时间是指从他把车送至维修中心到维修完毕所用的时间。

请转至：http://www.cnblogs.com/luyouqi233/p/7953480.html （POJ3686）

所以省选考了道原题emmmm……

#include<cstdio>
#include<iostream>
#include<queue>
#include<cstring>
#include<algorithm>
#include<cctype>
using namespace std;
typedef long long ll;
const int INF=1e9;
const int N=;
const int M=;
const int P=N+M*N+;
inline int read(){
    int X=,w=;char ch=;
    while(!isdigit(ch)){w|=ch=='-';ch=getchar();}
    while(isdigit(ch))X=(X<<)+(X<<)+(ch^),ch=getchar();
    return w?-X:X;
}
struct node{
    int nxt,to,w,b;
}edge[(N+N*N*M+N*M)*];
int head[P],cnt=-;
void add(int u,int v,int w,int b){
    cnt++;
    edge[cnt].to=v;
    edge[cnt].w=w;
    edge[cnt].b=b;
    edge[cnt].nxt=head[u];
    head[u]=cnt;
    return;
}
int dis[P];
bool vis[P];
inline bool spfa(int s,int t,int n){
    deque<int>q;
    memset(vis,,sizeof(vis));
    for(int i=;i<=n;i++)dis[i]=INF;
    dis[t]=;q.push_back(t);vis[t]=;
    while(!q.empty()){
    int u=q.front();
    q.pop_front();vis[u]=;
    for(int i=head[u];i!=-;i=edge[i].nxt){
        int v=edge[i].to;
        int b=edge[i].b;
        if(edge[i^].w&&dis[v]>dis[u]-b){
        dis[v]=dis[u]-b;
        if(!vis[v]){
            vis[v]=;
            if(!q.empty()&&dis[v]<dis[q.front()]){
            q.push_front(v);
            }else{
            q.push_back(v);
            }
        }
        }
    }
    }
    return dis[s]<INF;
}
int ans=;
int dfs(int u,int flow,int m){
    if(u==m){
    vis[m]=;
    return flow;
    }
    int res=,delta;
    vis[u]=;
    for(int e=head[u];e!=-;e=edge[e].nxt){
        int v=edge[e].to;
    int b=edge[e].b;
        if(!vis[v]&&edge[e].w&&dis[u]-b==dis[v]){
            delta=dfs(v,min(edge[e].w,flow-res),m);
            if(delta){
                edge[e].w-=delta;
                edge[e^].w+=delta;
                res+=delta;
        ans+=delta*b;
                if(res==flow)break;
            }
        }
    }
    return res;
}
inline int costflow(int S,int T,int n){
    while(spfa(S,T,n)){
    do{
        memset(vis,,sizeof(vis));
        dfs(S,INF,T);
    }while(vis[T]);
    }
    return ans;
}
int main(){
    memset(head,-,sizeof(head));
    cnt=-;
    int m=read(),n=read();
    int S=n+m*n+,T=n+m*n+;
    for(int i=;i<=n;i++){
    add(S,i,,);
    add(i,S,,);
    for(int j=;j<=m;j++){
        int z=read();
        for(int k=;k<=n;k++){
        int p=j*n+k;
        add(i,p,,z*k);
        add(p,i,,-z*k);
        }
    }
    }
    for(int i=n+;i<=n+n*m;i++){
    add(i,T,,);
    add(T,i,,);
    }
    printf("%.2f\n",costflow(S,T,n+m*n+)*1.0/n);
    return ;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者：luyouqi233。 　　　　　　　　　　　　　　+

+欢迎访问我的博客：http://www.cnblogs.com/luyouqi233/+

+++++++++++++++++++++++++++++++++++++++++++