bzoj千题计划139：bzoj2229: [Zjoi2011]最小割

http://www.lydsy.com/JudgeOnline/problem.php?id=2229

最小割树介绍：http://blog.csdn.net/jyxjyx27/article/details/42750833

http://blog.csdn.net/miaomiao_ymxl/article/details/54931876

#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>

using namespace std;

#define N 200
#define M 7000
const int inf=2e9;

int n;

int tot;
int front
,nxt[M<<1],to[M<<1],val[M<<1],from[M<<1];
int lev
,num
;
int path
;
int cur
;

int src,decc;

int a
,tmp
;
bool use
;

int dis

;

void read(int &x)
{
x=0; char c=getchar();
while(!isdigit(c))  c=getchar();
while(isdigit(c)) { x=x*10+c-'0'; c=getchar();  }
}

void add(int u,int v,int w)
{
to[++tot]=v; nxt[tot]=front[u]; front[u]=tot; from[tot]=u; val[tot]=w;
//to[++tot]=u; nxt[tot]=front[v]; front[v]=tot; from[tot]=v; val[tot]=0;
}

bool bfs()
{
queue<int>q;
for(int i=1;i<=n;++i) lev[i]=n;
q.push(decc);
lev[decc]=0;
int now,t;
while(!q.empty())
{
now=q.front();
q.pop();
for(int i=front[now];i;i=nxt[i])
{
t=to[i];
if(lev[t]==n && val[i^1])
{
lev[t]=lev[now]+1;
q.push(t);
}
}
}
return lev[src]!=n;
}

int augment()
{
int now=decc,flow=inf;
int i;
while(now!=src)
{
i=path[now];
flow=min(flow,val[i]);
now=from[i];
}
now=decc;
while(now!=src)
{
i=path[now];
val[i]-=flow;
val[i^1]+=flow;
now=from[i];
}
return flow;
}

int isap()
{
int flow=0;
if(!bfs()) return 0;
memset(num,0,sizeof(num));
for(int i=1;i<=n;++i) num[lev[i]]++,cur[i]=front[i];
int now=src,t;
while(lev[src]<n)
{
if(now==decc)
{
flow+=augment();
now=src;
}
bool advanced=false;
for(int i=cur[now];i;i=nxt[i])
{
t=to[i];
if(lev[t]==lev[now]-1 && val[i])
{
advanced=true;
path[t]=i;
cur[now]=i;
now=t;
break;
}
}
if(!advanced)
{
int mi=n;
for(int i=front[now];i;i=nxt[i])
if(val[i]) mi=min(mi,lev[to[i]]);
if(!--num[lev[now]]) break;
num[lev[now]=mi+1]++;
cur[now]=front[now];
if(now!=src) now=from[path[now]];
}
}
return flow;
}

void dfs(int x)
{
use[x]=true;
for(int i=front[x];i;i=nxt[i])
if(!use[to[i]] && val[i]) dfs(to[i]);
}

void solve(int l,int r)
{
if(l>=r) return;
src=a[l]; decc=a[r];
for(int i=2;i<=tot;i+=2) val[i]=val[i+1]=val[i]+val[i+1]>>1;
int flow=isap();
//    printf("%d %d %d\n",l,r,flow);
memset(use,false,sizeof(use));
dfs(src);
for(int i=1;i<=n;++i)
if(use[i])
for(int j=1;j<=n;++j)
if(!use[j])
dis[i][j]=dis[j][i]=min(dis[i][j],flow);
int i=l,j=r;
for(int k=l;k<=r;++k)
if(use[a[k]]) tmp[i++]=a[k];
else tmp[j--]=a[k];
for(int k=l;k<=r;++k) a[k]=tmp[k];
solve(l,i-1);
solve(j+1,r);
}

int main()
{
//    freopen("mincuto.in","r",stdin);
//    freopen("mincuto.out","w",stdout);
int T;
read(T);
int m,q;
int u,v,w,x;
int ans;
while(T--)
{
tot=1;
memset(front,0,sizeof(front));
memset(dis,127,sizeof(dis));
read(n); read(m);
for(int i=1;i<=n;++i) a[i]=i;
while(m--)
{
read(u); read(v); read(w);
add(u,v,w);
add(v,u,w);
}
solve(1,n);
read(q);
while(q--)
{
ans=0;
read(x);
for(int i=1;i<=n;++i)
for(int j=i+1;j<=n;++j)
if(dis[i][j]<=x) ans++;
cout<<ans<<'\n';
}
cout<<'\n';
}
return 0;
}

2229: [Zjoi2011]最小割

Time Limit: 10 Sec Memory Limit: 259 MB
Submit: 2429 Solved: 863
[Submit][Status][Discuss]

Description

小白在图论课上学到了一个新的概念——最小割，下课后小白在笔记本上写下了如下这段话： “对于一个图，某个对图中结点的划分将图中所有结点分成两个部分，如果结点s,t不在同一个部分中，则称这个划分是关于s,t的割。 对于带权图来说，将所有顶点处在不同部分的边的权值相加所得到的值定义为这个割的容量，而s,t的最小割指的是在关于s,t的割中容量最小的割” 现给定一张无向图，小白有若干个形如“图中有多少对点它们的最小割的容量不超过x呢”的疑问，小蓝虽然很想回答这些问题，但小蓝最近忙着挖木块，于是作为仍然是小蓝的好友，你又有任务了。

Input

输入文件第一行有且只有一个正整数T，表示测试数据的组数。 对于每组测试数据， 第一行包含两个整数n,m，表示图的点数和边数。 下面m行，每行3个正整数u,v,c(1<=u,v<=n,0<=c<=106)，表示有一条权为c的无向边(u,v) 接下来一行，包含一个整数q，表示询问的个数 下面q行，每行一个整数x，其含义同题目描述。

Output

对于每组测试数据，输出应包括q行，第i行表示第i个问题的答案。对于点对(p,q)和(q,p)，只统计一次（见样例）。

两组测试数据之间用空行隔开。

Sample Input

1

5 0

1

0

Sample Output

10

【数据范围】

对于100%的数据 T<=10,n<=150,m<=3000,q<=30，x在32位有符号整数类型范围内。

图中两个点之间可能有多条边