BZOJ 2229 ZJOI2011 最小割

Problem

BZOJ

Solution

这是一道可恨的结论题。结论是最小割不能相互分割，证明方法我也不造。
一共n-1个最小割，我们考虑分治，对于每次的割，会将点分为S集和T集，那么我们可以就此进行分治。
可能是我自带常数大的buff，洛谷上差点还T了？

Code

#include <cstring>
#include <cstdio>
#include <queue>
#define rg register
using namespace std;
const int maxn=160,maxm=3010,INF=0x3f3f3f3f;
struct data{int v,nxt,w;}edge[maxm<<1];
int z,n,m,q,p,dfn,a[maxn],t[maxn],vis[maxn],map[maxn][maxn];
int head[maxn],cur[maxn],dis[maxn];
queue<int> que;
struct max_flow{
    int bfs(int s,int t)
    {
        while(!que.empty()) que.pop();
        memset(dis,0xff,sizeof(dis));
        dis[s]=0;que.push(s);
        while(!que.empty())
        {
            int x=que.front();que.pop();
            for(int i=head[x];~i;i=edge[i].nxt)
              if(edge[i].w&&dis[edge[i].v]==-1)
              {
                dis[edge[i].v]=dis[x]+1;
                que.push(edge[i].v);
              }
        }
        return ~dis[t];
    }
    int dfs(int u,int v,int flow)
    {
        if(u==v) return flow;
        int rest=0,tmp;
        for(int &i=cur[u];~i;i=edge[i].nxt)
          if(dis[u]+1==dis[edge[i].v]&&edge[i].w)
          {
            tmp=dfs(edge[i].v,v,min(edge[i].w,flow-rest));
            rest+=tmp;edge[i].w-=tmp;edge[i^1].w+=tmp;
          }
        if(!rest) dis[u]=-1;
        return rest;
    }
    int dinic(int s,int t)
    {
        int res=0,tmp;
        while(bfs(s,t))
        {
            memmove(cur,head,sizeof(cur));
            while(tmp=dfs(s,t,INF)) res+=tmp;
        }
        return res;
    }
}g;
inline int min(int x,int y){return x<y?x:y;}
inline void insert(int u,int v,int w)
{
    edge[++p]=(data){v,head[u],w};head[u]=p;
    edge[++p]=(data){u,head[v],w};head[v]=p;
}
void input()
{
    scanf("%d%d",&n,&m);p=-1;
    memset(head,0xff,sizeof(head));
    memset(map,0x3f,sizeof(map));
    for(int i=1,u,v,w;i<=m;i++)
    {
        scanf("%d%d%d",&u,&v,&w);
        insert(u,v,w);
    }
}
void dfs(int x)
{
    vis[x]=dfn;
    for(int i=head[x];~i;i=edge[i].nxt)
      if(edge[i].w&&vis[edge[i].v]!=dfn)
        dfs(edge[i].v);
}
void binary(int l,int r)
{
    if(l>=r) return ;
    for(int i=0;i<=p;i+=2) edge[i].w=edge[i^1].w=(edge[i].w+edge[i^1].w)>>1;
    int cut=g.dinic(a[l],a[r]),L=l,R=r;
    ++dfn;dfs(a[l]);
    for(int i=l;i<=r;i++)
    {
        if(vis[a[i]]==dfn) t[L++]=a[i];
        else t[R--]=a[i];
    }
    for(int i=l;i<=r;i++) a[i]=t[i];
    for(int i=1;i<=n;i++)
      if(vis[i]==dfn)
        for(int j=1;j<=n;j++)
          if(vis[j]!=dfn) map[i][j]=map[j][i]=min(map[i][j],cut);
    binary(l,L-1);binary(R+1,r);
}
int main()
{
    #ifndef ONLINE_JUDGE
    freopen("in.txt","r",stdin);
    #endif
    scanf("%d",&z);
    while(z--)
    {
        input();
        for(int i=1;i<=n;i++) a[i]=i;
        binary(1,n);
        scanf("%d",&q);
        while(q--)
        {
            int x,ans=0;
            scanf("%d",&x);
            for(int i=1;i<=n;i++)
              for(int j=1;j<i;j++)
                if(map[i][j]<=x)
                  ans++;
            printf("%d\n",ans);
        }
        putchar('\n');
    }
    return 0;
}