BZOJ 2879 NOI2012 美食节 费用流

题目大意：给定n道菜和m个厨师，第i道菜需要p[i]份，第j个厨师做第i道菜需要时间t[i][j]，求最长总等待时间

一个厨师做的倒数第一道菜对答案的贡献是时间的一倍，倒数第二道菜对答案的贡献是时间的两倍，以此类推

厨师们怒了！发动符卡·禁忌『p重存在』！

将每个厨师拆成Σp[i]个点，每道菜向每个厨师的第i个点连一条流量为1，费用为时间的i倍，每个点向汇点连一条流量为1费用为0的边，跑最小费用最大流

这能跑？不会T到死？

动态加点即可

每个厨师初始只有一个点，如果这个厨师做完了最后一道菜，就添加一个点，这样可以保证厨师的点数是O(p+m)的

妈的T到死- - 一条增广路上只有一个厨师是做菜的- - 其它都是退菜的- - 因此要找到第一个厨师进行加点- - 这SB错误卡了一下午- -

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 100100
#define S 0
#define T 100099
#define INF 0x3f3f3f3f
using namespace std;
struct abcd{
	int to,flow,cost,next;
}table[1001001];
int head[M],tot=1;
int n,m,now,ans,a[50][110];
int belong[M],cnt[M],cost[M];
void Add(int x,int y,int z,int w)
{
	table[++tot].to=y;
	table[tot].flow=z;
	table[tot].cost=w;
	table[tot].next=head[x];
	head[x]=tot;
}
void Link(int x,int y,int z,int w)
{
	Add(x,y,z,w);
	Add(y,x,0,-w);
}
bool Edmonds_Karp()
{
	static int q[65540],flow[M],from[M];
	static unsigned short r,h;
	static bool v[M];
	int i,j;
	memset(cost,0x3f,sizeof cost);
	cost[S]=0;q[++r]=S;flow[S]=INF;flow[T]=0;//cost[T]=INF;
	while(r!=h)
	{
		int x=q[++h];v[x]=0;
		for(i=head[x];i;i=table[i].next)
			if( table[i].flow && cost[table[i].to]>cost[x]+table[i].cost )
			{
				cost[table[i].to]=cost[x]+table[i].cost;
				flow[table[i].to]=min(flow[x],table[i].flow);
				from[table[i].to]=i;
				if(!v[table[i].to])
					v[table[i].to]=1,q[++r]=table[i].to;
			}
	}
	if(!flow[T]) return false;
	ans+=flow[T]*cost[T];
	int x=0;
	for(i=from[T];i;i=from[table[i^1].to])
	{
		table[i].flow-=flow[T];
		table[i^1].flow+=flow[T];
		if(table[i].to<100000)
			x=belong[table[i].to];
	}
	Link(S,++now,1,0);
	for(j=1;j<=n;j++)
		Link(now,100000+j,1,a[j][x]*(cnt[x]+1));
	cnt[x]++;
	belong[now]=x;
	return true;
}
int main()
{

	#ifndef ONLINE_JUDGE
	freopen("2879.in","r",stdin);
	freopen("2879.out","w",stdout);
	#endif

	int i,j,x;
	cin>>n>>m;
	for(i=1;i<=n;i++)
	{
		scanf("%d",&x);
		Link(100000+i,T,x,0);
	}
	for(i=1;i<=n;i++)
		for(j=1;j<=m;j++)
			scanf("%d",&a[i][j]);
	for(j=1;j<=m;j++)
	{
		Link(S,++now,1,0);
		for(i=1;i<=n;i++)
			Link(now,100000+i,1,a[i][j]);
		cnt[j]++;
		belong[now]=j;
	}
	memset(cost,0x3f,sizeof cost);
	while( Edmonds_Karp() );
	cout<<ans<<endl;
	return 0;
}