BZOJ1070: [SCOI2007]修车(费用流的思维建模/EK+动态边)

题意：传送门

题解：参考传送门

构造一个图，左边n∗m个点，表示第j个技术人员修的倒数第i辆车，右边n个点，表示第k辆车，源点向左边的点(w=1,cost=0)的边，因为每个人同一时刻只能修一辆车，左边m∗n个点再分别向右边n个点连边，左边第(i,j)个点向右边第k个点(w=1,cost=i∗t[i][j])的边，相当于此时有i个人在等待着这辆车修完，代价就是等待的人数∗修这辆车的时间,右边的点向汇点连(w=1,cost=0)的边，表示每辆车只能被一个人修，容易看出，此时最大流一定是n，每一条w=1的流表示第j个人修的倒数第i辆车是第k辆车，那么最小费用除以n就是答案。

附上代码：

#include<bits/stdc++.h>
using namespace std;
const int maxn=60+5;
const int maxm=10+5;
int m,n,T[maxm][maxn];
 
const int MAX_V=1000;
const int INF=0x3f3f3f3f;
 
typedef pair<int,int>P;
 
struct edge{
    int to,cap,cost,rev;
    edge(int _to,int _cap,int _cost,int _rev):to(_to),cap(_cap),cost(_cost),rev(_rev){}
};
 
int u,v,w;
vector<edge>G[MAX_V];
int h[MAX_V];
int dist[MAX_V];
int prevv[MAX_V],preve[MAX_V];
 
void add_edge(int from,int to,int cap,int cost)
{
    G[from].push_back(edge(to,cap,cost,G[to].size()));
    G[to].push_back(edge(from,0,-cost,G[from].size()-1));
}
 int res,f;
void min_cost_flow(int s,int t)
{
    while(233){
        priority_queue<P,vector<P>,greater<P> >que;
        fill(dist,dist+n*m+n+5,INF);
        dist[s]=0;
        que.push(P(0,s));
        while(!que.empty()){
            P p=que.top();que.pop();
            int v=p.second;
            if(dist[v]<p.first){
                continue;
            }
            for(int i=0;i<G[v].size();i++){
                edge &e=G[v][i];
                if(e.cap>0&&dist[e.to]>dist[v]+e.cost+h[v]-h[e.to]){
                    dist[e.to]=dist[v]+e.cost+h[v]-h[e.to];
                    prevv[e.to]=v;
                    preve[e.to]=i;
                    que.push(P(dist[e.to],e.to));
                }
            }
        }
        if(dist[t]==INF){
            return ;
        }
        for(int v=0;v<n*m+n+5;v++){
            h[v]+=dist[v];
        }
        int d=INF;
        for(int v=t;v!=s;v=prevv[v]){
            d=min(d,G[prevv[v]][preve[v]].cap);
        }
        f+=d;
        res+=d*h[t];
        for(int v=t;v!=s;v=prevv[v]){
            edge &e=G[prevv[v]][preve[v]];
            e.cap-=d;
            G[v][e.rev].cap+=d;
        }
    }
    return ;
}
 
int main()
{
    scanf("%d%d",&m,&n);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            scanf("%d",&T[j][i]);
        }
    }
    for(int i=1;i<=m;i++){
        for(int j=1;j<=n;j++){
            for(int k=1;k<=n;k++){
                add_edge((i-1)*n+j,(n*m+k),1,j*T[i][k]);
            }
        }
    }
    for(int i=1;i<=(n*m);i++)add_edge(0,i,1,0);
    for(int i=(n*m+1);i<=(n*m+n);i++)add_edge(i,(n*m)+n+1,1,0);
    min_cost_flow(0,(n*m)+n+1);
    double ans=res*1.0/n;
    printf("%.2f\n",ans);
    return 0;
}

有一种优化方法：如果第j个人没修倒数第i辆车，他一定不会修倒数第i+1、i+2...辆车，所以初始的图左边只需要有m个点，表示每个人的倒数第1辆车，当图中左边点(j,i)到右边点k流了1时，再将左边的(j,i+1)和右边的k连上边（权值和之前说的一样），这样总点数是2+m+n+n，不超过131，总边数是2∗(m+m∗n+n)+2∗(n+2)∗n，不超过10000，EK党的胜利！

#include<bits/stdc++.h>
using namespace std;
const int maxn=65;
const int maxm=15;
const int maxe=1e4+5;
const int inf=0x3f3f3f3f;
struct edge{int u,v,w,c,next;}e[maxe];
int head[150],tot=1;
void add_edge(int u,int v,int w,int c)
{
    e[++tot].u=u;e[tot].v=v;e[tot].w=w;e[tot].c=c;e[tot].next=head[u];head[u]=tot;
    e[++tot].u=v;e[tot].v=u;e[tot].w=0;e[tot].c=-c;e[tot].next=head[v];head[v]=tot;
}
int m,n,t[maxn][maxm],ans;
int q[150],level[150],pth[150];
bool inq[150];
int belong[150],fin[maxm],use;
bool spfa()
{
    int s=0,t=0;
    memset(level,inf,sizeof(level));
    level[0]=0;
    q[t++]=0;inq[0]=true;
    while(s!=t){
        int u=q[s++%150];
        inq[u]=false;
        for(int i=head[u];i;i=e[i].next){
            if(e[i].w&&level[e[i].v]>level[u]+e[i].c){
                level[e[i].v]=level[u]+e[i].c;
                pth[e[i].v]=i;
                if(!inq[e[i].v]){
                    q[t++%150]=e[i].v;
                    inq[e[i].v]=true;
                }
            }
        }
    }
    return level[n+1]<inf;
}
int Edmond_Karp()
{
    for(int i=pth[n+1];i;i=pth[e[i].u]){
        --e[i].w,++e[i^1].w;
        if(e[i].u==0)use=e[i].v;
    }
    return level[n+1];
}
int main()
{
    scanf("%d%d",&m,&n);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++)scanf("%d",&t[i][j]);
    }
    for(int i=1;i<=n;i++)add_edge(i,n+1,1,0);
    int ptot=n+1;
    for(int i=1;i<=m;i++){
        belong[++ptot]=i,++fin[i];
        for(int j=1;j<=n;j++)add_edge(ptot,j,1,t[j][i]);
    }
    for(int i=n+2;i<=ptot;i++)add_edge(0,i,1,0);
    while(spfa()){
        ans+=Edmond_Karp();
        int k=belong[++ptot]=belong[use];++fin[k];
        for(int i=1;i<=n;i++)add_edge(ptot,i,1,t[i][k]*fin[k]);
        add_edge(0,ptot,1,0);
    }
    printf("%.2f\n",(double)ans/n);
    return 0;
}