BZOJ 1070 修车 费用流(拆点)

题意:同一时刻有N位车主到维修中心,告诉你a[i][j]第i个师傅修第j个人的车需要的时间.

现在需要安排这M位师傅所维修的车及顺序,使得顾客平均等待的时间最小.

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

2<=M<=9,1<=N<=60

左边定点i为师傅 右边点j为车主 花费为a[i][j] 跑流量为n的最小费用最大流.

wawawawaawa,一个师傅修多辆车时 等待时间需要累加!!

比如某师傅修了2,1,4 则2,1,4总共等待时间为3*a[2]+2*a[1]+a[4].怎么表示累加时间?

拆点:把每个师傅拆成n个点,每个点都与源点连流量1,费用0的边.
每个师傅拆出来的第i个点和每辆车相连,费用为i*a[x][y],代表:这是该师傅修的倒数第i俩车.(费用得到累加.)

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> ii;
const int N=2e3+20,inf=0x3f3f3f3f;
double n;
int m,s,t;
int cnt,head
;
struct node{
int u,to,next,flow,cost;
}e[N*N];
void add(int u,int v,int flow,int cost)
{
e[cnt].u=u,e[cnt].to=v,e[cnt].flow=flow,e[cnt].cost=cost;
e[cnt].next=head[u];
head[u]=cnt++;
}
void Add(int u,int v,int flow,int cost)
{
add(u,v,flow,cost);
add(v,u,0,-cost);
}
int pre
,inq
,f
,c
;
int SPFA(int s,int t)
{
memset(pre,-1,sizeof(pre));
queue<int> q;
for(int i=0;i<=t;i++)
c[i]=inf,inq[i]=0,f[i]=0;
q.push(s),c[s]=0,inq[s]=1,f[s]=inf;
while(!q.empty())
{
int u=q.front();
q.pop(),inq[u]=0;
for(int i=head[u];i!=-1;i=e[i].next)
{
int v=e[i].to,w=e[i].cost;
if(e[i].flow>0&&c[v]>c[u]+w)
{
c[v]=c[u]+w;
f[v]=min(f[u],e[i].flow);
pre[v]=i;
if(!inq[v])
inq[v]=1,q.push(v);
}
}
}
return f[t];
}
int a

;
int MCMF(int s,int t)
{
int ans=0,i,tmp;
while(tmp=SPFA(s,t))
{
for(int i=pre[t];i!=-1;i=pre[e[i].u])
{
int res=tmp*e[i].cost;
ans+=tmp*e[i].cost;
e[i].flow-=tmp;
e[i^1].flow+=tmp;
}
}
return ans;

}
int main()
{
int x;
cin>>m&g
4000
t;>n;
s=0,t=m*n+n+1;
memset(head,-1,sizeof(head));
cnt=0;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&a[i][j]);

for(int k=1;k<=n;k++)
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
Add((k-1)*m+j,m*n+i,1,k*a[i][j]);
for(int j=1;j<=n*m;j++)
Add(s,j,1,0);
for(int i=1;i<=n;i++)
Add(m*n+i,t,1,0);
double ans=MCMF(s,t);
ans=ans/n;
printf("%.2lf\n",ans);
return 0;
}