CZ市为了欢迎全国各地的同学，特地举办了一场盛大的美食节。作为一个喜欢尝鲜的美食客，小M自然不愿意错过这场盛宴。他很快就尝遍了美食节所有的美食。然而，尝鲜的欲望是难以满足的。尽管所有的菜品都很可口，厨师做菜的速度也很快，小M仍然觉得自己桌上没有已经摆在别人餐桌上的美食是一件无法忍受的事情。于是小M开始研究起了做菜顺序的问题，即安排一个做菜的顺序使得同学们的等待时间最短。小M发现，美食节共有n种不同的菜品。每次点餐，每个同学可以选择其中的一个菜品。总共有m个厨师来制作这些菜品。当所有的同学点餐结束后，菜品的制作任务就会分配给每个厨师。然后每个厨师就会同时开始做菜。厨师们会按照要求的顺序进行制作，并且每次只能制作一人份。此外，小M还发现了另一件有意思的事情: 虽然这m个厨师都会制作全部的n种菜品，但对于同一菜品，不同厨师的制作时间未必相同。他将菜品用1, 2, ..., n依次编号，厨师用1, 2, ..., m依次编号，将第j个厨师制作第i种菜品的时间记为 ti,j 。小M认为：每个同学的等待时间为所有厨师开始做菜起，到自己那份菜品完成为止的时间总长度。换句话说，如果一个同学点的菜是某个厨师做的第k道菜，则他的等待时间就是这个厨师制作前k道菜的时间之和。而总等待时间为所有同学的等待时间之和。现在，小M找到了所有同学的点菜信息: 有 pi 个同学点了第i种菜品（i=1, 2, ..., n）。他想知道的是最小的总等待时间是多少。
输入文件的第1行包含两个正整数n和m，表示菜品的种数和厨师的数量。 第2行包含n个正整数，其中第i个数为pi，表示点第i种菜品的人数。 接下来有n行，每行包含m个非负整数，这n行中的第i行的第j个数为ti,j，表示第j个厨师制作第i种菜品所需的时间。 输入文件中每行相邻的两个数之间均由一个空格隔开，行末均没有多余空格。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
int read()
{
int x=0,f=1; char ch=getchar();
while (ch<‘0‘ || ch>‘9‘) {if (ch==‘-‘) f=-1; ch=getchar();}
while (ch>=‘0‘ && ch<=‘9‘) {x=x*10+ch-‘0‘; ch=getchar();}
return x*f;
}
#define maxp 1000
#define maxn 100010
#define maxm 1000010
int n,m,S,T,ans,tot;int p[maxp],tt[maxp][maxp];
struct data{int from,next,to,cap,cost;}edge[maxm<<2];
int head[maxn],from[maxn],cnt=1;
void add(int u,int v,int w,int c)
{cnt++;edge[cnt].from=u;edge[cnt].to=v;
edge[cnt].next=head[u];head[u]=cnt;
edge[cnt].cost=c;edge[cnt].cap=w;}
void insert(int u,int v,int w,int c){add(u,v,w,c);add(v,u,0,-c);}
#define inf 0x7fffffff
bool visit[maxn];int dis[maxn];
bool spfa()
{
queue<int>q;
for (int i=S; i<=T; i++) dis[i]=inf;
q.push(S); visit[S]=1; dis[S]=0;
while (!q.empty())
{
int now=q.front(); q.pop();
for (int i=head[now]; i; i=edge[i].next)
if (edge[i].cap && dis[now]+edge[i].cost<dis[edge[i].to])
{
dis[edge[i].to]=dis[now]+edge[i].cost; from[edge[i].to]=i;
if (!visit[edge[i].to])
q.push(edge[i].to),visit[edge[i].to]=1;
}
visit[now]=0;
}
return dis[T]!=inf;
}
void MinCostMaxFlow()
{
int flow=inf,a,b,x;
for (int i=from[T]; i; i=from[edge[i].from])
{
flow=min(flow,edge[i].cap);
if (!edge[i].from) x=edge[i].to,a=(x-1)/tot+1,b=x%tot+1;
}
for (int i=from[T]; i; i=from[edge[i].from])
edge[i].cap-=flow,edge[i^1].cap+=flow,ans+=edge[i].cost*flow;
for (int i=1; i<=n; i++)
insert((a-1)*tot+b,m*tot+i,1,b*tt[i][a]);
}
void make()
{
S=0,T=100001;
for (int i=1; i<=tot*m; i++)
insert(S,i,1,0);
for (int i=1; i<=n; i++)
insert(m*tot+i,T,p[i],0);
for (int i=1; i<=m; i++)
for (int j=1; j<=n; j++)
insert((i-1)*tot+1,m*tot+j,1,tt[j][i]);
}
int main()
{
n=read(),m=read();
for (int i=1; i<=n; i++) p[i]=read(),tot+=p[i];
for (int i=1; i<=n; i++)
for (int j=1; j<=m; j++)
tt[i][j]=read();
make();
while (spfa()) MinCostMaxFlow();
printf("%d\n",ans);
return 0;
}
平常写费用流喜欢写zkw费用流，这题用基于spfa的费用流会更方便（于是YY了一下T了...于是找了个模板...）
