【bzoj 2521】 [Shoi2010] 最小生成树（网络流最小割）

2521: [Shoi2010]最小生成树

Time Limit: 10 Sec   Memory Limit: 128 MB
Submit: 429   Solved: 250
[ Submit][ Status][ Discuss]

Description

Secsa最近对最小生成树问题特别感兴趣。他已经知道如果要去求出一个n个点、m条边的无向图的最小生成树有一个Krustal算法和另一个Prim的算法。另外，他还知道，某一个图可能有多种不同的最小生成树。例如，下面图 3中所示的都是图 2中的无向图的最小生成树：

当然啦，这些都不是今天需要你解决的问题。Secsa想知道对于某一条无向图中的边AB，至少需要多少代价可以保证AB边在这个无向图的最小生成树中。为了使得AB边一定在最小生成树中，你可以对这个无向图进行操作，一次单独的操作是指：先选择一条图中的边 P1P2，再把图中除了这条边以外的边，每一条的权值都减少1。如图 4所示就是一次这样的操作：

Input

输入文件的第一行有3个正整数n、m、Lab分别表示无向图中的点数、边数、必须要在最小生成树中出现的AB边的标号。

接下来m行依次描述标号为1,2,3…m的无向边，每行描述一条边。每个描述包含3个整数x、y、d，表示这条边连接着标号为x、y的点，且这条边的权值为d。

输入文件保证1<=x,y<=N,x不等于y,且输入数据保证这个无向图一定是一个连通图。

Output

输出文件只有一行，这行只有一个整数，即，使得标号为Lab边一定出现最小生成树中的最少操作次数。

Sample Input

4 6 1
1 2 2
1 3 2
1 4 3
2 3 2
2 4 4
3 4 5

Sample Output

1

HINT

第1个样例就是问题描述中的例子。

1<=n<=500,1<=M<=800,1<=D<10^6

Source

day2

[ Submit][ Status][ Discuss]

【题解】【网络流最小割】

【由最小生成树（kruskal）可知，若一条边必出现在一个图的最小生成树中，当且仅当用小于当前边的边不可使这条边控制的两个点相连】

【由这一性质可得网络流建图方法：将原图中所有不大于指定边的权值的边的权值改为val[id]-val[j]+1建入图中，双向连边（由于不确定行走方向），然后以指定边连接的两个点作为源点和汇点跑最小割，并输出答案即可】

#include
#include
#include
#include
using namespace std;
struct map{
	int s,t,v;
}d[810];
int a[4010],nxt[4010],p[510],val[4010],tot;
int cur[510],dis[510],last[510],num[510];
int n,m,id,ans;
bool b[510];
inline void add(int x,int y,int v)
{
	tot++; a[tot]=y; nxt[tot]=p[x]; p[x]=tot; val[tot]=v;
	tot++; a[tot]=x; nxt[tot]=p[y]; p[y]=tot; val[tot]=0;
}
inline int  addflow(int s,int t)
{
	int s1=0x7fffffff,now=t;
	while(now!=s) 
	 {
	 	s1=min(s1,val[last[now]]);
	 	now=a[last[now]^1];
	 }
	now=t;
	while(now!=s)
	 {
	 	val[last[now]^1]+=s1;
	 	val[last[now]]-=s1;
	 	now=a[last[now]^1];
	 }
	return s1;
}
inline void bfs(int l)
{
    queueque;
	for(int i=1;i<=n;++i) dis[i]=n;
	memset(b,0,sizeof(b));
	que.push(l); dis[l]=0; b[l]=true;
	while(!que.empty())
	 {
	 	int u,v;
	 	u=que.front(); que.pop();
	 	v=p[u];
	 	while(v!=-1)
	 	 {
	 	 	if(!b[a[v]]&&val[v^1])
	 	 	 {
	 	 	 	b[a[v]]=1;
	 	 	 	que.push(a[v]);
	 	 	 	dis[a[v]]=dis[u]+1;
			   }
			v=nxt[v];
		  }
	 }
	return;
}
inline int isap(int s,int t)
{
	int sum=0,now=s,i;
	bfs(t);
	for(i=1;i<=n;++i) num[dis[i]]++;
	for(i=1;i<=n;++i) cur[i]=p[i];
	while(dis[s]