POJ 3189 Steady Cow Assignment

不错的一道题目，最大流和二分多重匹配都可以做

但是网络流我用EK或者Relabel-To-Front 都TLE了，网上说SAP不会超时，今天把二分多重匹配和SAP都学习了下

 

用二分匹配时，枚举区间很巧妙，设low为枚举等级区间的左极限，high为枚举等级区间的右极限

当这个区间可以满足要求时，说明可以缩小区间范围，low++

若不能，说明需要增加区间范围,high++

因为[low,high]满足时，[low,high+1]必然满足，[low+1,high]满足时，[low,high]也必然满足

这样，枚举次数就是O(b)了

 

另外，当数组大点，且数组初始化次数较多时，时间代价也是很大的，写二分匹配时，link数组开成1200×1200的，时间多花了300多MS

 

代码：

SAP版:

#include<iostream> using namespace std; #define MAX 1200 #define inf 9999999 #define min(x,y) (x)<(y)?(x):(y) struct node { int x,y,c,op,next; }g[120000]; int adj[MAX],cur[MAX],dis[MAX],num[MAX],pre[MAX];//pre[i],边(pre[i],i)在g中的标号 int rank[1200][25],barn[1200]; int n,b,max_flow,S,T,e; void add(int x,int y,int c) { e++; g[e].x=x; g[e].y=y; g[e].c=c; g[e].op=e+1; g[e].next=adj[x]; adj[x]=e; e++; g[e].x=y; g[e].y=x; g[e].c=0; g[e].op=e-1; g[e].next=adj[y]; adj[y]=e; } void relabel(int k) { int mm=T,t; cur[k]=adj[k];//当前弧指向首条邻接边 for(t=adj[k];t!=0;t=g[t].next) { if(g[t].c>0)//余量大于0 { mm=min(mm,dis[g[t].y]); } } dis[k]=mm+1; } void change() { int i,minf=inf;//该次最小增广量 for(i=T;i!=S;i=g[pre[i]].x) { minf=min(minf,g[pre[i]].c); } for(i=T;i!=S;i=g[pre[i]].x) { g[pre[i]].c-=minf; g[g[pre[i]].op].c+=minf; } max_flow+=minf; } void sap() { int i; max_flow=0; for(i=S;i<=T;i++) { cur[i]=adj[i]; dis[i]=0; num[i]=0; } num[0]=T+1; i=S; while(dis[S]<T+1)//T+1,节点总数 { //找节点i可以增广的路径(i,y) for(;cur[i]!=0;cur[i]=g[cur[i]].next)//cur[i]，节点i的当前弧在g中的标号 { if(g[cur[i]].c>0&&dis[i]==dis[g[cur[i]].y]+1) break; } if(!cur[i]) { num[dis[i]]--;dis[i]将要重标记 if(num[dis[i]]==0)//gap优化 break; relabel(i); num[dis[i]]++;//dis[i]已重标记 if(i!=S) i=g[pre[i]].x; } else { pre[g[cur[i]].y]=cur[i];//cur[i]为当前增广路径上的那条弧在g中的标号,即边(i,y)对应标号cur[i] i=g[cur[i]].y; if(i==T) { change();//最大流增加 i=S;//又重新从源点找增广路径 } } } } void range() { int i,j,k,l,a; for(l=0;l<b;l++) { for(i=1;i+l<=b;i++) { e=0; for(a=S;a<=T;a++) adj[a]=0; for(j=1;j<=n;j++) { add(S,j,1); for(k=i;k<=i+l;k++) { add(j,n+rank[j][k],1); } } for(j=n+1;j<=n+b;j++) { add(j,T,barn[j-n]); } sap(); if(max_flow==n) { printf("%d/n",l+1); return; } } } } int main() { int i,j; cin>>n>>b; for(i=1;i<=n;i++) { for(j=1;j<=b;j++) { scanf("%d",&rank[i][j]); } } for(i=1;i<=b;i++) scanf("%d",&barn[i]); S=0; T=b+n+1; range(); return 0; } 

 

二分匹配版：

#include<iostream> using namespace std; int link[25][1200],u[25]; int rank[1005][25],barn[25],num[25]; int n,b,low,high; int dfs(int t) { for(int i=1;i<=b;i++) { if(rank[t][i]<=high&&low<=rank[t][i]&&!u[i])//第t头牛住的第i号房间在等级范围内，且这次i还没被访问 { u[i]++; if(num[i]<barn[i]) { link[i][++num[i]]=t; return 1; } for(int j=1;j<=barn[i];j++) { if(dfs(link[i][j])) { link[i][j]=t; return 1; } } } } return 0; } bool hungary() { int i,ans=0; memset(link,-1,sizeof(link)); memset(num,0,sizeof(num)); for(i=1;i<=n;i++) { memset(u,0,sizeof(u)); if(!dfs(i)) return false; } return true; } void range() { int temp,ans; low=high=1; ans=99999; while(low<=high&&high<=b) { temp=high-low+1; if(hungary()) { ans=ans<temp?ans:temp; low++; } else high++; } printf("%d/n",ans); } int main() { int i,j,t; cin>>n>>b; for(i=1;i<=n;i++) for(j=1;j<=b;j++) { scanf("%d",&t); rank[i][t]=j;//第i头牛住的第j个房间的排名 } for(i=1;i<=b;i++) scanf("%d",&barn[i]); range(); return 0; }