POJ 3189 二分+Dinic
题意:
思路:
二分跨度 枚举最低座次
建图:源点向每头牛连边权为1的边 每头牛向当前枚举的B的区间这段连上边权为1的边 所有座次向汇点连边权为牛棚容量的边
判判流量是不是等于n
一开始写得是直接枚举答案,在这个答案下枚举座次最低的值 T了。。
(然后我看了一发数据 )
数据中有这样一个点:1000 20 balabala
答案是20
这就很蛋疼了 我枚举了200次 略多 (虽然在我电脑上还是能1s之内跑出来 但是POJ评测机并不快啊)
我就把答案的那个20二分了一下 变成了6
枚举大概几十次 嗯 A了 938ms
//By SiriusRen
#include <queue>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 1055
int n,b,map[N][25],B[25],e=1050,answer;
struct Dinic{
int first[N],next[60*N],v[60*N],w[60*N],vis[N],tot;
void add(int x,int y,int z){w[tot]=z,v[tot]=y,next[tot]=first[x],first[x]=tot++;}
bool tell(){
memset(vis,-1,sizeof(vis));
queue<int>q;q.push(0);vis[0]=0;
while(!q.empty()){
int t=q.front();q.pop();
for(int i=first[t];~i;i=next[i])
if(w[i]&&vis[v[i]]==-1)
q.push(v[i]),vis[v[i]]=vis[t]+1;
}
return vis[e]!=-1;
}
int zeng(int x,int y){
if(x==e)return y;
int r=0;
for(int i=first[x];~i&&y>r;i=next[i])
if(vis[v[i]]==vis[x]+1&&w[i]){
int t=zeng(v[i],min(w[i],y-r));
w[i]-=t,w[i^1]+=t,r+=t;
}
if(!r)vis[x]=-1;
return r;
}
int flow(){
int ans=0,jy;
while(tell())while(jy=zeng(0,0x3fffffff))ans+=jy;
return ans;
}
bool solve(int begin,int end){
memset(first,-1,sizeof(first)),tot=0;
for(int i=1;i<=n;i++)
for(int j=begin;j<=end;j++)
add(i,1000+map[i][j],1),add(1000+map[i][j],i,0);
for(int i=1;i<=b;i++)add(1000+i,e,B[i]),add(e,1000+i,0);
for(int i=1;i<=n;i++)add(0,i,1),add(i,0,0);
return flow()==n;
}
}dinic;
int main(){
scanf("%d%d",&n,&b);
for(int i=1;i<=n;i++)
for(int j=1;j<=b;j++)
scanf("%d",&map[i][j]);
for(int i=1;i<=b;i++)scanf("%d",&B[i]);
int l=0,r=b;
while(l<=r){
int Mid=(l+r)>>1;
for(int i=1;i+Mid<=b;i++)
if(dinic.solve(i,i+Mid)){
answer=Mid,r=Mid-1;
goto ed;
}
l=Mid+1;
ed:;
}
printf("%d\n",answer+1);
}