poj 3189

很久没有一次a的感觉啦。关键是建模，建一个有b+2个点的网络流，s为0，t为b+1，1--b为各个棚。二分枚举答案。s到各棚的容量为当前可行区间内，可选择此棚的牛的个数，各棚到t的容量为割棚可容纳牛的个数。数据量不大，直接用的EK.

#include <iostream>
#include <cstdio>
#include <queue>
using namespace std;
const int maxn=1001;
const int maxe=25;
const int inf=200000000;
int n,b,tot,s,t;
int flow[maxe][maxe],cap[maxe][maxe],ra[maxn][maxe],c[maxe];
void build(int l,int r)
{
	for(int i=1;i<=b;i++) cap[i][t]=c[i];
	for(int i=l;i<=r;i++)
	{
		for(int j=1;j<=n;j++)
			cap[s][ra[j][i]]+=1;
	}
}
bool maxflow()
{
	queue<int> q;
	int d[maxe],p[maxe],tflow=0,a[maxe];
	for(;;)
	{
		memset(d,-1,sizeof(d));
		d[s]=0;a[s]=inf;
		q.push(s);
		while(!q.empty())
		{
			int x=q.front();q.pop();
			for(int i=s;i<=t;i++)
			{
				if(d[i]==-1&&cap[x][i]>flow[x][i])
				{
					a[i]=min(a[x],cap[x][i]-flow[x][i]);
					d[i]=d[x]+1;
					p[i]=x;
					q.push(i);
				}
			}
		}
		if(d[t]==-1) break;
		tflow+=a[t];
		int u=t;
		while(u!=s)
		{
			flow[p[u]][u]+=a[t];
			flow[u][p[u]]-=a[t];
			u=p[u];
		}
	}
	if(tflow<n) return false;
	else return true;
}
int solve()
{
	int high=b+1,low=1,mid;
	while(high-low>0)
	{
		bool flag=false;
		mid=low+(high-low)/2;
		for(int i=1;i<=b-mid+1;i++)
		{
			memset(flow,0,sizeof(flow));
			memset(cap,0,sizeof(cap));
			build(i,i+mid-1);
			if(maxflow())
			{
				flag=true; 
				break;
			}
		}
		if(flag) high=mid;
		else low=mid+1;
	}
	return high;
}
int main()
{
	while(cin>>n>>b)
	{
		s=0;t=b+1;
		int i,j;
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=b;j++)
			{
				scanf("%d",&ra[i][j]);
			}
		}
		for(i=1;i<=b;i++) scanf("%d",&c[i]);
		cout<<solve()<<endl;
	}
	return 0;
}