BZOJ 1532: [POI2005]Kos-Dicing|网络流

二分答案网络流

#include<set>
#include<map>
#include<ctime>
#include<queue>
#include<cmath>
#include<cstdio>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define inf 1e9
using namespace std;
int sc()
{
	int i=0; char c=getchar();
	while(c>'9'||c<'0')c=getchar();
	while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
	return i;
}
int dis[20010],q[20010];
int head[20010],lst[98888],nxt[98888],c[98888];
int u[10001],v[10001];
int n,m,S,W,tot;
void insert(int x,int y,int z)
{
	lst[++tot]=y; nxt[tot]=head[x]; c[tot]=z; head[x]=tot;
	lst[++tot]=x; nxt[tot]=head[y]; c[tot]=0; head[y]=tot;
}
bool BFS()
{
	for(int i=1;i<=W;i++)dis[i]=0; dis[S]=1;
	int l=1,r=2;q[1]=S;
	while(l<r)
	{
		int x=q[l++];
		for(int i=head[x];i;i=nxt[i])
		    if(c[i]&&dis[lst[i]]==0)
		    {
		    	dis[lst[i]]=dis[x]+1;
		    	q[r++]=lst[i];
		    }
	}
	return dis[W];
}
int dfs(int x,int f)
{
	if(x==W)return f;
	int ww=0,w;
	for(int i=head[x];i;i=nxt[i])
	    if(c[i]&&dis[lst[i]]==dis[x]+1)
	    {
	    	w=dfs(lst[i],min(c[i],f-ww));
	    	ww+=w; c[i]-=w,c[i^1]+=w;
	    	if(ww==f)return f;
	    }
	if(ww==0)dis[x]=0;
	return ww;
}		
bool jud(int x)
{
	int ans=0; 
	memset(head,0,sizeof(head));	tot=1;
	for(int i=1;i<=n;i++) insert(S,i,x);
	for(int i=1;i<=m;i++) insert(u[i],i+n,1),insert(v[i],i+n,1),insert(i+n,W,1);
	while(BFS())ans+=dfs(S,inf);
	//cout << x << " " << ans << endl;
	return ans==m;
}	
int main()
{
	n=sc(),m=sc();
	S=n+m+1,W=S+1;
	for(int i=1;i<=m;i++) u[i]=sc(),v[i]=sc();
	int l=1,r=m,ans=r;
	while(l<=r)
	{
		int mid=l+r>>1;
		if(jud(mid))ans=mid,r=mid-1;
		else l=mid+1;
	}
	cout<<ans;
	return 0;
}


你可能感兴趣的:(网络流)