BZOJ 2561: 最小生成树 最小割

2561: 最小生成树

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 2025  Solved: 967
[Submit][Status][Discuss]

Description

　给定一个边带正权的连通无向图G=(V,E)，其中N=|V|，M=|E|，N个点从1到N依次编号，给定三个正整数u，v，和L (u≠v)，假设现在加入一条边权为L的边(u,v)，那么需要删掉最少多少条边，才能够使得这条边既可能出现在最小生成树上，也可能出现在最大生成树上？

Input

　　第一行包含用空格隔开的两个整数，分别为N和M；
　　接下来M行，每行包含三个正整数u，v和w表示图G存在一条边权为w的边(u,v)。
　　最后一行包含用空格隔开的三个整数，分别为u，v，和 L；
　　数据保证图中没有自环。

Output

　输出一行一个整数表示最少需要删掉的边的数量。

Sample Input

3 2
3 2 1
1 2 3
1 2 2

Sample Output

1

HINT

对于20%的数据满足N ≤ 10，M ≤ 20，L ≤ 20；
对于50%的数据满足N ≤ 300，M ≤ 3000，L ≤ 200；
对于100%的数据满足N ≤ 20000，M ≤ 200000，L ≤ 20000。

就是看删掉多少边权小于/大于L的边能使得U,V不连通

最小割上

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>
using namespace std;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=20100,M=200100,inf=0X3f3f3f3f;

int n,S,T;

int ecnt=1,last[N];
struct EDGE{int to,nt,val;}e[M<<1];
inline void readd(int u,int v,int val)
{e[++ecnt]=(EDGE){v,last[u],val};last[u]=ecnt;}
inline void add(int u,int v)
{readd(u,v,1);readd(v,u,1);}

int flow;

int d[N],q[N];

bool bfs()
{
	memset(d,0,sizeof(d));
	register int head=0,tail=1,u,i;
	q[head]=S;d[S]=1;
	while(head<tail)
	{
		u=q[head++];
		for(i=last[u];i;i=e[i].nt)if(e[i].val&&!d[e[i].to])
		{d[e[i].to]=d[u]+1;q[tail++]=e[i].to;}
	}
	return d[T];
}

int dfs(int u,int lim)
{
	if(u==T||!lim)return lim;
	register int i,res=0,tmp;
	for(i=last[u];i;i=e[i].nt)
	if(e[i].val&&d[u]+1==d[e[i].to])
	{
		tmp=dfs(e[i].to,min(e[i].val,lim));
		lim-=tmp;res+=tmp;e[i].val-=tmp;e[i^1].val+=tmp;
		if(!tmp)d[e[i].to]=-1;if(!lim)break;
	}
	return res;
}

void dinic()
{while(bfs())flow+=dfs(S,inf);}

int U[M],V[M],VAL[M];

int main()
{
	n=read();int m=read();
	register int i,val;
	for(i=1;i<=m;++i){U[i]=read();V[i]=read();VAL[i]=read();}
	S=read();T=read();val=read();
	for(i=1;i<=m;++i)if(VAL[i]<val)add(U[i],V[i]);
	dinic();
	ecnt=1;memset(last,0,sizeof(last));
	for(i=1;i<=m;++i)if(VAL[i]>val)add(U[i],V[i]);
	dinic();
	print(flow);puts("");
	return 0;
}
/*
3 2
3 2 1
1 2 3
1 2 2

1
*/