网络流/费用流
跟 BZOJ 1070 修车 几乎是一道题，只是这题“要修的车”（即菜）多了很多……几乎是从$n$变成了$n^2$，所以建图的时候就得动态加点……
也就是说，当一个厨师已经确定了他的后k道菜的时候，再加入表示倒数第k+1道菜的节点。
我个xx一开始写成了每次加一层……不 T 才怪……
加点过程看min_cost_flow（即mcf）这个函数吧，每次通过找流过来的节点，来确定要加的点。
/**************************************************************
Problem: 2879
User: Tunix
Language: C++
Result: Accepted
Time:7416 ms
Memory:72872 kb
****************************************************************/ //BZOJ 2879
#include<cmath>
#include<vector>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
#define CC(a,b) memset(a,b,sizeof(a))
using namespace std;
int getint(){
int v=,sign=; char ch=getchar();
while(!isdigit(ch)) {if(ch=='-') sign=-; ch=getchar();}
while(isdigit(ch)) {v=v*+ch-''; ch=getchar();}
return v*sign;
}
const int N=,M=,INF=~0u>>;
const double eps=1e-;
/*******************template********************/
int n,m,ans,tot,tim[][],p[];
struct edge{int from,to,v,c;};
struct Net{
edge E[M];
int head[N],next[M],cnt;
void ins(int x,int y,int z,int c){
E[++cnt]=(edge){x,y,z,c};
next[cnt]=head[x]; head[x]=cnt;
}
void add(int x,int y,int z,int c){
ins(x,y,z,c); ins(y,x,,-c);
}
int from[N],Q[M],d[N],S,T,ed;
bool inq[N],sign;
bool spfa(){
int l=,r=-;
F(i,,T) d[i]=INF;
d[S]=; Q[++r]=S; inq[S]=;
while(l<=r){
int x=Q[l++];
inq[x]=;
for(int i=head[x];i;i=next[i])
if(E[i].v> && d[x]+E[i].c<d[E[i].to]){
d[E[i].to]=d[x]+E[i].c;
from[E[i].to]=i;
if (!inq[E[i].to]){
Q[++r]=E[i].to;
inq[E[i].to]=;
}
}
}
return d[T]!=INF;
}
void mcf(){
int x=INF,a,b;
for(int i=from[T];i;i=from[E[i].from]){
x=min(x,E[i].v);
if (E[i].to==T){
a=E[i].from%tot+;
b=(E[i].from-)/tot;
}
}
for(int i=from[T];i;i=from[E[i].from]){
E[i].v-=x;
E[i^].v+=x;
}
ans+=x*d[T];
if (a!=)
F(i,,n)
add(i,b*tot+a,,a*tim[i][b]);
}
void init(){
n=getint(); m=getint(); cnt=ed=; tot=;
F(i,,n) {p[i]=getint(); tot+=p[i];}
S=; T=tot+tot*m+;
F(i,,n) F(j,,m) tim[i][j]=getint();
F(i,,n) add(S,i,p[i],);
F(i,tot+,tot+tot*m) add(i,T,,);
F(i,,n)
F(j,,m)
add(i,j*tot+,,tim[i][j]);
ans=;
while(spfa()) mcf();
printf("%d\n",ans);
}
}G1;
int main(){
#ifndef ONLINE_JUDGE
freopen("2879.in","r",stdin);
freopen("2879.out","w",stdout);
#endif
G1.init();
return ;
}