题目大意:
有n种菜,需求为p_i。有m名厨师,第j名厨师做第i种菜的时间为t_ij。各厨师独立做菜,最小化上菜等待时间之和。
简单题解:
费用流(跪翱犇)。
一名厨师做的最后一个菜的时间被算了一次,第二个两次,依次类推。使用这种方式来统计答案。
初始时,菜向厨师连边,表示最后一次做菜。
有n种菜,需求为p_i。有m名厨师,第j名厨师做第i种菜的时间为t_ij。各厨师独立做菜,最小化上菜等待时间之和。
简单题解:
费用流(跪翱犇)。
一名厨师做的最后一个菜的时间被算了一次,第二个两次,依次类推。使用这种方式来统计答案。
初始时,菜向厨师连边,表示最后一次做菜。
对当前增广的厨师,新建一个点表示其前一次做菜,此时费用应为两倍;再次新建即为三倍,依次类推。
注意last数组可以省很多的时间,即记录第一个没有满流的地点,形成一个凸包。
#include <cstdio>
#include <queue>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define maxm 3000005
#define maxn 100010
#define INF 0x3f3f3f3f
struct node
{
int u,v,f,c,next;
}e[maxm];
int k,head[maxn],pre[maxn],dist[maxn],vis[maxn],last[maxn];
int next[maxn],cnt[maxn],to[maxn],cost[maxn],flow[maxn];
int en=1,s,t,maxflow,mincost,m,n,num,sum;
int save[500][500],need[500];
void add(int u,int v,int f,int c)
{
en++;
e[en].u=u;
e[en].v=v;
e[en].c=c;
e[en].f=f;
e[en].next=head[u];
head[u]=en;
en++;
e[en].u=v;
e[en].v=u;
e[en].c=-c;
e[en].f=0;
e[en].next=head[v];
head[v]=en;
}
int q[maxm];
int spfa()
{
int i,u,v;
for(i=0;i<=t;i++)
dist[i]=INF;
dist[s]=0;
vis[s]=1;
int H=0,T=0;
T++;q[T]=s;
while(H<T)
{
H++;
u=q[H];
for(i=head[u];i;i=e[i].next)
{
v=e[i].v;
if(e[i].f>0&&dist[u]+e[i].c<dist[v])
{
dist[v]=dist[u]+e[i].c;
pre[v]=i;
if(!vis[v])
{
vis[v]=1;
T++;
q[T]=v;
}
}
}
vis[u]=0;
}
}
void add2()
{
int v,x,y,a,b;
int maxf=INF;
for(v=pre[t];v;v=pre[e[v].u])
{
e[v].f--;
e[v^1].f++;
}
mincost+=dist[t];
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%d",&need[i]);
sum+=need[i];
}
for(int i=1;i<=n;i++)
add(s,i,need[i],0);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
scanf("%d",&save[i][j]);
add(i,n+j,1,save[i][j]);
}
}
t=n+m+sum+1;
for(int i=1;i<=m;i++)
{
cnt[i]=1;
add(n+i,t,1,0);
last[i]=en-1;
}
for(int j;sum;sum--)
{
spfa();
add2();
for(j=1;j<=m && e[last[j]].f;j++);
++cnt[j];
for(int i=1;i<=n;++i)
add(i,n+m+sum,1,save[i][j]*cnt[j]);
add(n+m+sum,t,1,0);
last[j]=en-1;
}
printf("%d\n",mincost);
return 0;
}
/*
3 10 2 3 3 2
5
6
7
*/