码迷,mamicode.com
首页 > 其他好文 > 详细

【分块答案】【最小割】bzoj1532 [POI2005]Kos-Dicing

时间:2014-12-20 18:07:58      阅读:112      评论:0      收藏:0      [点我收藏+]

标签:

引用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;
}

  

【分块答案】【最小割】bzoj1532 [POI2005]Kos-Dicing

标签:

原文地址:http://www.cnblogs.com/autsky-jadek/p/4175677.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!