BZOJ2229: [Zjoi2011]最小割 分治最小割

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2229

至于什么是分治最小割就懒得在这里写了，网上的讲解应该很多，总之做n次最大流就可以得出n个点两两之间的最小割了，计算答案即可

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue> 
using namespace std;
#define LL long long
const int N=205;
const int M=10005;
const LL inf=1e12;
int n,m,cnt=1,to[M],nxt[M],lj[N],pre[N];
LL w[M],sum[N][N],f[N][N];
void ins(int f,int t,LL p) {cnt++,to[cnt]=t,nxt[cnt]=lj[f],lj[f]=cnt,w[cnt]=p;}
void add(int f,int t,LL p) {ins(f,t,p),ins(t,f,p);}
int d[N],S,T,TT;
queue<int>Q;
bool bfs()
{
	memset(d,0,sizeof(d));
	d[S]=1;
	Q.push(S);
	while(!Q.empty())
	{
		int x=Q.front();
		Q.pop();
		for(int i=lj[x];i;i=nxt[i])
		if(w[i]&&!d[to[i]])
		{
			d[to[i]]=d[x]+1;
			Q.push(to[i]);
		}
	}
	if(d[T]) return true;
	return false;
}
LL dfs(int x,LL v)
{
	if(x==T||v==0) return v;
	LL ret=0;
	for(int i=pre[x];i;i=nxt[i])
	if(d[to[i]]==d[x]+1)
	{
		LL f=dfs(to[i],min(w[i],v));
		pre[x]=i;
		w[i]-=f,w[i^1]+=f;
		v-=f,ret+=f;
		if(v==0) break;
	}
	if(ret==0) d[x]=-1;
	return ret;
}
bool b[405][405];
int tim;
void solve(int x,int y,int now)
{
    S=x,T=y;
    int pos=0;
	LL ans=0;
    for(int i=1;i<=n;i++)
    for(int j=i+1;j<=n;j++)
    if(sum[i][j])
	{
		pos+=2;
		w[pos]=w[pos^1]=sum[i][j];
    }
	while(bfs())
	{
		for(int i=1;i<=n;i++) pre[i]=lj[i]; 
		ans+=dfs(S,inf);
	}
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	if(d[i]&&!d[j]||!d[i]&&d[j]) f[i][j]=min(f[i][j],ans);
	int xx=0,yy=0;
	int t1=++tim;
	for(int i=1;i<=n;i++)
	if(b[now][i]&&d[i])
	{
		b[t1][i]=true;
		if(i!=S) xx=i; 
	}
	int t2=++tim;
	for(int i=1;i<=n;i++)
	if(b[now][i]&&!d[i])
	{
		b[t2][i]=true;
		if(i!=T) yy=i;
	}
	if(xx) solve(xx,x,t1);
	if(yy) solve(yy,y,t2);
}
void clear()
{
	memset(f,127/3,sizeof(f));
	memset(sum,0,sizeof(sum));
	memset(lj,0,sizeof(lj));
	memset(b,0,sizeof(b));
	cnt=1;
	tim=0;
}
int main()
{
	scanf("%d",&TT);
	while(TT--)
	{
		clear();
		scanf("%d%d",&n,&m);
		int x,y;
		LL z;
		for(int i=1;i<=m;i++)
		{
			scanf("%d%d%lld",&x,&y,&z);
			sum[x][y]+=z;
			sum[y][x]+=z;
		}
		for(int i=1;i<=n;i++)
		for(int j=i+1;j<=n;j++)
		if(sum[i][j]) add(i,j,sum[i][j]);
		for(int i=1;i<=n;i++) b[0][i]=1;
		solve(1,n,0);
		int q;
		scanf("%d",&q);
		while(q--)
		{
			scanf("%d",&x);
			int ans=0;
			for(int i=1;i<=n;i++)
			for(int j=i+1;j<=n;j++)
			if(f[i][j]<=x) ans++;
			printf("%d\n",ans);
		}
		if(TT) puts("");
	}
}