有N个订单和M个机器,给出第i个订单在第j个机器完成的时间Mij,每台机器同一时刻只能处理一个订单,机器必须完整地完成一个订单后才能接着完成下一个订单。问N个订单完成时间的平均值最少为多少。
把每个工厂分成n个,每个代表某时刻,某时刻的工厂只能执行一个。这样解释就很合理了。
但是费用怎么求呢。
对于某个工厂,T=(z1) + (z1+z2)+ (z1+z2+z3)
变形得到 T=1*z3 + 2*z2 + 3*z1 !!!
也就是倒数第i个进场,他的花费就是i*z
所以对于某个工厂,某个玩具在不同时刻的花费就确定了。z,z*2,z*3,z*k类推。
有人可能无法理解,觉着这样如果两个玩具都选择了同一时刻的工厂怎么办。
这是行不通过的,每一时刻的工厂还连接着汇点,流量是1,这代表着只能走一个,所以不会出现问题
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
using namespace std;
const int N=4200;
const int MAXE=535000;
const int inf=1<<30;
int head[N],s,t,cnt,n,m;
int flow=0,ans=0;
int d[N];
int pre[N];
bool vis[N];
int q[MAXE+79000];
struct Edge
{
int u,v,next;
int c,w;
}edge[MAXE];
void addedge(int u,int v,int c,int w)
{
edge[cnt].u=u;
edge[cnt].v=v;
edge[cnt].w=w;
edge[cnt].c=c;
edge[cnt].next=head[u];
head[u]=cnt++;
edge[cnt].v=u;
edge[cnt].u=v;
edge[cnt].w=-w;
edge[cnt].c=0;
edge[cnt].next=head[v];
head[v]=cnt++;
}
int SPFA()
{
for(int i=0;i<=t;i++)
d[i]=inf;
memset(pre,-1,sizeof(pre));
memset(vis,0,sizeof(vis));
d[s]=0;
int l=0,r=0;
q[r++]=s;
vis[s]=1;
while(l<=r)
{
int u=q[l++];
vis[u]=0;
for(int j=head[u];j!=-1;j=edge[j].next)
{
int v=edge[j].v;
if(edge[j].c>0&&d[u]+edge[j].w<d[v])
{
d[v]=d[u]+edge[j].w;
pre[v]=j;
if(!vis[v])
{
vis[v]=1;
q[r++]=v;
}
}
}
}
if(d[t]==inf)
return 0;
return 1;
}
void MCMF()
{
flow=0;
while(SPFA())
{
ans+=d[t];
int u=t;
int mini=inf;
while(u!=s)
{
if(edge[pre[u]].c<mini)
mini=edge[pre[u]].c;
u=edge[pre[u]].u;
}
flow+=mini;
u=t;
while(u!=s)
{
edge[pre[u]].c-=mini;
edge[pre[u]^1].c+=mini;
u=edge[pre[u]].u;
}
}
}
int mp[160][160];
int main()
{
int cas;
scanf("%d",&cas);
while(cas--)
{
scanf("%d%d",&n,&m);
s=0;
ans=0;
flow=0;
t=n+n*m+1;
memset(head,-1,sizeof(head));
cnt=0;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&mp[i][j]);
for(int i=1;i<=n;i++) //s到每一个玩具,都制作一个
addedge(s,i,1,0);
for(int i=1;i<=n;i++)//对于当前玩具
for(int j=1;j<=m;j++)//对于当前工厂
for(int k=1;k<=n;k++)//当前工厂的每一时刻
{
addedge(i,n+(j-1)*n+k,1,mp[i][j]*k); //话费为v*k
}
for(int i=1;i<=m;i++)//工厂
for(int j=1;j<=n;j++)//每一时刻
{
addedge(n+(i-1)*n+j,t,1,0);//每一工厂的每一时刻都与汇点链接
}
MCMF();
printf("%.6lf\n",ans*1.0/n);
}
return 0;
}