【分块答案】【最小割】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 100301 int 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; }
轉載于:https://www.cnblogs.com/autsky-jadek/p/4175677.html
總結
以上是生活随笔為你收集整理的【分块答案】【最小割】bzoj1532 [POI2005]Kos-Dicing的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: HTML5表单的创建及与PHP的交互
- 下一篇: JNI/NDK开发指南(三)——JNI数