[SCOI2007]修车【排队总时间最小化问题 费用流】

题目描述：

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

题目分析：

要求总等待时间最小。
一个技术人员依次修理各个顾客的车，如果按顺序表示修理这个车对总时间花费的贡献，则没法很好地解决后效性的问题。然后我们发现，如果用车$i$是在第$j$个技术人员手中倒数第$k$个修理的车来表示它对时间的贡献，则可以表示为$k\ast t[i][j]$。

于是把每个技术人员拆成$n$个点，表示在他的手中倒数第几个修理。然后每个顾客向每个技术人员的每个点连边，费用为$k\ast$维修时间。源点向顾客连边，容量为1，费用为0，技术人员的每个点向汇点连边，容量为1，费用为0,。

每个技术人员的n个点表示倒数第几个修理，那么可能过程中会有不连续的情况，但是这是肯定不优的（因为倒数越前花费越小），所以最后结果一定是连续的，即可行的。

Code：

#include<cstdio>
#include<deque>
#include<cstring>
#include<algorithm>
#define maxn 605
#define maxm 70005
using namespace std;
const int inf = 0x3f3f3f3f;
int n,m,S,T;
int fir[maxn],cur[maxn],nxt[maxm],to[maxm],c[maxm],w[maxm],tot=1;
inline void line(int x,int y,int z,int v){
	nxt[++tot]=fir[x],fir[x]=tot,to[tot]=y,c[tot]=z,w[tot]=v;
	nxt[++tot]=fir[y],fir[y]=tot,to[tot]=x,c[tot]=0,w[tot]=-v;
}
namespace Costflow{
	int dis[maxn],ans;
	bool vis[maxn],inq[maxn];
	deque<int>q;
	bool SPFA(){
		memset(dis,0x3f,(T+1)<<2);
		dis[T]=0,q.push_back(T);
		while(!q.empty()){
			int u=q.front();q.pop_front(),inq[u]=0;
			for(int i=fir[u],v;i;i=nxt[i]) if(c[i^1]&&dis[v=to[i]]>dis[u]+w[i^1]){
				dis[v]=dis[u]+w[i^1];
				if(!inq[v]) {inq[v]=1;if(!q.empty()&&dis[v]<dis[q.front()]) q.push_front(v);else q.push_back(v);}
			}
		}
		return dis[S]!=inf;
	}
	int aug(int u,int augco){
		if(u==T) {ans+=dis[S]*augco;return augco;}
		vis[u]=1;
		int need=augco,delta;
		for(int &i=cur[u];i;i=nxt[i]) if(c[i]&&!vis[to[i]]&&dis[u]==dis[to[i]]+w[i]){
			delta=aug(to[i],min(need,c[i]));
			c[i]-=delta,c[i^1]+=delta;
			if(!(need-=delta)) break;
		}
		vis[u]=0;
		return augco-need;
	}
	int solve(){
		ans=0;
		while(SPFA()) memcpy(cur,fir,(T+1)<<2),aug(S,inf);
		return ans;
	}
}
int main()
{
	scanf("%d%d",&m,&n),S=0,T=n+m*n+1;
	for(int i=1,x;i<=n;i++) for(int j=1;j<=m;j++){
		scanf("%d",&x);
		for(int k=1;k<=n;k++) line(i,n+(j-1)*n+k,1,k*x);
	}
	for(int i=1;i<=n;i++) line(S,i,1,0);
	for(int i=1;i<=m*n;i++) line(n+i,T,1,0);
	printf("%.2f\n",Costflow::solve()*1.0/n);
}