首页 > 代码库 > 【分块答案】【最小割】bzoj1532 [POI2005]Kos-Dicing
【分块答案】【最小割】bzoj1532 [POI2005]Kos-Dicing
引用zky的题解:http://blog.csdn.net/iamzky/article/details/39667859
每条S-T路径代表一次比赛的结果。最小割会尽量让一个人赢得最多。
因为二分总是写挂,所以写了分块答案,比暴力枚举好像快不了多少。
#include<cstdio>#include<cstring>#include<algorithm>#include<queue>#include<cmath>using namespace std;#define INF 2147483647#define MAXN 20011#define MAXM 100301int v[MAXM],cap[MAXM],en,first[MAXN],next[MAXM];int d[MAXN],cur[MAXN];queue<int>q;int n,m,S,T;void Init_Dinic(){memset(first,-1,sizeof(first)); en=0; S=0; T=n+m+1;}void AddEdge(const int &U,const int &V,const int &W){v[en]=V; cap[en]=W; next[en]=first[U]; first[U]=en++;v[en]=U; next[en]=first[V]; first[V]=en++;}bool bfs(){ memset(d,-1,sizeof(d)); q.push(S); d[S]=0; while(!q.empty()) { int U=q.front(); q.pop(); for(int i=first[U];i!=-1;i=next[i]) if(d[v[i]]==-1 && cap[i]) { d[v[i]]=d[U]+1; q.push(v[i]); } } return d[T]!=-1;}int dfs(int U,int a){ if(U==T || !a) return a; int Flow=0,f; for(int &i=cur[U];i!=-1;i=next[i]) if(d[U]+1==d[v[i]] && (f=dfs(v[i],min(a,cap[i])))) { cap[i]-=f; cap[i^1]+=f; Flow+=f; a-=f; if(!a) break; } if(!Flow) d[U]=-1; return Flow;}int max_flow(){ int tmp=0,Flow=0; while(bfs()) { memcpy(cur,first,(n+m+5)*sizeof(int)); while(tmp=dfs(S,INF)) Flow+=tmp; } return Flow;}int us[10001],vs[10001];void Rebuild(const int &x){ Init_Dinic(); for(int i=1;i<=m;++i) { AddEdge(S,i,1); AddEdge(i,us[i]+m,1); AddEdge(i,vs[i]+m,1); } for(int i=1;i<=n;++i) AddEdge(i+m,T,x);}int main(){ scanf("%d%d",&n,&m); for(int i=1;i<=m;++i) scanf("%d%d",&us[i],&vs[i]); int sz=sqrt(m); int last=0; for(int i=1;last<=m;i+=sz) { Rebuild(i); if(max_flow()>=m) { for(int j=last+1;j<=i;++j) { Rebuild(j); if(max_flow()>=m) { printf("%d\n",j); return 0; } } } last=i; } return 0;}
【分块答案】【最小割】bzoj1532 [POI2005]Kos-Dicing
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。